FINAL REPORT 


(NflSi-Cfi- 120580) SESEAECH STUDY OK 875-150^ 

STABILIZATION AND CONTSOL; KOBE BN 
SAMPLED-DATA CONTROL THEORY. DESIGN OF TBE 
LABGE SPACE TELESCOPE SYSTEM Final fieport Unclas 

v.(jJ s i: e M s Research Lab. , Champaign, 111.) G3/37 Q 8 Q 55 > 


RESEARCH STUDY ON STABILIZATION AND CONTROL 
MODERN SAMPLED-DATA CONTROL THEORY 


SYSTEMS RESEARCH LABORATOR 


P.O. BOX 2277, STATION A 
3206 VALLEY BROOK DRIVE 
CHAMPAIGN, ILLINOIS 61620 



PREPARED FOR GEORGE C. MARSHALL SPACE FLIGHT CENTER 


HUNTSVILLE, ALABAMA 



1-75 


F I NA L_ REPORT 

RESEARCH STUDY ON STABILIZATION AND CONTROL 
— MODERN SAMPLED -DATA CONTROL THEORY 


SUBTITLE: 


DESIGN OF THE 
LARGE SPACE TELESCOPE 
SYSTEM 


January 1, 1975 NAS8-29853 


BY B.C. KUO 

G. SINGH 


PREPARED FOR GEORGE C. MARSHALL SPACE FLIGHT CENTER 

HUNTSVILLE. ALABAMA 
CONTRACT NAS8 - 29853 DCN 1 2-40 23018 


SYSTEMS RESEARCH LABORATORY 

P.O.BOX 2 277, STATION A 
CHAMPAIGN. ILLINOIS 61820 



11 




TABLE OF CONTENTS 


Page 

1. Prediction by Numerical Methods of Self- Sustained Oscillations 

in a Two-Axis Model of the Nonlinear Continuous-Data LST System 1 

1-1. Introduction 

1-2. The Continuous-Data Two-Axis LST System Model and Its 
Stability Equation 

1-3. Prediction of Self-Sustained Oscillations by the 
Approximation Method 

1-4. A Direct Method of Predicting Self-Sustained Oscillations 
in the Two-Axis LST System 

1- 5. Exact Solution of the Stability Equation by Numerical- 

Iterative Techniques 

2. Prediction by Numerical Methods of Self-Sustained Oscillations 

in a Two-Axis Model of the Nonlinear Sampled-Data LST System 62 

2- 1. Introduction 

2-2. The Sampled-Data Two-Axis LST System Model and Its 
Stability Equation 

2- 3. Exact Solution of the Stability Equation of the Two- 

Axis Sampled-Data LST System by Numerical-Iterative 
Techniques 

3. Study of Unequal -Ampl itude Oscillations in Two-Axis Coupled 

Nonlinear Systems 93 

3- 1. Introduction 

3 - 2 . Conditions of Unequal -Ampl itude Oscillations in the 
Continuous-Data Two-Axis Coupled LST System 

3-3. A Simple Two-Axis Coupled Nonlinear System 

3-4. A Simplified Multiple-Loop Two-Axis Coupled Nonlinear 
Sys tern 



ABSTRACT 


The objective of this investigation is to study the conditions 
of self-sustained oscillations in a two-axis model of the nonlinear 
LST system. 

The describing function of the CMG frictional nonlinearity of 
the LST system is used for the analysis, and both the continuous-data 
and the discrete-data models of the simplified LST control system are 
used. 

When two nonlinear systems are coupled together, the conventional 
graphical method using describing function can no longer by used for 
the prediction of limit cycles. A numerical-iterative method is 
described and used for the analysis of the two-axis system. Approxi- 
mation methods as well as the direct plotting of the stability equation 
are also used in the study. 

It is shown that although the dynamics of the two axes are 
identical, the amplitudes of self-sustained oscillations in the two 
axes may in principle be different. Analysis shows that for the LST 
system, most of the LST systems are of equal amplitudes but with 180- 
degree phase shift. 

The techniques described in this report can be extended to non- 
linear systems with more than two axes without too much difficulty. 
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1. Prediction by Numerical Methods of Self-Sustained Oscillations in 
a Two-Axis Model of the Nonlinear Continuous-Data LST System 

1-1 . Introduction 

It has been demonstrated in [1] that the methods of continuous 
and discrete describing function analysis can be applied to predict 
the existence of self-sustained oscillations in the single-axis model 
of the LST system with nonlinear CMG friction characteristics. 

Furthermore, it has been shown in [2] that the stability equations 
as a result of the describing function analysis may be solved by a 
numerical-iterative technique instead of the usual graphical methods. 

With an appropriate guess of the initial condition, the numerical 

method is found to be quite effective in leading to a convergent solution 

rapidly. 

In the two-axis model of the LST system, the system contains 
two nonlinearities, and the general form of the stability equation in 
the continuous-data case is 

1 + Gg( jw)N(A) + G c (jw)N 2 (A) = 0 (1-1) 

where Gg(jui) and G^(jto) denote transfer functions which depend on the 
linear elements of the coupled systems, and N(A) represents the describing 

A A 

function of the CMG nonlinearity. In general, Gg(ju)) and G c (jw} are 
functions of frequency at, and N(A) is a function of the amplitude A of 
the assumed sinusoidal input to the nonl inearity. 

The usual graphical method cannot be used to solve for the values 



2 


of to and A for self-sustained oscillations in Eq. (1-1), due to the 
2 

N (A) term in the equation. However, if the system parameters are 

A 

such that the linear term Gg(ju>)N(A) dominates over the quadratic 

term, or vice versa, then by neglecting the smaller quantity in Eq. (1-1), 

an approximate solution may still be obtained graphically. 

A direct but more time consuming approach of solving Eq. (1-1) would 
be to calculate the terms in Eq. (1-1) for a wide range of values of w 
and A until all combinations which satisfy the equation are found. 

Still a third alternative of solving Eq. (1-1) is to use the 
numerical-iterative method reported in [2]. If an appropriate initial 
guess can be made, the numerical -iterative method should lead to the 
exact solutions in an efficient manner. 

All three of the above-mentioned methods have been applied to the 
two-axis continuous-data LST system. It is found that the graphical 
method with approximation and the direct method can provide useful 
information to the final solutions, and the numerical-iterative method 
is effective in arriving at the exact solution on a digital computer. 

The results of these studies are reported in the ensueing sections. 

1-2. The Continuous-Data Two-Axis LST System Model and Its Stability 
Equation 

Figure 1-1 Shows the signal flow graph representation of the 
continuous-data single-axis LST system. The nonlinear friction of the 
CMG is modeled by the branch with the gain N. If two such models are 
coupled together at the output stage through output torsional coupling, 
the two-axis model of Fig. 1-2 results. This representation of a 
continuous-data two axis LST system may not be a rigorous one from the 




Figure 1-1. The simplified single-axis continuous-data 
LST system model . 




Figure 1-2. The simplified two-axis continuous-data LST system model 
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structural standpoint. However, the purpose of the present study is to 
develop analytical techniques which are applicable to multi-coupled 
nonlinear systems. It is conjectured that for a two-axis nonlinear 
system, the describing function analysis will generally lead to a 
stability equation of the form of Eq. (1-1). Therefore, for these 
above-mentioned reasons, the model shown in Fig. 1-2 is considered to 
be adequate. The techniques developed in this report can be applied 
to the prediction of self-sustained oscillations in any continuous- 
data two-axis nonlinear system that is amenable to the describing 
function method. 

With reference to Fig. 1-2, it can be seen that there are seven 
individual loops, with the following loop gains: 

NG 0 


= - -r (’-3) 


B 1 

= - g 2 g 3 

0-4) 

B 2 

= - g 2 g 3 

0-5) 

C 1 

= " G 1 G 2 G 3 G 4 H 

0-6) 

C 2 

= - 

0-7) 

D = 

2 2 

0-8) 


where G-j , G 2 , G 3 , and G^ are transfer functions as defined in Fig. 1-2. 
All system parameters and variables, including the nonl inear describing 
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function N, are consistent with those defined in the single-axis LST 
model reported previously [1,2]. The describing function N{A) is a 
function of the amplitude of the input sinusoid of the CMG nonlinearities. 
In the present case, since the two axes are identical, and the couplings 
are symmetric, it is assumed that the amplitudes of the input signals 
at the two nonlinearities are identical. 

The characteristic equation of the coupled system in Fig. 1-2 is 

A = 0 (1-9) 


where 


A = 1 - (A-| + ^2 + B] + + C] + C 2 + D) 

+ {A^A^ + A-jBg + + A-jD + B-|A 2 + ^1^2 + ^1^2 + B 1 B 

+ C-|A 2 + B2 + C 1 C 2 + A^D + BgD) 

- (A 1 A 2 D + A 1 B 2 D + B^D + B^D) (1-10) 

The last equation is simplified if we define 


A » 

A 1 

“ A 2 





B = 

B 1 

= B 2 



[V 

-11) 

C = 

C 1 

C 2 





Eq. 

(1 

-10) t 

becomes 




A = 

1 

- 2 (A 

+ B + C) 

- D + (A 2 + B 2 + C 2 + 2AB + 2AC + 2BC 





+ 2AD 

+ 2BD) - 

(A + B)V 

(1- 

-12) 


or 



7 


A = 1 - 2(A + B + C) - D + (A + B + C) 2 - (A + B) 2 D 


Substituting Eqs. (1-2) through (1-8) in Eq. (1-13) yields 


1 + -T- + 2G 2 G 3 + 2G 1 G 2 G 3 G 4 H ‘ K 2 G 4 


G 3 2 2 2 2 2NG 9 G 9 

+ 4 n + + (GjGgGgG^r + -—- 1 


2NG 1 G 2 G 2 G 4 H 2 , 2NG 3 G 4 K 2 

+ 5 + 2G 1 G 2 G 3 G 4 H 5 


„ ccc 2 r 2 N ' G 3 G 4 K 2 2NG 2 G 3 G 4 K 2 

“ 2G 2 G 3 G 4 K 2 T ? 

r 2 r 2 r 2 2 
G 2 G 3 G 4 K 2 


The characteristic equation in Eq. (1-9) can be written as 


A - 1 + G^ + GgN + G^N = 0 


2„2 , „ 2„2 


G A 2G 2 G 3 + 2G 1 G 2 G 3 G 4^ K 2 G^ + G 2 G^ + ( ) 


+ 2G 1 G 2 G 2 G 4 H - 2G 2 G 3 G 2 K 2 - G 2 G 2 G 4 K 2 


!G 3 2G..G 2 2G,G 2 g2g 4 H 

* 1 1 1 t * - 

S S 5 

2G 2 G 2 G 2 K 2 


2G 3 6 4 K 2 


0-17) 
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G„ = -4 - 


g ! g I k 1 


(1-18) 


For stability analysis, the characteristic equation of Eq. (1-15) 
is written as 


1 + GgN + G q N = 0 


(1-19) 


where 


A 



g b 

1 + g a 


( 1 - 20 ) 


Equation (1-19) is a function of the system frequency w and the 
input amplitude A. Since this equation is defined in the complex 
plane, it represents a set of two nonlinear equations with two unknowns. 

When these nonlinear equations have a solution for w and A, it represents 
a condition of self-sustained oscillations for the system. The oscillations 
may be stable or unstable; therefore, the solutions must always be 
checked for stability. 

1-3. Prediction of Self-Sustained Oscillations By The Approximation 
Method 

It was mentioned in Sec. 1-1 that the graphical method of predicting 
self-sustained oscillations can still be applied to the two-axis 1ST 
system if one of the last two terms in Eq. (1-1) can be neglected. In 
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other words, the following two conditions may exist: 
1. |G B ( jco)N(A) | » |G c (j W )N 2 (A)i 
Then Eq. (1-1) may be approximated by 


1 + G B ( jto)N(A) = 0 (1-22) 

and the condition of self-sustained oscillations is found from 
the following equation: 

G B (jw) = ‘ NW (1 " 23) 

2. |G B (jo>)N(A)| « iG c ( J oj)N 2 (A) | 

Then Eq. (1-1) is approximated by 

1 + G c (jw)N 2 (A) = 0 (1-24) 

and the condition of self-sustained oscillations is found from 

G r (jw) = - 0-25) 

C N 2 (A) 

Therefore, the graphical solutions involve the plotting of the 

/\ f-. r\ 

curves for Gg(joj) and -1/N(A), or G^(jto) and-l/fr(A), as the case 
may be. 

Figure 1-3 illustrates the plots of |Gg(jw)N(A) | and |G C ( jto)N 2 (A) | 
versus u for various values of A, and |N(A)j versus A. The following 
parameters are used for the LST system: 

J v = 10 5 , J Q = 2.1 , Kp = 216 , Kj = 9700 



DECIBELS 



Figu 
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100 
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H = 600 , K q = 5758.35 , K ] = 1371.02 , l< 2 = 100 , K 3 = 3 

Y = 1.38 x 10 7 

The last two parameters, K 2 and are the coefficients of coupling 
between the two axes. 

The curves in Fig. 1-3 give information on the ranges of to and 
A in which the approximation of Eq. (1-22) or Eq. (1-24) is valid. 

Figure 1-4 illustrates the regions in the to versus A plane in which the 
two approximations are valid. The cross-hatched area represents the 
region in which no approximation can be made, and the graphical method 
cannot be used. The criterion of >20 db and <-20 db is used for 
magnitude comparison for significance. 

The results of Fig. 1-4 show that the graphical approximation 
method is valid for the following ranges of w and A: 

1 . Use 1 + G C N 2 = 0 : 

A <_ 10” 7 w < 40 rad/sec 

and 60 < to _< 80 rad/sec 

2. Use 1 + GgN = 0 : 

A > 6 x 10~ 5 to >_ 600 rad/sec. 

For the region of validity of case 1 above, Fig. 1-5 shows the 
~ 2 

plots of Gq and -1/N- . The heavy portions of the curves indicate the 
parts which are valid for stability analysis. Similarly, Fig. 1-6 
shows the case when the equation 1 + G^N = 0 may be used for approxi- 

A 

mation. Again, the intersection between the heavy portions of the G D 

D 

and the -1/N curves would indicate the possibility of self-sustained 
oscillation in the two-axis LST system. Since there are no intersections 
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between the valid curves in Figs. 1-5 and 1-6, the system under study 
is stable for the parameter values used. This result will be sub- 
stantiated by the two other methods which are discussed in the following 
sections. 

1-4. A Direct Method of Predicting Self-Sustained Oscillations in 

The Two-Axis 1ST System 

A direct but time-consuming method of predicting self-sustained 
oscillations in the two-axis continuous-data LST system is to calculate 
all the terms in the following stability equation for a wide range of 
values of u» and A: 

1 + Gg{ joj)N(A) + G c (joj)N 2 (A) = 0 (1-26) 

Although the direct method would be quite tedious, and it is 
possible that with a given selected increment of variation for the 
values of to and A, the solution runs may still miss the exact solutions, 
the results in general will give insight to the appropriate guess on 
the initial values for the numerical-iteration method. Therefore, it 
is enlightening to investigate this direct approach at this point. 

Figures 1-7 through 1-9 show the plots of Gg(jco)N(A) + G c ( jco)N 2 (A) 
for y = 1.38 x 10 , and all the system parameters as given previously, 
and for three different sets of values for the coupling coefficients 
Kg and K^. Note that for the three cases illustrated, all the trajectories 
for various combinations of A and w do not intersect the critical point 
which is at 0 db and -180 degrees. Thus the system is stable. However, 
the results show that as the value of Kg is decreased, the trajectories 














DEGREES 


Figure 1-8. 






20 


-60 


-3< 


DEGREES 

Figure 1-9. 
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get closer to the critical point. Figure 1-10 shows the trajectories 
when l<2 = 0.1 and = 0.1. It is noticed that several trajectories 
are very close to the critical point. Figure 1-11 shows a magnified 
version of Fig. 1-10 around the critical point. The figure shows that 
the trajectories for A = 5 x 10 10 5 x 10'^, are all very close to 

the critical point. This gives indications that there may be more than 
one solution. It will be shown in the next section by the numerical- 
iterative method that there are indeed two solutions for the stability 
equation. One is at A = 5.9867 x 10 ^ and u> = 1.88 rad/sec, and the 
other is at A = 5.07397 x 10 ^ and w = 4.1086 rad/sec. These solutions 
must still be checked for stable or unstable equilibrium solutions. 

1 _ 5 . Exact Solution of the Stability Equation by Numerical-Iterative 

Techniques 

In this section the stability equation developed in Sec. 1-2 is 
solved numerically for its exact solutions. The numerical method 
utilized has been described in [2] and is found to be quite effective 
for the two-axis LST system. 

The stability equation, Eq. (1-19), can be written as 

1 + Og(jco)N(A) + Gg(ju)N^(A) = 0 (1-27) 

Define 

/\ 

G B (jw) = G r1 + jG n (1-28) 

A 

G c (jw) = G R2 + jG I2 (1-29) 

N(A) = N r1 + jN n 


(1-30) 





DEGREES 


Figure 1-11 
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N 2 (A) = N r2 + (1-31) 

Where G R -j * ^11 * ^R2> ^12* W r -| 5 ^11 * ^R2’ Nj 2 are dll real 
quantities. 

If Eqs. (1-28) through (1-31) are substituted into Eq. (1-27), 
it becomes 

1 + (G R1 + jG Il )(N Rl +jN Il ) + (G R2 + J' G I2 )(N R2 + jN I2 } = 0 0-32) 

When the real and imaginary parts are separated, Eq. (1-32) yields 

A = A r + jAj = 0 (1-33) 

with 

A R = 1 + G R1 N R1 " G I1 N I1 + G R2 N R2 " G I2 N I2 = 0 (1-34) 

and 

A I = G I1 N R1 + G R1 N I1 + G I2 N R2 + G R2 N I2 = 0 (1-35) 

Further simplification is possible in Eqs. (1-34) and (1-35) 
if we recognize that N R2 + jNj 2 = (N R1 + jN^ ) but is not necessary 
since the solution of Eqs. (1-34) and (1-35) will be performed on a 
digital computer. 

Equations (1-33) and (1-34) represent two equations in the two 
variables w and A. As before [2], let us define 

f w 1 


x 


A 


(1-36) 
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Then, Eqs. (1-34) and (1-35) can be written as 
F(x) = 0 


0 - 37 ) 


0 - 38 ) 


The Newton-type quadratically convergent numerical method 
described in [2] can now be directly applied to this two-variable 
system. 

To initiate the iterations, an initial solution is needed. The 
results of the direct calculation method, presented in the previous 
section provide an adequate initial solution. As in that section, the 
parameters which need to be varied are Kg and Kg. Once a solution for 
a given value of IC, and K 3 is known, Kg and/or Kg could be changed 
slightly and the new solution is obtained by using the old solution 
as the initial guess for the new system. 

With K 2 and K 3 set to very small values, the system is decoupled, 
and the solution of the single-axis LST system could be used as an 
initial guess for that of the two-axis LST system. 

If Kg = 0.1 and Kg = 0.1, the results from the direct calculation 
method shown in Figs. 1-10 and 1-11 indicate that two possible solutions 
exist. The approximate amplitudes and frequencies of the solutions are 

Solution No. 1 to = 5 rad/sec 

A = 5 x io' 6 rad 
to = 2 rad/sec 
A = 1 x 1 o~ 6 rad 


Solution No. 2 
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With these approximate solutions as initial conditions, the results 
of the numerical iteration procedure are obtained and tabulated in the 
first set of iterations in Fig. 1-12 and Fig. 1-13. These iterations 
indicate that the actual solutions are 

Solution No. 1 w = 4.1086 rad/sec 

A = 5.0739 x 10 " 6 rad 

Solution No. 2 w = 1.88 rad/ sec 

A = 5.99 x 10 -7 rad 

Figures 1-12 and 1-13 also show the change in the solutions when 
l <2 and Kg are reduced to - K 3 - 0.01 and = Kg = 0 . 001 , respectively. 

As mentioned previously, the initial guess for a new set of values of 
K 2 and Kg is the solution with the previous values of K 2 and Kg. 

With Kg fixed at 0.001, if K 2 is increased (Figs. 1-14 and 1-15), 
the solutions will change, until Kg = 8 for solution no. 1 and K 2 = 7 
for solution no. 2, beyond which no solution is obtained. It is 
interesting to note that the two solutions move closer together with 
increasing K 2 and almost merge into each other before disappearing 
altogether. The result is analogous to lowering the G ( j oj) curve or 
raising the -1/N curve in the single-axis case. 

It can be concluded that for Kg = 0.001, Kg should be less than 
8 for a sustained oscillation to occur. In fact, with K 2 = 8 (solution 
1) and K 2 = 7 (solution 2) the roots of the stability equations do not 
appear to have converged adequately, and well defined roots occur only 
for lower values of Kg. A maximum of 20 iterations are attempted before 
the numerical process is terminated. 

If Kg is kept fixed at 0.2 and K 2 is varied, the iterations in 
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Figs. 1-16 and 1-17 show that can be increased only up to = 3 
for solutions to be obtained. 

On the other hand, can be kept fixed and varied; these results 
are shown in Figs. 1-18 through i-25. In Figs. 1-18 and 1-19 K 2 is 
fixed at 0.001, and the maximum value of for a solution to the stability 
equation is approximately K 3 = 0.4 for solution no. 1 and K 2 = 0.3 for 
solution no. 2. With K 2 = 0.1, 1.0 and 5.0, the maximum values of are 
0.4 (solution 1 and 2), 0.3 (solution 1 and 2) and <0.1 (solution 1 and 2), 
respecti vely. These results are shown in Figs. 1-20 through 1-25. 

Table 1-1 summarizes all the results. In the table, the solutions 
for K 2 = ^3 = 0 are the solutions with the single-axis LST system. It is 
seen that when K 2 and K 3 are very small, the solutions of the two-axis 
system approach those of the single-axis system. 

The data in Table 1-1 can be used to plot the region in the K 2 - 
plane in which sustained oscillation can exist. Figure 1-26 shows this 
region. For values of K 2 and outside the crosshatched area, no 
sustained oscillations should exist. 

The motion of the roots with variable X 2 and are plotted in 
Fig. 1-27, using the data in Table 1-1. The points marked K 2 = 0, 

= 0, represent the single-axis LST system solutions. The solid closed 
regions represent the regions of root movements where reliable convergence 
was obtained. The dotted lines include the regions where reliable 
convergence of the iterative scheme was not obtained with 20 iterations, 
and, which can be considered as a boundary area between existence and non- 
existence of solutions. Note that in Fig. 1-27 the two solutions are 
farthest apart when K 2 = = 0. Figure 1-27 also shows how the two 



solutions tend to merge just before the solutions disappear. 



Table 1-1 



ROOT #1 


ROOT n 


FREQUENCY AMPLITUDE FREQUENCY AMPLITUDE 


5.99 x 


3.95** 


5.46** 


x IQ" 6 ** 


4.74 


3.81 x 10' u 
3.72 x lo" 6 
7.13 x ](f 6 ** 
4.17 x 10" 6 
3.64 x 10~ 6 

2.86 x 10" 6 

5.57 x 10" 6 
5.54 x 10" 6 


3.97** 


4.58 x 

4.45 x 

4.45 x 

4.46 x 

4.59 x 
5.79 x 
7.17 x 
8.69 x 
1.04 x 
1.25 x 
1.53 x 
2.06 x 


5.79 x 
5.93 x 
7.71 x 
1.03 x 
1.43 x 
2.57 x 
1.25 x 
1.30 x 


3.75** 2.25 x 

2.39 1 .04 x 

2.79 1.37 x 

3.28 1.95 x 

1.75 4.45 x 

1.76 4.56 x 


3.2** 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N-*N*G£=0 . 

GAMMA* 1 .3800011 07 TD= l.QOOQOD-Ol JV= 1.00000D 05 JG= 2.100000 00 
KP= 2.160000 02 KI= 9 .70 0 0 CIO '03 H= 6.000000 02 K0= 5.753350 03 
K 1 = 1.371020 03 K2= 1.000000-01 K3= 1.000000-01 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 .5000 0000 OOD 

01 

0 .5 000 0 000 0011-05 

1 

0.47337629720 

01 

0.5433 07 13490- 05 

£ 

0.45265635841) 

01 

0.5454 333 4 73D ~ 0 5 

3 

0.4375323004H 

01 

0 . 5 3 6 U 0 1 5 5 3 6D “ 0 5 

4 

0 . 427 097596213 

01 

0 .5259907533D-05 

5 

0 .42 032325730 

01 

0.51 83 1 87 33 2 D - 0 5 

6 

0 .4161 3 2*3 85 5 r 

01 

0.51 32 0 0 0856D- 05 

7 

0.4 13628645 ID 

01 

0.510019 083 1 D- 05 

3 

0 .41213130 OSD 

01 

0 .50806654730—05 

9 

0.41 1 157345 OD 

01 

0.50661877200-05 

10 

0.412522728311 

01 

0.5118 03 06 0 OD- 05 

11 

0 .41 15433437D 

01 

0 . 5 0994373980-05 

12 

0.41 04 4339 76 D 

01 

0 .50751567700-05 

13 

0.41 0335 1 055D 

01 

0 .50750903830-05 

14 

0.41 08 776 58 3D 

01 

0.50744499310-05 


ND= 1 FREQUENCY* 4.1086 00 OORRD-'SEC AMPLITUDE* 5.073970-06 NIT* 14 


LST SYSTEM-NUMERICAL SDLUTIDM DF 1+N*G1+N-*N*G2=0 . 

GAMMA* 1.380000 07 TO* 1.000000-01 JV= 1.000000 05 JG= £.100000 00 

KP= 2 . 16-0000 02 KI = 9.700000 03 H= 6.000000 02 K0= 5.75S35D 03 

Kl= 1.371 020 03 K2= 1.0 0 0 0 00-02 K3= 1.0000 00- 02 

ITERATION FREQUENCY AMPLITUDE 


0 

0.41 085972 03D 

01 

0 . 5 0739665840- 05 

1 

0.41688201450 

01 

0 . 524828886 OD- 05 

2 

0 .42055317850 

01 

u .536435353211-05 

3 

0. 4226335 066D 

01 

0 . 543225973? D- 05 

4 

0.4237936472H 

01 

0 - 547 066358 1 D- 05 

5 

0 .424450S200H 

01 

0 . 5 492 764 6 59 D— 05 

6 

0.4246721 233 D 

01 

0.5477545295D-05 

-y 

0 . 4253298522 D 

01 

0.5511 02375 OD- 05 

3 

0.48 5 3 8 1 8 8 6 8 D 

01 

0.551 754 0846D- 05 

9 

0.425357313211 

01 

0.55226 0793 OD- 05 

10 

0 .486869652 1 D 

01 

0 .552 053333 OD- 05 

11 

0. 4255358 006B 

01 

0 . 55 124774 52B- 05 

18 

0.425524647 ID 

01 

0 .551 63 3 9666 D— 05 


ND= 1 FREQUENCY* 4.25522D OORAD^SEC AMPLITUDE* 5.521140-06 NIT* 12 


LST SYSTEM-NUMERICAL SOLUTION DF 1 +N*G 1 +N*N*G2= 0 . 
GAMMA* 1.330000 07 T0= 1.000000-01 JV= 1.000000 05 

KP* 2.160000 02 K I = 9.700000 03 H= 6.000000 02 

Kl= 1.371020 03 K2- 1.000000-03 k'3= 1. 000000-03 


JG= 2.100000 00 
K 0= 5.753350 03 


ITERflTim 

FREQUENCY 


AMPLITUDE 

0 

0.4 25521 7824 D 

01 

0.5521 1 4S485D- 05 

1 

0 -426 067 6424 D 

01 

0 .554279531 3D- 05 

o 

L_ 

0. 42565348 16D 

01 

0 .5327396968D-G5 


0.42526 13 07 OD 

01 

0.548 067562211-05 

4 

0 .486 1346964 D 

01 

0. 549651 763 OD- 05 

5 

0 .42666934 5 9 D 

01 

0.55361 355 09B- 05 

6 

0 .426953 054 4D 

01 

0 .5557835647D-05 

7 

0 .42690423 1SB 

01 

0 .556 1545739D-05 


ND= 1 FREQUENCY* 4.269020 OORAL.-SEC AMPLITUDE* 5.566450-06 NIT* 7 


Figure 1-12 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1+N*N^G£=0. 

GAMMA® 1.330000 07 TO= 1 .00 000D- 01 .JV= 1.00 GOOD 05 JG= 2. 10 GOOD 00 

KP® 2.16 0 0 OD 02 KI= 3. 70 GOOD 03 H= 6.00000D 02 KO= 5.75835D 0-3 

K 1 = 1.371 02 D 03 K2= 1 . 0 00 0 OD- 0 1 K3= 1 . 0 0 0 0 OD- 0 1 

ITERATION FREQUENCY AMPLITUDE 


u 

0 . 2 

0 0 0 0 fi o 0 0 OD 

01 

0 . 1 0 0 0 0 0 0 0 0 0D- 05 

1 

0.1 

349233 024 D 

01 

0 . 62428486 U2D- 06 


0 . 1 

35295 0575 D 

01 

u - 5yyy98 117* OD— Ob 


0.1 

367 1 39237D 

01 

i_) * P 9 y 0 bc6c 1 4D- U b 

4 

0 . 1 

374829545D 

01 

U . 5 9 6 c 8 7 2 8 0 6 D — 0 6 

C T 

0 . 1 i 

37889 023 ID 

01 

0 . 5979 1 38 586 D- 06 

6 

0.1 

38 0 1 29 0 1 1 D 

01 

0 . 5979328 084D- 06 

-j 

( 

0 ■ 1 

38 Oc 0 1 346D 

01 

0 . 5986 1 539 1 OD- 06 

8 

0.1 

38033435 OD 

01 

0 . 59867 05996D- 06 


ND= 1 FREQUENCY® 1 .38Q42D OORfiD-'SEC AMPLITUDE® 5.98812D-G? HIT® 3 


LST SYSTEM-NUMERICAL SOLUTION OF l+N*bl+N*N*b£=0 . 

GAMMA® 1 .33 GOOD 07 TO= 1.00000D-01 JV= 1.00000D 05 JG= 2.10000D 00 

KP® 2.16 0 0 OD 02 KI= 9. 70 GOOD 03 H= 6. GO GOOD 02 KG® 5 . 75S35D 03 

Kl= 1.37102D 03 K2= 1.0 000 OH- 02 K3® 1.00000D-02 

ITERATION FREQUENCY . AMPLITUDE 


0 

0.1 

y y u 4c 1 y 0 r D 

01 

0.59831 13404D-06 

1 

0.1 

816575499D 

01 

0.52047436 0 8 D— 06 

C. 

0.1 

78779 111 6D 

01 

0 .4874 09 03 3 2D- 06 

•“» 

0 . 1- 

774 071 43 OD 

01 

0 . 4 72 0 024797D— 06 

4 

0.1 

76731 0255D 

01 

0 . 46453 3 3 4 9 C B “ 0 6 

5 

0.1 

763726929D 

01 

0.46 05878423D - 06 

b 

0.1 

75337421 7D 

01 

0 .45257934 09D- 06 

*7 

J 

0.1 

7583 12 03 OD 

01 

U .461 38776 USD- 06 

3 

0.1 

758882 1 66D 

01 

0.45 3 3 095 3 8 £ D - 0 6 

9 

0.1 

7596 0682 7D 

01 

0 . 4 5776. 1 £ 065D- 06 

10 

0.1 

759617287D 

01 

0 .4581 187707D-06 


ND= 1 FREQUENCY® 1 .75963P OORRD/'SEC AMPLITUDE® 4.58129D-07 NIT= 10 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N+61 +N*N^G2=0 . 

GAMMA® 1 .33 GOOD 07 TD= 1 . 00000D-01 JV= 1 . 0 00 GOD 05 JG = 2.10000D 00 

KP= 2.16000D 02 KI = 9.70000D 03 H= 6.00000D 02 K0= 5.75S35D 03 

Kl= 1.37102D 03 K£= 1.00000D-03 K3= 1.00000D-03 

ITERATION FREQUENCY AMPLITUDE 

0 0 . 1759634 01 2D 01 0 .45312S6762D-06 

1 0.1 754793 059D 0 1 0.451 5838346D- 06 

2 0.1752 0 96824 D 01 0 .4481 19431 4D- 06 

3 0.1 7437 1 6 09 3D 0 1 0 .43233665 1 2D- 06 

4 0.1 7447852 1 9H 0 1 0 .4403 1 76 1 04D- 06 

5 0. 1746624395D 01 0 .44232703S6D-06 

6 0. 17471 72576D 01 0 .4436675832D-06 

7 0 . 1 74 75 95835 D 01 0.4442578923D-06 

8 0 . 1 747828656 E 01 0 .444676586 1 D-06 

ND* 1 FREQUENCY® 1 .74799D OORfiD/SEC AMPLITUDE® 4.44921P-07 NIT® 8 


Figure 1-13. 
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LST SYSTEM-NUMERICAL SOLUTION OF 1 +N-*G 1 +N*N^S2= 0 . 

GAMMA* 1.3S000D 07 TO* 1 .OOOOOD-Ol JV~ 1.00000D 05 J6= 2.100001) 00 
KP= £.160 0 0 D 02 KI= 9. 7000 OH 03 H= 6. 00 GOOD 02 K0= 5.75835B 03 
Kl= 1.37102D 03 K2= 1.0 000 OH- 02 K3= 1.000000-03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 . 4 £ 6 9 0 £ 9 7 6 8 P 

01 

0 . 556644825 OD- 05 

1 

U .4260430902D 

01 

0.5563841 178D- 05 

o 

L_ 

0 . 4 c! 6> I - * 4 o 4 c! 6* l! D 

01 

0 . 556 4 35936 OD - 05 

o 

0 . 427236375 1 D 

01 

0 . 5564759926D- 05 

4 

0. 42725 12787D 

01 

0 . 557 1 7 03289D - 0 5 

5 

0.427270 092 OD 

01 

0.557 143743 ID- 05 


NB= 1 FREQUENCY* 4.27286D OORRD^SEC AMPLITUDE* 5.57162B-06 NIT= 5 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N^N*G2=0 . 

6AMMA= 1 .380001* 07 TO* 1. OOOOOD-Ol JV= l.OOOOOD 05 JG= 2.10 0 0 OD 00 
KP= 2.1600 Oil 02 K I = 9.70000D 03 H= 6. 00 GOOD 02 K0= 5 . 75S35D 03 
Kl= 1 .37102D 03 K2= 1.00 GOOD- 01 K3= 1.00 GOOD- 03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0.42723604 04D 

01 

U . 557 1617 0 63 D- 05 

1 

0 ■ 4 2 776 3 7 36 3D 

01 

0 . 556 1 946369D- 05 

2 

0.4275391 190D 

01 

0 . 555 05S5436D- 05 

3 

0.4275463586D 

01 

0. 554527 1994D-05 


MD= 1 .FREQUENCY* 4.27574D OORflD^SEC AMPLITUDE* 3.54293D-06 NIT* 3 


LST SYSTEM-NUMER I CAL 
GAMMA- 1 .33 GOOD 07 TD 
K'P- 2. 1600 OD 02 K I 
Kl* 1.37102D 03 K2 


□LUTION OF 1+N«G1+N+N*G2=0. 

1 . 0 0 0 0 0 1) - 1 j 1 J V = 1 . 0 0 0 0 0 Ii 0 5 

9.7 0 0 0 0 DOS H = 6 . 0 0 0 0 0 D 0 2 

1.00 0 0 OD 0 0 K3= 1.0000 OD- 03 


JG= 2. 100 GOD 00 
K0= 5.75835D 03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0.4 

27573661 OD 

01 

0 .55429348740- 05 

1 

0.4 

335235297D 

01 

0.54816631 1 2D- 05 

c. 

0.4 

32448 1 02 6 D 

01 

0 . 5 3 6 5 1 0 4 0 8 5 D - 0 5 

3 

0.4 

316383079D 

01 

0 . 53 03 1 6 03 1 1 D- 05 

4 

0.4 

i 1 c J t c! 9 8 1 D 

01 

0. 5272793 022D- 05 

5 

0.4 

-> 1 0 coo 6 6 3D 

01 

0. 525743243 OD- 05 

* 

0.4 

30911471 OD 

01 

0 . 52493 02355D- 05 

— <y 

f 

0.4 

3 0762285 OD 

01 

0.5243785475D-05 

P 

0.4 

3024866 15D 

01 

0 . 52369 G4245D- 05 

9 

0 .4 

3 07 06- 045 3D 

01 

0 . 5237 0 03046D- 05 

10 

0.4 

306932S43D 

01 

0 . 52369 1 5323 D- 05 


ND= 1 FREQUENCY* 4.306/6B OOPAD^SEC AMPLITUDE* 5.23683D-06 NIT* 10 


Figure 1-1 4a . 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1+N*N»G£*0. 

GAMMA* 1.380000 07 TD* 1 . 00 00 Oil- 01 .JV= 1 . 00 00 OH 05 ,JG= £.100000 00 

KP= 2.160001' 0 2 KI= 9.70000 D 03 H= 6.000000 0£ K 0= 5.758350 03 

Kl= 1 .371 Oc'B 03 K£= 1.000000 00 K3= 1. 000000-03 

ITERATION FREQUENCY AMPLITUDE 


0 

0 .430 00 00 00 OB 

01 

LI 

.58000000000-05 

1 

0 .43058318860 

01 

0 

. 58883039010-05 

£ 

0.43 0 5 £ 8 3 8 5 0 0 

01 

0 

.58884534390-05 

•j 

0 . 4 3 0 87 4 4 y 7 0 0 

01 

0 

.58457487880-05 

4 

0 .43036343530 

01 

0 

.5843841 0150-05 


NIi= 1 FREQUENCY* 4.303650 OORAB.-'SEC AMPLITUDE* 5. £41850-06 HIT* 4 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N+G1+N^N^G£=0 . 
GAMMA* 1.380000 07 TO* 1.000000-01 JV= 1.000000 05 

K P = £.16 0 0 0 0 0 £ K I = 9.7 0 0 0 0 0 0 3 H= 6.0 0 0 0 0 0 0 £ 

K 1 = 1 . 37 1 0£D 03 K£= £ . 0 0 0 0 00 0 0 K3= 1.000 0 OB- 03 

ITERATION FREQUENCY AMPLITUDE 


0 

0.43 08645 193D 

01 

0 .58418508890-05 

1 

0. 4-37596 00 OSD 

01 

0.5 185863 19 OB- 05 

IZ. 

n . 4 :':^, n fi S ij nit 

01 

0 .50473814180-05 

*J» 

■-A 

0 .43499847780 

01 

0 .49764865880-05 

4 

0.43449053800 

01 

0.49419138490-05 

5 

0. 43483 13858B 

01 

0.49844759580-05 

6 

0.434 0863 1070 

01 

0 . 49 1 58476 050- 05 

—f 

i 

0 .43398843060 

01 

0 . 4 9 0994755 00— 05 

y 

0 . 4339 0689830 

01 

0 .49065868830-05 


NB= 1 FREQUENCY* 4.338100 OORAB^SEC AMPLITUDE* 4. 9 04 150-06 HIT* 8 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N*G 1 +M^N-*6£= 0 . 

GAMMA* 1.380000 07 TO* 1. 000000-01 JV= 1.000000 05 

KP= £.16 0000 02 K I * 9.7 000 00 03 H= 6.000 0 00 0£ 

K 1 * T . 37 1 0£D 03 K£= 3 . 0 0 0 0 OB 0 0 K3= 1.000 0 00- 03 

AMPLITUDE 
0. 4904 15£ 1080-05 
0 . 486 1 0 09 1 £ 80 - 05 
0.471 1 09335£P-05 
0.46363018150-05 
0 .459 9 9 4 8 9 £70- 0 5 
0.45815783940-05 
0 .4571 7785810- 05 
0 . 4566 1331 050- 05 
0 . 45685484 1 70-05 

NO* 1 FREQUENCY* 4.361990 OORftD^SEC AMPLITUDE* 4.559970-06 HIT* 3 


ITERATION FREQUENCY 

0 0.43381008730 01 

1 U .44 U 794 733 00 01 

'*■ £ 0 . 4 3 8 7 3 3 3 4 9 0 0 01 

3 0.43755069490 01 

4 0 . 4369644 0780 01 

5 0.43666188310 01 

6 0 . 4 3 6 4 9 188 £ 5 0 0 1 

7 0 .43638155070 01 

8 0.43689394910 01 


JG= £ . 1 000 OB 00 
KO* 5.75335B 03 


JG= £.100 0 OB 0 0 
KO* 5. 75 83 5 B 03 


Figure 1-1 4b . 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1 +N*N*G2=0 . 

GAMMA* 1 .38000X1 07 TO* 1.00000D-01 .JV* l.OOOOOD 05 JG= 3.100 0 OX 00 

KF= 2.16 0 0 0 D 02 KI= 9.70000D 03 H= 6. 00 00 Oil 02 K0= 5.75835D 03 

Kl= 1 .37108D 03 K8= 4. 00 0 0 Oil 0 0 K3= l.OOOOOD- 03 

ITERATION FREQUENCY AMPLITUDE 


0 

0 . 436 1 9381 75 D 

01 

0 . 4559965684D- 05 

1 

U • 44 Jc! 44 0 6 7 D 

01 

0’ . 452593 06 1 61 ) - 05 

•Tt 

t_ 

0. 440563 1554D 

01 

0 .436 09 42 0 28 D - 0 5 

3 

0 .439 1262 01 6D 

01 

0 .428 0822 026D- 05 

4 

0.4384243852D 

01 

0 .42421 08 964 D- 05 

5 

0 .4380628732D 

01 

0 .4222494599D-05 

6 

0.43 7859 033 OD 

01 

0.421 192392 ID- 05 

\ 

0.4377394674D 

01 

0 .42058939 8 9 D - 0 5 

3 

0.437665 07 05D 

01 

0 .42 022 1826 ID- 05 


ND= 1 FREQUENCY* 4.37615D OORAD/SEC AMPL ITUDE= 4.19983D-06 NIT* 8 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1 +N*N*G2=0 . 
GAMMA* 1 .380000 07 TO* 1 . 000000-01 JV* l.OOOOOD 05 

KP= 2.16 0 0 00 02 KI= 9.70000D 03 H= 6.00000D 02 

Kl = 1 .37102D 03 K2= 5.00000D 00 K3= 1.00000B-03 


ITERATION 

FREQUENCY 


0 

0 .4376 14861 3D 

Cl 

1 

0 . 444626869 OD 

01 

2 

0.441 0 72339 3 D 

01 

O 

0. 439231 6296D 

01 

4 

0 .43834979 1 4 D 

01 

5 

0.43 73966 024 D 

01 

6 

0.4376399 04 ID 

01 

7 

0.43743983 OSD 

01 

y 

0 . 4374 03635 OD 

01 

9 

0. 4373567 032D 

01 


AMPLITUDE 
0.4139 3 3 0 0 2 7 D - 0 5 
0 .4172529128D-05 
0 . 398623937 6I> — 05 
0.38973 1285 ID- 05 
0 . 3854826622D— 05 
0 . 38332 1418 OD- 05 
0 . 3 82 1 3 92659 D - 05 
0 .381 45369 04D- 05 
0.381 035 0253D- 05 
0.38 07689 16 ID- 05 


JG= 2.1000 OD 00 
K0= 5.75835D 02: 


ND= 1 FREQUENCY* 4.37333D OORAD-'SEC AMPLITUDE* 3.80595D-06 NIT* 9 


LST SYSTEM-NUMERICAL SOLUTION DF 1 +N*G 1 +N*rHG2= 0 . 

GAMMA* 1 .38 00 OD 07 TO* l.OOOOOD- 01 JV* l.OOOOOD 05 JG= 2.1 000 OD 0 0 
KP= 2.16000D 02 K I = 9.70000D 03 H= 6.00000D 02 KC- 5.75835D 03 
Kl= 1.37102D 03 K2= 6. 0 0 GOOD 0 0 K3= l.OOOOOD- 03 


ITERATION 

0 

1 

FREQUENCY 
0 .4373325369D 
0 .4443683 083D 

01 

01 

0. 
0 . 

AMPLITUDE 

3805950594D-U5 

»■ y 6 u j 

c, 

0 .439371601 4D 

01 

0. 

356394471 ID- 05 

O 

ij . 4 3 6 b 9 Li 6 5 b y If 

01 

0 - 

3456 0 1 5384D- 05 

4 

0 . 4 3 5454 7 5 2 7 D 

01 

0. 

34 05764 07 Op- 05 

zz 

_l 

0 .434S162746D 

01 

0. 

33 799163 6 5 D— 0 5 

A 

0 .434 448564 4D 

01 

0 . 

33 b 5 3 8 £ £ 7 2 D - 0 5 

f 

0 .4342266743D 

01 

0 . 

3356649 144B- 05 

y 

0 .4341 14 7983D 

01 

0. 

335 1365017 D- 05 

9 

0 .434 08331 14D 

01 

0. 

j' y 4 y c64 1 y y D — u 5 


ND= 1 FREQUENCY* 4.34081D OORAD-'SEC AMPLITUDE* 3.34647D-06 NIT* 9 


Figure 1 -14c. 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N^N+G£=G . 
GAMMA® 1 .38 00 OL 07 TO= 1.0000GD-01 .JV= 1 . 00000D 05 

KP= £.16 0 0 0 D OS K I = 9.70 00 OB 03 H= 6.00 0 0 OP 03 

Kl= 1 .371 02D 03 K£= 7.000001' 00 K3= 1 . 000 OOP- 03 

ITERATION FREQUENCY AMPLITUDE 

0 0 .4340313381 D 01 0 .3346473373D- 05 

1 0 .4403645 049P 0 1 0 .3335789753D- 05 

3 0 . 4.385 039 1 4 6 U 0 1 0 . 3 0 1 359785711- 05 

3 0.436270179711 01 0 .2888394047D-05 

4 0 . 4 £33 1 822 OSD 0 1 0 . 2738939 0S6D- 05 

5 0.42 1653 06 16D 01 0.268392506511-05 

6 0.42 0 49353 0 6 D 01 0 .2655545194D-05 

7 0.41 9789 62 93 D 0 1 0 . 2634 1 0015 OD— 05 

8 0.418875132711 01 0 . 26 1 26 19774D- 05 

9 0 .41 89724 044D 01 0 .2606441657P-05 

10 0 .419121 1536D 01 0 .2604446450C-05 


JG= 2.10000D 00 
K0= 5.75335D 03 


MD= 1 FREQUENCY = 4.19137D OORAD-'SEC AMPLITUDE® 2. 6 043 OD- 06 NIT= 10 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N#G1+N*N*G2=0 . 
GAMMA= 1 .38000D 07 T0= 1 . 00000D-01 .JV= l.OOOOOD 05 

KP= 2. 16 GOOD 02 KI= 9.70000D 03 H= 6. 00 GOOD 02 

Kl= 1.37102D 03 K2= S.OOOOOD 00 K3= 1.000D0D-03 


ITERATION 

FREQUENCY 

0 

0.41 91 366586D 

1 

0 .432033 056 7 D 

■'i 

C 

0 .3385654674D 

o 

U . 9 9 cl i o o 9 S A B 

4 

0.41 0121 1254D 

5 

0.4049 0 6 8157 D 

6 

0 . 42 1 865 D368D 

■ 

0 . 532854 1 03 9 D 


U .49 ’J 'd >: 7 4 0 9 7 Ii 

9 

0 .45 y u 7 4 cl 9 9 0 D 

10 

0.42726331 17D 

11 

0.4 726 01 27 COD 

12 

U .4642462 0 3 2 D 

13 

0.41 02385 1 6 1 D 

14 

0.4511 00094 OD 

15 

0 .4 18504943 ID 

16 

0.433792 154 ID 

17 

0.45847 1929 ID 

18 

0.441 03 1 3 1 36 D 

19 

0 . 39 0 1 53372 OD 

20 

0 . 8 9 51 -3 5 5 -3 «3 9 D 


AMPLITUDE 

0 1 0 .26 0 4 3 0 0 9 6- 5 D — 0 5 

0 1 0 . 275544992 OD- 05 

01 0. 1537227546D-05 

01 0.1 7770982 03B- 05 

0 1 0 .201 8 1 97669D- 05 

01 0. £133 7 5 27 1 53)- 05 

0 1 0 .24428 OS 0 1 6.D- 05 

0 1 0.5121 5428 03D- 05 

0 1 0.4 09 0396 1 6 1 D- 05 

01 0. £893763479 D- 05 

0 1 0.23 3 8 9 3 0 9 5 9 D - 0 5 

01 0 . 3527236 0 1 3D- 05 

0 1 0 • c 8 U 6-8 14 2 8 r’ D — U 5 

0 1 0.171 26 07528B— 05 

0 1 0 .23 07657 1 0 1 D- 05 

0 1 0.214 399524 OD- 05 

0 1 0 . 2566335 733D- 05 

0 1 0 . 3 1 7 2396 06 4 D— 05 

0 1 0 . 25954 1 76 1 60- 05 

0 1 0.1 45709 02 02D- 05 

0 1 0.15 0 0 0364 79D- 05 


Jb= 2. 10 GOOD 00 
K0= 5.75835D 03 


NB= 1 FREQUENCY® 3.95136D OORAD'SEC AMPLITUDE® 1 .50004D-06 


NIT® 20 


Figure 1 - 1 4d . 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N*N*G2=0. 

GAMMA 3 1 .330000 0? TO= 1 . 0000 OB- 01 JV= l.OOOOOD 05 J6= S' .1 00 0 OD 00 
KP= 2.16000B 02 K I 3 9.70 0 0 OD 03 H= 6.00000B 02 K0= 5.75835D 03 
K 1 = 1 .371021' 03 K2= 1.00 GOOD- 02 K3= 1 . 0 00000-03 


ITERATION 

FREQUENCY 


AMPLITUDE 

o 

0 . 1747994849D 

01 

0 .444921 0692B-06 

1 

0 . 1750970353D 

01 

fj ,44 5 3 *Z' cl y s’ u 6' 

’ d 

0.1751 578924D 

01 

0.445861 0524 D- 06 

3 

0 . 1751663855B 

01 

0 . 4463 09 084 OD- 06 

4 

0 • 1 i 5 1 i 6* d 6 O i‘ 35 

01 . 

0 . 4463992 1 73B- 06 


NB= 1 FREQUENCY 3 1 .751830 OGRAB-'SEC AMPLITUDE 3 4.46504D-Q7 NIT 3 4 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1+M*N*G2=0 . 

GAMMA 3 1 .38 GOOD 07 T0= 1.00000D-01 JV= l.OOOOOD 05 JG= 2.10000D 00 
KP= 2 . 1 GOO OB 02 K I = 9.70000D 03 H= 6.00000D 02 K0= 5.75835D 03 
Kl= 1 .37102D 03 K2 3 1.00000D-01 K3= 1.00000D-03 


ITERATION 


FREQUENCY 


AMPLITUDE 

0 

0 . 

1751834176D 

01 

0 .4465 041 1 OSD- 06 

1 

0 . 

1 769294 195D 

01 

0.45101 01543D-06 

c. 

0 . 

1 7789 04699D 

01 

0 . 4546263 062D— 06 

3 

0 . 

17840095360 

01 

0 .456764443 ID- 06 

4 

0 . 

178666671 2D 

01 

0 .4577447044D-06 

5 

0 . 

1786379223D 

01 

0 . 45797 07749D- 06 


ND= 1 FREQUENCY 3 1 .78701D OORAD/'SEC AMPLITUDE 3 4.58193B-07 NIT 3 5 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N*N*G2=0 . 

GAMMA 3 1.38000D 07 TO 3 1 . OOOOOD-Ol JV 3 l.OOOOOD 05 JG= 2.10000D 00 
KP= 2. 16 00 OS 02 KI= 9.70000D 03 H 3 6.000000 02 K0= 5.75335D 03 
Kl= 1 .37102D 03 K2= l.OOOOOD 00 K3= 1 . 00 OOOD- 03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 .1787 0 1 0 0 1 2D 

01 

0 . 4 58 1888 08 IB- 06 

1 

U • 1 S 6 6* 4 6- b 1 o d D 

01 

0 .3314511 173D-06 

d 

0 . 1 953999 1 2 OD 

01 

0 .4 14834897 OB- 06 

o 

'W . 

0 . 2 02 1 29 085 3D 

or 

0. 4325286553 D- 06 

4 

0 .20645287 09B 

01 

0 ■ 52643980 04B- 06 

er 

•J 

0 .20891 01 053 D 

01 

0 . 55 1 4356 079D- 06 

6 

0.21 01978469D 

01 

0 . 56443 003 OSD— 06 

-y 

i 

U .21 0826652 1 D 

01 

0 .5706690202D-06 


0.21 1 0689 098 D 

01 

0.57351 03622D- 06 

9 

0.21 1 1477387D 

01 

0 .575567825711—06 

1 0 

0.21 12135359 D 

01 

0 .5764S54872B-06 

1 1 

0.21 1 249 1 8 1 3D 

01 

0 . _ 1 1 t c. 1 r* c iz 6 c. D — U 6 

12 

U . 2 1 1 27 688 1 9 D 

01 

0 . 5776432494D— 06 


ND= 1 FREQUENCY 3 2.11296D OORRD-'SEC AMPLITUDE 3 5.77976D-07 NIT= 12 


Figure 1 - 1 5a . 



LST SYSTEM-NUMERICAL SOLUTION OF 1 *■ N *G 1+ N*N*G2=0 . 

G AMM A = 1.3800CD C7 rn = l.QOOOOO-OI JV= 1.000000 05 JG= 2.100QQD 00 


KP= 2.160000 
KL* 1.37L02D 

C 2 K I = 9. 

03 K2= l. 

700 COD 
000 OOD 

C3 H= 6. OOOOOD 02 

00 K 3= i • 000 OOD- 03 

K0= 5 • 75 8350 03 

iteration 

0 

FREQUENCY 
0.21 000000 OOP 

01 

amplitude 

0.58000000000-06 


1 

0.21058192380 

Oi 

0.5 76440 82 1 8 D-C6 


2 

0.21 106834450 

Oi 

0.57779394370-06 


3 

0.21 134006370 

01 

0 . 5 784 88 53 5 8D- 06 


4 

0. 21 13458576D 

01 

0. 5787551958 D- 06 



ND = I FREQUENCY* 2.113520 OORAO/SEC AMPLITUDE* 5. 788290-07 NIT* 4 


LST SYSTEM-NUMERICAL SOLUTION OF 
GAMMA* 1.30000D 07 TO* 1. OOOOOD- 

1*N*GX+N*N*G2=0. 

01 J V* 1.000000 05 

JG= 2.100000 

00 


KP= 2.160000 C2 K I = 9. 

700000 

03 H * 6. OOOOOD 02 

KQ= 5,758350 

03 


K 1= 1.371020 03 K 2= 2. 

00 0 000 

00 K 3= 1.000 000-03 




ITERATION 

FREQUENCY 


AMPLI TUDE 




0 

0.21 135223680 

01 

0.57382942720-06 




I 

0. 219878 3766D 

01 

0.4 7457342560-06 




2 

0.22862730950 

01 

0. 5541505311 D- 06 




3 

0.23516917170 

01 

0.6233 7062680-06 




4 

0. 2392340682C 

01 

0.6677152962D-06 




5 

0.24145437180 

01 

0.6 9194946710-06 




6 

0.24255743730 

01 

0.70374038480-06 




7 

0 . 24 303 515400 

01 

0.7C88785615 0-06 




8 

0.24318808540 

01 

0.7 120244100D-06 




9 

0.24328393040 

01 

C. 7 1389960070-06 




10 

0.24334026530 

01 

C. 71510600900-06 




11 

0.243 39 3 54860 

01 

€.7 15938 66780-06 




12 

0.24342705380 

01 

0.71653880260-06 





N D* I FREQUENCY- 2.434530 00RAD/5EC AMPLITUDE* 7. I6991D-07 NIT* 12 


LST SYSTEM-NUMERICAL SOLUTION CF 
GAMMA* 1 .380 OOD 07 TO* l.OOOOOD- 
KP= 2.160000 02 K I = 9.700000 

L+N*G1+N*N*G2=0. 

01 JV= 1. OOOOOD 05 

03 H= 6.000000 02 

JG= 2.100000 
K 0= 5.75835D 

00 

03 

K 1 = 1.371020 

03 K2* 3.' 

oooooc 

00 K 3* 1. 003000-03 



ITER AT ION 

FREQUENCY 


AMPLI TUDE 



0 

0.24345273020 

01 

C. 7 1699083600-06 



1 

0.25307037250 

01 

0.67539179420-06 



2 

0.26 l 50496680 

01 

C. 74730169010-06 



3 

0.26 718 1895 70 

01 

C. 8C575802490-06 



4 

C. 27038492980 

01 

0.839567 19440-06 



5 

0.27198198200 

oi 

C. 85575344710-06 



6 

0.27266560500 

OI 

C. 86176561730-06 



7 

0.27281734610 

0 1 

0 . 8 6 39006 04 90-06 



8 

0.272891911 70 

01 

C. 86589763960-06 



9 

0.27296032510 

01 

0.86694999620-06 



10 

0.27300454510 

01 

C. 86793407340-06 



11 

0.27304218240 

01 

C.86858662080-C6 




NO* 1 FREQUENCY* 2.730710 OOPAD/SEC AMPLITUDE* 8.691500-07 NIT* II 


Figure 1-1 5b. 





1ST SY5TEM-N)UV C R 1CAL SOLUTION OF l + N*G 1+N*N* G2=0 . 

GAMMA* 1 . 3 60 COD C7 T0 = 1.000000-01 JV= l.OOOOOD 05 JG= 2*100000 00 

KP= 2. 160000 02 K I = 9.700000 03 H = 6.00Q00D 02 K0= 5.75B35D 03 

Kl= 1.371020 03 K2 = 4.00000D 00 K3= 1.000000-03 

. iteration frequency amplitude 

0 0.27307112410 01 0.8691 49 8890 D- 06 

1 0.28327612720 01 C . 869056 29920-06 

2 0.29 140 72 3 6 90 01 C. 940941 74470-06 ; - 

3 0.29652075090 01 0.9 S3123 0333D-06 

4 0. 29927995950 01 0 . 1 02 1 78 1 452 D-Q5 

5 0. 30063065 160 01 0. 1C34369396D-05 

6 0.30111433820 01 Q . 1 0 3300 19 3B D-05 

7 0.30116053060 01 0. 1 C4028 546 1 0-05 

8 0.30 12294-6 L4D 01 0. 104091 6602 0-05 

9 0. 30 126669600 01 0.10419866650-05 

10 0.30130411490 01 0 . 1 C42 5 8 56 02 D-05 

ND= 1 FREQUENCY = T.01333D OORADTS’EC AMPLITUDE 1.043210-06 N IT= 10 


LST SYSTEM-NUMERICAL SOLUTION OF 
GAMMA* 1.38CT0D C7 T0= l.OOOOOD- 
KP= 2. 160000 02 K 1= 9.700000 

1+N*G1*N*N*G2=0 . 

01 JV= l.OOOOOD 05 

03 H= 6. 0000 OD 02 

JG= 

K0= 

2.10000D 00 
5.75835D ‘03 

K 1= 1. 371020 

03 K2= 5 • 000 GOD 

00 K 3= 1.000QOD-03 



ITERATION 

FREQUENCY 


amplitude 



0 

0.30 133297290 

01 

0.1 043207423D-05 



1 

0.31 158859480 

Cl 

0.1C6599 22250—05 



2 

C. 31957166540 

01 

0.1 1425688390-05 



3 

0.3245957530C 

01 

0.1 1987198770-05 



4 

0.327 35744 7 ID 

01 

0.1230191890 D-05 



5 

0.32 86982993D 

01 

0 . 1 244 194422D— 05 



6 

0.32899643840 

01 

0. 12456812450-05 



7 

C.329Q445550D 

01 

0. 12487315320-05 



8 

0.3291 33 30090 

01 

0. 12493769450-05 



9 

0.32917659820 

01 

0 . L 2508 26 1 72 D-05 



10 

0. 32922447670 

01 

0.12515Q35020-05 


kM T _ i ri 


< N0= 1 FREQUENCY = 3.292610 OORAO/SEC AMPLITUDE^ 1.252440-06 N1T= 10 


Figure 1-1 5c . 



L ST SYSTEM-NUMERICAL SOLUTION OF" 
GAMMA= 1 . 38C T CD 07 TQ = U00030D- 
K P = 2.160000 02 K I = 9.700C0D 

1+N*G1+N*N*G2=0 . 

01 JV= l.OOOOOD 05 

03 H = 6.00000D 02 

JG= 2.10C00n 
K0= 5.758350 

00 

03 

K I s 1.371020 

03 K 2= 6 • i 

COOCOD 

00 K3= 1.000000-03 



ITEP AT ION 

FREQUENCY 


AMPLI TUDE 



0 

0.32926099030 

01 

0.125244 1749 C- 05 



1 

0,33931 182040 

01 

0. 1288015344 D-0 5 



2 

0.341 664 8CC4D 

01 

0.13821226540-05 



3 

0.35324191450 

01 

0. 1 4555935280-05 



4 

0. 3564585758D 

01 

0.14992048020-05 



5 

0.35782924460 

01 

0.15150643390-05 



6 

0.35799595170 

“ 01 

0 . 1 5 1789 12160-05 



7 

0.35810681380 

01 

0. 15209709710-05 



8 

0.3582153412D 

01 

0. 1 £230940250-05 



9 

0.35 8299 9 1 1 3D 

01 

0.1525 1788380-05 



10 

0.35 837 76 544D 

01 

0.1 52682 78690-05 



11 

0. 3584436003D 

01 

0. 1528349560 D-05 




ND= 1 FP EO'JENC Y = 3.58503D CORAD/SEC AMPLITUDE 3 1. 529650-06 NIT = 11 


LST SYSTEM-NUMERICAL SOLUTION OF 1 + N* G 1 + N* N* G2 = 0 . 

6AWMA= 1.38C00D 07 10= i.OGOCGD-Ol J V= l.OOOOOD 05 JG= 2. 10000D 00 

KP— 2 « 1 60 COD 02 Kj_= 9.70000D 03 H= 6.000000 02 K0= 5.75835D 03 

K 1 = 1.37102D 03 K2= 7.C00OGD 00 K3= 1.00000D-03 


ITER AT I CN 

FREQUENCY 


AMPLI TUDE 

0 

0.35050294550 

01 

0.15256501160-05 

1 

0.36835499180 

01 

0. 15758265830-05 

2 

0.37823023 170 

01 

0.1 7157825650-05 

3 

0.38603849890 

01 

0. 18458471720-05 

4 

0.39100843630 

01 

0.19306176490-05 

5 

0.39215573620 

01 

0.1 5473093000-05 

6 

0. 392795 50 lOn 

01 

0. 15590148480-05 

7 

0.393 1729440D 

01 

0. 19900152340-05 

8 

0.39497523 19D 

0 1 

0 . 2 C 14 02 1481 D-05 

9 

0.39509208140 

01 

0.2 024 24 24 2 ID- 05 

10 

0.39547837670 

OL 

0. 20274338040-05 

11 

0.3956024832D 

01 

0.2C3302U970-05 

■ 12 

0 . 395788905 7D 

01 

0 .20357573530-05 

13 

0.395914 2055D 

01 

0.20396LS541D-05 

14 

0.39605362450 

01 

0.2C4246 5906D-05 

15 

0.39617 190700 

01 

0.20455586640-05 

16 

0.39628987570 

01 

0.2C482462340-05 

17 

0.39639796620 

01 

0.2 C50905635Q-05 

18 

0.396502080 10 

01 

0. 20533612810-05 

19 

0 .39660004090 

01 

0.2C55723118D-05 

20 

0 . 396693 6 16 1 D 

01 

0.2057952294D-05 


Np= J. FRFOU£NCY= 3.96694D OORAD/SEC AMPLITUDE 3 2.057950-06 NIT= 20 


Figure l-15d. 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1+N*N*G£=0. 

6 ft M M A = 1.38 u 0 0 D 07 T □ = 1 . 0 0 0 0 0 D - 0 1 J V = 1 .0 0 0 0 0 0 05 ‘ JG* £ . 1 0 0 0 0 D 0 0 
KP= 2.160 00D 02 KI= 9.700000 03 H= 6.00000D 02 K0= 5.75835D 03 
Kl* 1 .37102D 03 K2= 1.0 0 GOOD 0 0 K3= £.000000-01 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 . 39 0 0 0 0 0 0 0 0 B 

01 

0.42 0 0 0 0 0 0 0 OD- 05 

1 

0.391 1764448D 

01 

0.41 98 1 54369D- 05 

2 

0.391 04454650 

01 

0.41 856 1 02960-05 

3 

0.39 0 7 2 8 8 8 6 2 D 

01 

0.41 755 050 1 7D- 05 

4 

0 .39 097 1 59460 

01 

0 .41742228590-05 


NO 3 1 FREQUENCY= 3.31 004D OORRD/'SEC AMPLITUDE 3 4.1737SD-06 NIT= 4 


LST SYSTEM-NUMERICAL SOLUTflON OF 1 +N^G1+N*N#G£=0 . 


GAMMA 3 1.38 00 OB 07 TD 3 1. 

00000D- 

01 -JV 3 1.0 000 OB 05 

JG= 2.1000 OD 

KP= 2. 1600 on 

0. 

2 K 1 3 9. 

700000 

0 3 H= 6. 0 0 00 OB 02 

K0= 5.758350 

K 1= 1 .371020 

0 

J K8= 8 . 

00 GOOD 

00 K3= 2. 0000 OB- 01 


ITERATION 


FREQUENCY 


AMPLITUDE 


0 

n . 

391 0 038835 D 

01 

U • 4 1 1 8 1 ' o 1 4220— U 5 


1 

0. 

39443586 120 

01 

0 . 4 03429 1 9270- 05 


O 

u* 

0. 

39 02 4 9399 7 B 

01 

0 .38345292230— U5 


3 

0 • 

3880015397B 

01 

0.37384458390-05 


4 

o . 

386932556 7 B 

01 

0.369295 1 6350- 05 


5 

0. 

386396 'Jo U 1 B 

01 

0.36701442950-05 



0. 

3861 01 87230 

01 

0 . J 6 5 7 8 6 56650 — U 5 

• 

i 

0 - 

3859279764D 

01 

0.36507819460-05 


p 

0. 

33581 943390 

01 

0 . 36464 1 02530-05 


9 

0. 

38574841540 

01 

0.36435442240-05 



ND= 1 FREQUENCY 3 3.S5700D OORRB'SEC AMPLITUDE 3 3. 641570-06 NIT* 9 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N*G 1 +N*N*G2= 0 . 

GAMMA 3 1 .330 COD 07 TO 3 1.00000D-01 JV* 1.00000D 05 JG 3 2.10000D 00 
KP= 2.160 0 OB 02 KI = 9.70000D 03 H= 6.00000D 02 K0= 5.75S35D 03 
K 1 3 1 .371 02 B 03 K2= 3 . 0 0 0 0 OD 0 0 K3= £ . 0 0 0 0 0D- 0 1 


ITERATION 

FREQUENCY 



AMPLITUDE 

0 

0.38569978820 

01 

0 

.36415 6 9 6 0 8 0 - 0 5 

1 

0.33681518770 

01 

n 

.34463564510-05 

cl 

0.37951401390 

01 

0 

.31782072220-05 

3 

0.37447438500 

01 

0 

. 3 0 8 3 o 8 5 8 1 8 0 — 0 5 

4 

0.37202850470 

01 

0 

.295 0 6940830- 0 5 

tr 

._i 

0.37085265180 

01 

0 

.2914 3 9 9 7 0 3 0 - 0 5 

K 

0 . 3 7 0 0 7 5 6 5790 

01 

0 

.28919134620-05 

7 9 

0 . 3 6 9 6 8 9 3 5 8 7 0 

01 

o 

.23792088110-05 

s 

0.36894108510 

01 


.28646457450-05 

9 

0 .36924498290 

01 

o 

.28629051 160-05 


ND= 1 FREQUENCY 3 3.692(40 OORRD/SEC AMPLITUDE 3 2.86133D-06 NIT 3 9 


Figure 1-16. 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G L+N*N*G£=0 . 

GAMMA- 1 .38 GOOD 07 TD= 1 . 00 GOOD- 0 1 JV= 1 . 0000GB 05 JG= 2- 1 000 OB DO 
KP= 2 .1600 OB 02 K I = -9.70 0 0011 03 H* 6. 00 00 OB 02 K0= 5.75S35D 03 
Kl= 1 .371020 03 K2= 1.000000 00 K3® 2.000000-01 


ITERATION 

FREQUENCY 



AMPLITUDE 

0 

0.2400000 0 0 OD 

01 

o 

. l 0 00 00 0 0000-05 

1 

0.24094794410 

01 

n 

. 1 0467125640-05 

2 

0.2399290724D 

01 

0 

. 1 033304 1690-05 

3 

0.S394536273D 

01 

o 

. 1 035725356B- 05 

4 

0.2390:56496711 

01 

0 

. 1 0313639430-05 

c 

. J 

0.2395 07685 OD 

01 

0 

• 1 046663 0 1 70— 05 

6 

0.23984 04 13SB 

01 

0 

. 1 051367 0190-05 

♦' 

024 n<W 3 9 82611 

01 

0 . 

1 0424853980-05 

8 

0.2-392532346 D 

01 

0 

.10313138350-05 

9 

0.23929 0461 3D 

01 

0 

. 1 0373983360-05 

10 

0.2393401 1 18D 

01 

0 

. 1 04067441 10-05 

11 

0. 23936 03834D 

01 

0 

. 1 040720365D- 05 


ND= 1 FREQUENCY* 2.393390 OORAD^SEC AMPLITUDE 3 1.041350-06 NIT- 11 


LST SYSTEM-NUMERICAL SDLUTIDM CF 1+M*G1+N*M*G£=0 . 

GAMMA- 1.380000 07 TO- 1.00000D-01 . JV- 1.000000 05 JG= 2.100000 00 
KP* 2.160000 02 K I = 9.700000 03 H= 6. 00 00 00 02 K 0= 5.753350 03 
Kl= 1.371020 03 K2= 2.000000 00 K3= 2.000000-01 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 . 239 3886 1350 

01 

0 . 1 0413517910-05 

1 

G.° 5455332730 

01 

0.11151902020-05 

2 

0.26539495670 

01 

0 . 12224496630-05 

3 

0.27300435520 

01 

0.12983743120-05 

4 

0.27687369050 

01 

0 .13413839240-05 

5 

0.27876453620 

01 

0 .13612073700-05 

6 

0.27945522190 

01 

0.13666903380-05 

7 

0. 279533262 OD 

01 

0 . 1363716 0 OSD- 05 

3 

0 .27961369060 

01 

0.13700211050-05 


NO- 1 FREQUENCY* £.796790 OORAB/'SEC AMPLITUDE* 1 .371270-06 NIT= 3 


LST SYSTEM -NUMERICAL SDLUTIDM DF 1+N*G1+N*N*G£=0 . 
GAMMA* 1.330000 07 TO 1.000000-01 JV= 1.000000 05 

K P= 2.160000 02 K I * 9.700000 03 H= 6.000000 02 

Kl* 1.371020 03 K2= 3.000000 00 K3= £.000000-01 


JG= 2-1 00 ODD 00 
K0= 5.758350 03 


ITERATION 

FREQUENCY 

0 

0.2796785954D 

1 

0 .29578095821) 

2 

0.3082893 0200 

3 

0 .317240 09310 

4 

0 . 3 c 2 9 o 1 U 9 0 8 B 

5 

0.32552671090 

6 

0.3260443823D 

—> 

r 

0 . 3 c! 6 4 IJ 1 37260 

S 

0 ,32667987780 

9 

0 .32692864 110 

10 

0.32713712760 

11 

0 .327323 03270 

12 

0 .32743559300 

13 

0 .327ee£SflD6D 

15 

0 . J c! 7 y 3 Uc r 2 1 D 

16 

0 .32793762540 

17 

0.32303577690 

13 

0.32817582330 


AMPLITUDE 

01 0 .13712719940-05 
01 0. 15 029 14739 B- 05 
01 0. 16542507510-05 
01 0 . 1 77936632 OB- 05 
01 0. 13654335140-05 
01 0.19 0 03 037390- 05 
01 0 . 19083475590-05 
01 0 . 1 91670727 1 D— 05 
0 1 0.1 9235 0978 OD- 05 
01 0 . 19287213530-05 
0 1 0.19 33362 076 0-05 
0 1 0.1 93 73 06384 0-05 
0 1 0 . 194 032676 OD- 05 
01 0.19439235580-05 
0 l 0 .19492348150-05 
01 0 . 19515184270-05 
01 0.19535996010-05 
01 0.19555080650-05 


NO- 1 FREQUENCY- 3.28259D OORAD^SEC AMPLITUDE* 1.957260-06 NIT* 13 


Figure 1-17. 
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LST SYSTEM-NUMERICAL SOLUTION 'DF 1+N^G1+N*N*G£=0 . 

GAMMA® 1.33 GOOD 07 TD= 1 . 0 0 0 0 OD- 0 1 J7= 1.0 000 OD 05 JG= 2.100000 00 
KP= £.160000 02 KI= 9.700000 03 H= 6.000000 02 K0= 5.753350 03 
Kl= 1.371020 03 K2= 1.000000-03 K3= 1.000000-03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 ,430000 0 00 00 

01 

0 ■ 55 0 U iJ 0 0 0 0 OD- 05 

1 

0.4271 0360 58 D 

01 

0 . 55 08 161 796D- 05 

■™l 

0.4265736694D 

01 

0 .5527424953D-05 

3 

0 . 42 6 y 7 5 0 6 4 c D 

01 

0 . 5552 1 97589D- 05 

4 

0 .4290753 07 ID 

01 

0 . 5645724 078D- 05 

er 

0 .427801 5463D 

01 

0 . 56 04 1 79332D- 05 

6 

0.4271 0731 96 D 

01 

•0 .558248041 SB- 05 

7 

0.42724761 7 OD 

01 

0 . 55792 1 6 0 1 7D- 05 


ND= 1 FREQUENCY® 4.273810 0 ORAD-'SEC AMPLITUDE® 5.578840-06 MIT= 7 


LST SYSTEM-NUMERICAL SOLUTION OF l+N+Gl+N*N*b£=0 . 

GAMMA® 1.380000 07 TC® 1.00000D-01 - JV= 1.000000 05 JG= 2.10000D 00 
KP= £.160000 02 K I = 9.700000 03 H= 6. 00 GOOD 02 KQ= 5 . 75835B 03 
Kl= 1.371 02D 03 K2= 1.00000D-03 K3= 1.00000D-02 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0.427331 1433D 

01 

0 . 557S338989B- 05 

1 

0 .42651 6528 OD 

01 

0 .556156£337D-05 

. £ 

0.42 56391 9 08 D 

01 

0,553985473411-05 

3 

0 .4253571 91 5D 

01 

0 .553766*56511-05 


ND= 1 FREQUENCY® 4.25947D OCR AD. •-'SEC AMPLITUDE® 5. 5374 OD- 06 NIT= 3 


LST SYSTEM-NUMERICAL SOLUTION DF 1+N-*G1+N^N*G£=0 . 

GAMMA® 1 .38000D 07 TO® 1.00000D-01 JV= l.OOOOOD 05 JG= 2.10000D 00 
KP= 2. 16 GOOD 02 K I = 9. 70 GOOD 03 H= 6.00000D 02 K0= 5.75835D 03 
Kl= 1.37102D 03 K2= l.OOOOOD— 03 K3= 1.00 GOOD- 01 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

U .425947337 1 D 

01 

0. 5537397 133D- 05 

1 

0.41 962 12759D 

01 

0 . 53853233 09D- 05 

*_ 

0.41 54989 1 28D 

01 

0 . 526532 05 09D- 05 

3 

0.41301 48775D 

01 

0 . 5 1 9 09493 1 6D- 05 

4 

0.41 152 1633 OD 

01 

0 .51465662950-05 

cr 

0.41 07 03 1 7 78 D 

01 

0.51c' c c 6 6 3 1 D— iJ 5 

6 

0.41 092 8 6 4 2 2 D 

01 

.512 0752653D- 05 


ND= 1 FREQUENCY® 4 . 1 1 028D GGRAD-'SEC AMPLITUDE® 5.12006D-06 HIT® 6 


Figure 1-1 8a. 
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LSI SYSTEM-NUMERIAAL SOLUTION OF 1+N*G1+N*N*G2=0. 

GAMMA® 1 .3800011 0? TO® 1 . L'i 0 0 0 CD - 0 1 JV= 1 . 00 00 OP 05 JG= £.10 GOOD 00 

KP= £ . 1 SOO CD 0£ KI= 9.70 0 0 0 D 03 H- 6. 0000 OH 0£ K0= 5.753351) 03 

K 1 = 1.371 020 03 K£= 1 . 0 0 0 000- 03 K= £ . 0 0 0 0 00- 0 1 


ITERATION 


FREQUENCY 


AMPLITUDE 


90 

0.41 1 028451 40 

01 

0.51300573570-05 

1 

0.403301 988 ID 

01 

0 .49409485700-05 


0.39304746390 

01 

0 .4799 1 383650- 05 

o 

■_* 

0.39483978050 

01 

0 .471 04448530-05 

4 

0 . 39295355330 

01 

0 .465678591 10— U 5 

5 

0.39303408150 

01 

0.46314186740-05 

6 

0 .39336047060 

01 

0 .46386 028680— U5 


ND= 1 


FREQUENCY® 


'70 OORflD^SEC 


AMPLITUDE® 4.637660-06 NIT= 6 


LST SYSTEM-NUMERICAL SOLUTION DF 1+N^61+N*N*G2®0. 

GAMMA® 1 .380000 07 T0= 1.000000-01 .JV= 1.000000 05 -JG® £.100000 00 
KP= £.160000 0£ KI= 9.700000 03 H= 6.000000 0£ K0= 5.75S35D 03 
K 1= 1 .371 0£D 03 K£= 1.00000-03 K3= 3.000000-01 


ITERATION 


FREQUENCY 


AMPLITUO 

0 

0 . 

39237667770 

01 

0 .46376636930-05 

1 

0. 

3 8 3 2 6 8 9 £ 9 3 D 

01 

0 . 443 09736 1 90-05 

2 

0. 

37709777010 

01 

0.4371 9358320- 05 

-^1 

0 . 

37318353090 

01 

0.41703164810-05 

4 

0. 

37 067333 1 30 

01 

0.41 0 6969 0 4 0 0 - 0 5 

rr 

0. 

37 0£ 0438470 

01 

0 . 4 09 04 08377 0-05 

6 

0. 

37019500790 

01 

0.40847537970-05 


0, 

370451 17770 

01 

0.40846149610-05 


ND= 1 FREQUENCY® 3.704950 OORAD^SEC AMPLITUDE® 4.085730-06 HIT® 7 


LST SYSTEM-NUMERICAL SOLUTION OF 1+H*G1+N*N*6£=0 . 
GAMMA® 1.380000 07 TD = 1. 000000-01 JV= 1.000000 05 

KP= £.160000 0£ KI= 9.700000 03 H= 6.000000 0£ 

K 1 = 1.371 0 £ 0 0 3 K £ = 1 .00 0 0 0 0 - 0 3 K 3= 4 . 0 0 0 0 0 0-01 


ITERATION 


FREQUENCY 


0 

0 . 

37049481360 

01 

0 - 

4 

1 

0 . 

359 06 0 03 OiBIi 

01 

0 . 

z 

c 

0 . 

35035157780 

01 

0. 

3 

o 

Ij . 

34516064410 

01 

0. 

£ 

4 

0 * 

34124206450 

01 

0 . 

s 

CT 

0 - 

341 64935480 

01 

0. 

£ 

f. 

0. 

34 08329 0 080 

01 

0 . 

£ 


0 , 

34 1 223446 1 0 

01 

0 . 

£ 

y 

0 - . 

34123943870 

01 

0. 

£ J 


AMPLITUDE 
0857£68£4D- 05 
8 5 £ 3 £ 6 0 0 5 D — 0 5 
6 5 7 £ U 654 1 0 — U 5 
5333051 450- 05 
4367563480— 05 
4331691850-05 
4 1 87469390—05 
4 1 8 4 3 6 8 570 - 0 5 
4 171£9£7£D-G5 


JG= £.100000 00 
K'O® 5.758350 03 


N0= 1 FREQUENCY® 3.41234D OORRO/SEC AMPLITUDE® 3.416410-06 NIT® 8 


Figure l-18b. 



4 j 


GAMMA® 
KP= 2 . 
K 1 = 1 . 


GAMMA® 
KP= 2. 

Kl= 1 . 


TEM-MUMERICAL SOLUTION OF 

1 +M*G 1 +H-*N+G 2 = 0 . 

1 .33 

noon • 07 TO= 1 . 

0 00 0 0 D - 

01 JV= 1.000 

lb 00 on 02 KI= 9 . 

7 0 0 0 OD 

03 H= 6.0000 

371 02 

H 03 K 2 = 1 . 

0 0 GOOD- 

03 K 3 = 5.000 

AT I ON 

FREQUENCY 


AMPLITUDE 

0 

0 . 34 1 234 1 33 1 D 

01 

0 . 34164 0 5 2 0 3 0 - 0 5 

1 

0 . 3241 06482 OD 

01 

0 . 3 0 3 3542 9 370 - 0 5 

o 

0 . 3 077393747 D 

01 

0 . 27339204010-05 


0 . 23049952971 ! 

0 1 

0 . 221 0 7925 6 00 - 9 5 

4 

0 . 292 059 0 09 2 D 

01 

0 . 23945924140-05 

c r 
J 

0 . 2947 S 96527 D 

01 

0 .£ 7467997 £ 4 Ii —05 

6 

0 . 3223599934 D 

01 

0 . 3 1 5 3 39991 3 0 - 0 5 

t 

0.3 0397343 140 

01 

0 . 2 6 9 3 2 3 33770 - 0 5 

o 

U - cl W c 1 4 6 i cl i - 6 Jj 

01 

0 . 22470343270-05 

9 

0 - 3 0 U 1 55264 OD 

01 

ij m d 5 4 cl y 1 6 6 5 5 D — U 5 

10 

0 . 309770 0472 D 

01 

0 .272 0 4 3 3 3 3 4 0 - 0 5 

11 

0.27495064330 

01 

0 . 1 99 0339 06711 - 05 

12 

0 .29430331640 

01 

0 . 24039329210-05 

13 

0 .32957239420 

01 

0 . 31301395430- 05 

14 

0.31596592220 

01 

0 . 3 0 692752820 - 0 5 

15 

0 .30226119770 

01 

0 . 26798983350-05 

16 

0 .d 85893 05570 

01 

0 . 23242329560-05 

17 

0 .4634155030 

0 1 0 

.• 338359 1 5440 - 05 

IS 

0 .33203920300 

01 

0 - 3 0 576 - 0 6 1 5 £ D — U 5 

19 

0 .3141 3335720 

01 

0 . 276563466 u D — 0 5 

20 

0.29412203530 

01 

0 . 2394984 0450- 05 

FREQUENCY® 2.94122 

D OORAD 

.-•SEC AMPLITUDE 1 

'STEM- 

NUMERICAL SOLUTION DF 

1 +N* 6 1 +N*N-*G£= 0 

1 . 33 

0000 07 TO® 1 . 

0 0 0 0 OD- 

0 1 JV= 1 . 0 0 0 

16 0 0 OB 02 KI = 9 . 

7 0 0 0 0 D 

03 H= 6.0000 

37102 

D 03 K 2 = 1 . 

oooooo- 

03 K 3 = 6.000 

AT ION 

FREQUENCY 


AMPLITUDE 

0 

0.29412203530 

01 

0 . 2394 9 S 4 0 4 5 0 - 0 5 

1 

0.31633142680 

01 

0 . 362719591 00-05 

£ 

0.270330761 ID 

01 

0 . 2 0 5 2 9 9 0 6280 - 0 5 

3 

0 .31605070420 

01 

0 . 2 8 6 3 9 0 92290 - 0 5 

4 

U . £ 695 12494 4 D 

01 

0 . 2 0 1304 3 2 3 2 D - 0 5 

5 

0 .31155443340 

01 

0 . 27636719450-05 

6 

0 . 8 b 338 U 9 7330 

01 

0 . 1 937533 08 10- 05 

—y 

r 

0 . 299936532 1 0 

01 

0 . 25679907640- 05 


0 .19733495540 

01 

0 . 6 4 8 3 5 4 3 0 940 - 0 6 

9 

0 .20535601340 

01 

0 9 0 9766 3 0 430 - 0 6 

1 0 

0 .225631 13600 

01 

0 . 12295239990-05 

u 

0 .24350314170 

01 

0 . 15921056760-05 

12 

0.29521225370 

01 

0 . 21 998709990- 05 

13 

— U * 5 7 6 o’ r=< 7 y 7 6 9 D 

01 

0 . 13192165590-04 

5 *763 

630 0 0 -1 .31 

922 D- 05 


14 

fi „ *999999996 D 

0 0 

ft . 3 3 8 8 8 8 8 q 8 A 0 0 0 

15 

0 . 1 0 0 0 0 0 0 0 0 OD 

01 

0 . 1 00000000 OD 01 


i 0‘ 
02 


JG= 2.10 0 0 Oil 00 
K 0 = 5. 758350 03 


-JG= 

K0= 


10 00 00 00 
753350 03 


HD® 1 


FREQUENCY® 1 . 00 C 00 D OORAD-'SEC AMPLITUDE® l.DOOOOD 00 NIT® 15 


Figure 1-1 Sc 
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LST SYSTEM-NUMER I CAL SOLUTION OF 1+N*G1+N^N*G2=0. 

GAMMA® 1.38 GOOD 07 TO= 1 . 0 0 GOOD- 01 JV= 1.0 0 GOOD 05 .JG= 8 . 1 000 CD 0 0 
KP= 2.16000D 02 KI = 8. 7 0 GOOD 03 H= 6.000001* 02 K0= 5 . 758351* 03 
Kl= 1.37102B 03 K2= 1 . 0000 OB- 03 K3= 1.G0000D-03 


ITERATION 


FREQUENCY 


AMPLITUDE 

0 

0 . 

1 7 0 0 0 0 0 0 0 0 D 

01 

0.4500 0 0 0 0 0 OD- 06 

1 

o. 

1719691894D 

01 

0 . 4377777686D-06 

•■s 

L. 

0. 

173381721 OD 

01 

U .44 0 1 73 07 66D- U6 

o 

0. 

1 741 65639 OD 

01 

0 . 4426 0 1 8322D- 06 

4 

0. 

1 74589 1 63 1 D 

01 

0 .444 1 635449D-06 

CT 

■Jt 

0 . 

1 748065537D 

01 

0 .4447812824D-06 


0 - 

1748141543D 

01 

0 .4451 726996B-06 


ND= 1 . FREQUENCY® 1.74823B OORRDxSEC AMPLITUDE® 4. 452540-07 NIT® 6 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N.G1+N*N*G2=Q. 

GAMMA® 1 .38000D 07 TO® 1.00000D-01 JV= 1.00000D 05 JG= 2.10000D 00 
KP= 2 . 1600 00808 KI = 9.70Q0DD 03 H= 6. 00 00 OH 02 K 0= 5 . 75835 D 03 
Kl= 1.371 OSD 03 Kc= 1 . 0 00 GOD- 03 K3= 1.0 0000D-02 


ITERATION 

FREQUENCY 


5 AMFLITIJDE 

0 

0 .174822728 OD 

01 

0 . 4452536775D- 06 

1 

0 . 1752126778B 

01 

0 . 45 0279625 OD- 06 

2 

0\ 1 75365 0 1 9 OD 

01 

0 .45297748390- 06 

3 

0.1754421 161 L 

01 

0 . 4544379 0 07D- 06 

4 

0 . 1 754838359D 

01 

0.4552138844D-06 

• *r 

-J 

0. 17551 72 026D 

01 

0 .4556845266D-06 

l£ 

0. 1755362634D 

01 

0.4559779957D-06 


ND= 1 FREQUENCY® 1 .75550D OORRD-'SEC AMPLITUDE® 4.56178B-07 NIT® 6 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N*61+N*N*G2®Q. 

GAMMA® 1 .38 GOOD 07 TO® 1 . 0000 OH- 01 .JV= 1.0 0 GOOD 05 JG= 2.1 0 00 OD 00 
KP= 2.16 0 0 OD 02 KI = 9 .70 GOOD 03 H= 6.00000D 02 K0= 5.75835B 03 
Kl= 1.371 OSD 03 K2= 1.00000D-03 K3= 1.00000D-01 


ITERATION 


FREQUENCY 


AMPLITUDE 

0 

0,1 

7554971 9 ID 

01 

0.4561 7849 94 B- 06 

1 

0 . 1 

791779346D 

01 

0.51 016 12925D-06 

£ 

0.1 

813712215D 

01 

0 .542 01 8421 4D- 06 

•n 

0 . 1 

826 03 1811 D 

01 

0 .559583 0487D- 06 

4 

0 - 1 

C' -'c6c 1 r - 8 1 D 

01 

U • 5 o 7 o 3 7 3 U 4D— V 8 

cr 

-J 

0. 1 

835348 0 0 ID 

01 

ij .57 3 8 8 0 3 5 €• 0 D — 0 6 

ti 

0 . 1 

83642 082 OD 

01 

U .57687422 3 8 D - 0 6 

T7 

l 

0 . 1 

837388 0 02 D 

01 

0.57801 7 1524D-06 


0.1 

o z'' { 8 c' 3 6 cz *Z' D 

01 

0.57891 074070-06 

9 

0 . 1 

838 1 7472 1 D 

01 

0 . 5793698723D- 06 


ND= 1 FREQUENCY® 1 .83840D OORflIL'SEC AMPLITUDE® 5.79735D-07 NIT® - 9 


Figure 1-1 9a. 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G 1 +N*N^G2=0 . 

GAMMA 3 1 .380000 07 TO= 1.000000-01 ,JV= 1.00000D 05 JG= £. 10 0 0 GO 00 
KP= £.16 0 0 0 D 0£ KI= 9.70 0 0 0 D 03 H= 6.000000 0£ K0= 5.75835D 03 
K 1 = 1.371 0£II 03 K'2= 1 . 0 0 0 0 0D- 03 K3= £.0000 0D- 0 1 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 . 1 8334 0 0 034 D 

01 

0.579735429 IB- 06 

1 

0 . 1886569754B 

01 

0 . 65 1 1 742 1 56 B— 06 

o 

0 .191 74334 OSD 

01 

0 . 695722 0257B- 06 

3 

0 . 193654654£B 

01 

0 .7211 374 0 1 1 B- 06 

4 

0 .1 92794548 OD 

01 

0 .74994 184330 > 06 

5 

0 . 1 9388228860 

01 

0.746801 0454B-06 

£ 

0.194751 51 77D 

01 

0 .75003480810-06 

—> 

0 . 19523442 39 D 

01 

0 .75210338960-06 

s 

0 . 1 954077725D 

01 

0. 7522564 047B- 06 


ND= 1 FREQUENCY 3 1.954140 OORAD^SEC AMPLITUDE 3 7.529320-07 NIT = 3 


LST SYSTEM-NUMERICAL SOLUTIOF OF 1+N*G1+N*N*G£=0 . 

GAMMA 3 1 . 3 3 0 0 0 D 07 T D = 1 . 0 0 0 0 0 D - 0 1 J V= 1 . 0 0 0 0 0 D 05 J G = £ . 1 0 0 COD 0 0 
KP= £.16 0 0 0 D 0£ KI= 9.70000D 03 H= 6.000 0 0 D 0£ K0= 5.75835D 03 
Kl= 1.371 OSD 03 K£= 1 .00000D-03 K3= 3.00000D-01 


ITERATION 
3 0 

FREQUENCY 
0 . 1 954 1 36225B 

01 

AMP I TUBE 

0 . 75233247440— 06 

1 

0 . £ U 1 5669259D 

01 

0 .84372498900-06 

2 

0 .20598126840 

01 

0 .9 03421 £ 1 990-06 


0 .2 063822449D 

01 

0 • 96c' 08 07 45 30— U6 

4 

0 .2 0345732350 

01 

0 . 97 1 92254400- 06 

5 

0 . £ 0964 1 4 1 1 i_B 

01 

0 .97 8 Ei £ 3 0 9 8 90—06 

6 

0.21 025421410 

01 

0 . 98228 0 076 1 0- 06 

7 

0.21 040713970 

01 

0 . 98239542 1 90— 06 


ND= 1 FREQUENCY® £ . 1 042 OD OORAD' Q'EC AMPLITUDE 3 9.83224D-05 NIT= 7 


Figure 1-1 9b 
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LST SYSTEM-NUMERICAL SOLUTION OF 1 +N+G1 +N*N*G2=0 . 

GAMMA* 1 .38 00 CD 07 TQ= 1 . 00 0000-01 .JV= 1 .0 0 00 GO 05 J6= 2. 1 000 00 00 
KP= 2.160000 02 KI = 9.700000 03 H= 6.000000 02 K0= 5.758350 03 
K 1 = 1 . 371 020 03 K2= 1 . 0 0 0 0 00-0 1 K3= 1.000 0 OO- 0 1 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0.40 0 0 0 0 0 0 0 0 0 

01 

0 „ 5 n n n n o ij n n oil 

1 

M • 4 Ll 9 o 4 c‘ 4 c* 9 D 

01 

0 .5060816759D 

o . 

0 .4096 0028750 

01 

0 .50751469700 

■-i 

o 

0.41058525050 

01 

0 .50770691710 

4 

0,41 0994741 OO 

01 

0 .50768339340 


NO* 1 FREQUENCY* 4.109960 OORAD-'SEC AMPLITUDE* 5.076300-06 NIT* 4 


LST SYSTEM-NUMERICAL SOLUTION OF 1 + N*G1+N*N-»G2= 0 . 

GAMMA* 1.380000 07 TO* 1.000000-01 JV= 1.00 GOOD 05 JG= 2.100000 00 
KF'= 2.160000 02 KI = 3.700000 03 H* 6.000000 02 KO* 5.758350 03 
Kl= 1.371020 03 K2= 1. 000000-01 K3= 2.000000-01 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0.41 099637520 

01 

0 .50763027390-05 

1 

0 • 4 08 1 i 09 r 1 cD 

01 

0 .489728 1 3150- 05 

2 

0.39800087040 

01 

0.47553569810-05 

O 

O . 3 9 4 8 2 3 751 3 O 

01 

0 . 46664386 1 1 O- 05 

4 

0.39286515960 

01 

0 .46125205820-05 

5 

0 .39200003700 

01 

0 . 45885 1 86420- 05 

6 

0 .39222927 6b0 

01. 

0 .45856361430-05 


MB* 1 FREQUENCY* 3.923430 QORRD/SEC AMPLITUDE* 4.534660-06 NIT* 6 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N»G1+N*N*62*0. 

GAMMA* 1.330000 07 TO* 1.000000-01 .JV= 1.000000 05 JG= 2.1 00000 00 
KP= 2.160000 02 K I = 9.700 000 03 H* 6.00000D 02 K0= 5.758350 03 
Kl= 1.371020 03 K2= 1.000000-01 K3= 3.000000-01 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0.39234340060 

01 

0 . 458466259 1 0- 05 

1 

0.38305731880 

01 

0 . 4333986 1050- 05 

£ 

0 - 3 7 6 7 o 9 1 7 0 1' 

01 

0 .422171 1 1540-05 

3 

U m ■ d * 0 1‘‘ z' 9 c! 1 D 

01 

0.41174951780-05 

4 

U ■ Zt i 0 1 IJ - 1 C w* Z< 5 £i 

01 

0.40524041570-05 

5 

0 . 3 6 9 7 3 3 439 3 D 

01 

0.40373892040-05 

6 

0. 369645 1845B 

01 

0 . 4 0 3 0 7 8 8 8 570- 0 5 

y 

i 

0 - i 0 Uc r fc'5y r D 

01 

0 .40317933200-05 

y 

0 . 37 0 06S43S 0D 

01 

0 .40343951960-05 


NO* 1 FREQUENCY* 3.701570 OORAD/SEC AMPLITUDE* 4.036110-06 NIT*. 8 


Figure 1-20a. 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+H*N^G£=0. 
GAMMA- 1.38 GOOD 07 TD= 1 . 0 0 00011-0 1 JV= 1.0 0 GOOD 05 

K P= £.16 0 0 0 D 0 £ K I = 9.70 0 0 0 D 0 3 H = 6 . 0 0 0 0 0 D 0 £ 

K 1= 1.371 OSD 03 K£= 1 . 00 00 00-01 K3= 4.0QQQQD-01 


I TERRI I ON FREQUENCY AMPLITUDE 


37 0 1 567 1 7 1 Ii 


0 . 4 036 1 08745D- 05 


1 

0 . 358£3 ! 65870 

01 

0. 

37934£64£6D- 05 

c! 

0.34965581680 

01 

0 . 

3589545£850- 05 

o 

0.34355119730 

01 

0 - 

34477807170-05 

4 

0.339893 00380 

01 

0 . 

3 3 5 4 3,7 £ £ 3 6 0 - 0 5 

er 

j 

0.33981930570 

01 

0 . 

33588979960-05 

6 

0.3398 1473860 

01 

0 . 

334 063 14740- 05 

7 

0.3395177657D 

01 

0 . 

3 3 3 9 6 913790- 0 5 


•J6= £.10 GOOD 00 
K0» 5.75335D 03 


HD= 1 FREQUENCY- 3.39539D OORRD-'SEC AMPLITUDE- 3.33865L-06 NIT= 7 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N*N*G£=0 . 

GAMMA- 1.3S000D 07 TO= 1 . 00000D-01 JV- 1 . 0 0 0000 05 

K p = £.16 0 0 0 D 0 £ K I = 9.70 0 0 0 D 0 3 H = 6 . 0 0 0 0 0 D 0 £ 

Kl= 1.37102D 03 K£= 1.000000-01 K3- 5.000000-01 

ITERATION FREQUENCY AMPLITUDE 

0 0 . 3 3 9 5 3 S 5 £ 6 3 0 0 1 0 . 3338654 6740- 05 

1 0 .3209 1 1 35£ OD 01 0 .£976£ 1 49760- 05 

£ 0 . 3 0 0 S 6149350 0 1 0 .£56 03 1 38650- 0 5 

3 0.11 354 1 33c 1 Ii 0 1 - 0 . 84 973759 04 D- 0 6 

1 1.185410 00 -8.497380-07 

4 0.1 0 0 0 0 0 0 0 9 8 0 0 1 0 . 9 9 9 9 9 9 9 0 3 8 0 0 0 

5 0.10 0 0 0 0 0 0 0 OO 01 0.1 0 0 0 0 0 0 0 0 0 0 0 1 

ND= 1 FREQUENCY- 1.000000 OORAD^SEC AMPLITUDE- l.QOQOQD 00 NIT- 5 


Jb= £.100000 00 
K0- 5.758350 03 


Figure 1 -20b . 
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LST SYSTEM-NUMERICAL SOLUTION OF 1 +N*G 1 +N*N^G2= 0 . 
GAMMA® 1 .38 00 Oil 07 TD= 1.0 000 OB- 01 JV= l.OOOOOD 05 

K P= 2.16 0 0 0 B 0 2 K I = 9.70 0 0 0 D 0 3 H = 6 . 0 0 0 0 0 B 0 2 

Kl = 1 .371 02 B 03 K2® 1.00000D-01 K3= 1.00000B-G1 


iteration 

0 

1 


FREQUENCY 
0 . 1 S 0 0 0 0 0 0 0 0 B 
0.132993527311 
0.135465 06 33 B 
0 . 13633794 03 B 
0.1 375561 283B 
0 . 1 3792 5 3743 D 


01 

01 

01 

01 

01 

01 

01 


ND= 1 


FREQUENCY® 1.330 0SB OOF: Alt. •'SEC 


AMPLITUDE 
6 o 0 0 0 0 0 0 0 0 1' - 0 6 

5796 3 224 3 4 D - U6 
537191657211-06 
5 9 3 0 4 3 6 3 5 4 B - 0 6 
596324 t_t 7 4 3 D — U 6 
5973 3 64 735 D - U6 
59786 095 1 2D- 06 


Jh= 2. 10 00 OB 00 
K0= 5 .7583 5 B 03 


AMPLITUDE® 5 .9339911-01 


NIT® 6 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N*N^G2=0. 

GAMMA = 1 .38 00 OB 07 TD= 1 . 0 00 0 OB- 01 JV= l.OOOOOD 05 

KP= 2. 16 00 OB 02 KI= 9. 70 0 0 OB 03 H® 6.00 000B 02 

K 1 = 1 . 37 1 02B 03 K2 = 1 .00 0 0 OB- 0 1 K3= 2 . 0 0 0 0 OB- 0 1 


J6= 2.1 000 OB 00 
K0= 5.75335D 03 


ITERATION 


FREQUENCY 


AMPLITUDE 


0 

0. 1 330 0779 17D 

01 

0 

.5933992937B-G6 

1 

O.l 98941 2767D 

01 

0 

. 672325 1343B- 06 

c 

0 . 1 '5*6 1 c4 1 bbcD 

01 

0 

. 7 1 36786352B- 06 


0.1931 174 4 32 D 

01 

0 

.74528 2 3 2 4 6 B - 0 6 

4 

0.1 979633 157D 

01 

IJ 

. 77 1 5334579D- 06 

5 

0 . 199 00 7 372 OD 

01 

o 

. 77 47937864D- ij6' 

£ 

0. 1 996 13 0359 B 

01 

0 

.7774268861B-06 

“7 

0 . 1999 099 796D 

01 

0 

. 7783768 1 35B- 06 

8 

0 . 199920766 IB 

01 

0 

. r 7368y72b'5B- 06 


NB= 1 


FREQUENCY® 1 .99932B OORAB.-'SEC 


AMPLITUDE® 7.78826B-07 


NIT® 


— o 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N*6 1 +N*N*62= 0 . 
GAMMA® 1 .38 00 OB 07 TD= 1.00000D-01 JV= l.OOOOOD 05 

KP= 2.16 0 0 OB 02 KI = 9. 7 000 OB 03 H= 6. 00 00 OB 02 

Kl= 1.37102E 03 K2- 1.00000D-01 K3= 3.00000D-01 


JG= 2. 1000 OB 00 
K0= 5 . 75335B 03 


ITERATION 


FREQUENCY 


AMPLITUDE 


0 

0 . 1 99932387611 

01 

U • r • * J b lJ b iJ 5 C'4 I) U b 

1 

0 . 2 062672 759 B 

01 

0 .87343191 04 B- 06 

cl 

0.21 09 1 22 09 1 D 

01 

0 . 9 3 6824 9 6 U 1 B - U 6 

O 

0.21 1536 739 3D 

01 

0 • 1 0 0 7 9 2 0 0 '39 D - 0 5 

4 

U .cl i'cc'5 j'cc 1 B 

01 

0 . 1 0 U9 IJ U68 04 D— 05 

er 

0.21 456656 78B 

01 

0.101 6 0 09 07 3D- 05 

6 

0 . 2 1 52649862B 

0 1 

0 . 1 02 0 05 1 049 D- 05 

*7 

0.2 15501307 ID 

01 

0 . 1 020575384 D- 05 

y 

0.21551231 39B 

01 

0 . 1 02 1 4 386*89 D— 05 


NB= 1 


FREQUENCY® 2.15534B 0 0 RAD-" SEC AMPLITUDE® 1.02156D-06 


NIT® 8 


Figure 1 -2 1 a . 
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LST SYSTEM-NUMERICAL SOLUTION OF 1 + N*G 1 +N*N^G2= 0 . 


GAMMA* 

1 .33 

0 

o 

OB 

07 TO* 

1 . 

0 0 

0 

0 

KP 

— C. • 

1 6 0 0 0 

B 


02 

K I = 


70 

0 

0 

K1 

= 1 . 

»■’ 1 U cl 

D 


0 3 

K2= 

1 . 

0 0 

0 

o 


ITERATION 




FREQUENCY 






0 


o 

■ cl 

1 553392 

6 1 D 

0 

1 




1 


0 

■ c. 

2442251 

91 D 

0 

1 






o 

• c. 

03537 08 

35D 

0 

1 




O 


n 

• c. 

1802474 

8 OD 

0 

i 




4 


0 

w c 

d t 1 1 c c. 6 

14B 

0 

1 




tzr 

-> 


0 

• *2. 

3239593 

9 OD 

o 

1 




6 


o 

• C 

36 1 0756 

24 D 

0 

1 




"7 

i 


o 

* a 

3771848 

15B 

o 

1 




8 


0 

- cl 

3824993 

2SB 

o 

1 




3 


0 

• c 

383 05 06 

92B 

0 

1 



10 0.23837717141) 01 


0D-01 JV= 1.0 000 Oil 05 

0 D 03 H= 6 . 0 0 0 0 0 0 02 

on- o i K 3 = 4.0000 on-o 1 

AMP L ITU PE 
0 . 1 02 1 562747B- 05 
0.11 53 9885 1 7D— U 5 
0.1261 87 0534 Ii- 05 
0 . 122574526 IB- 05 
0 . 1 29 093631 9 B- 05 
0 .134 -3 0 3 915 3 D - 0 5 
0.1 373743 172B- 05 
0 . 1333691 12 OB- 05 
0 . 1 392 1 66 0 03 B- 05 
0 . 1 393 8 93 795 B- 05 
0 .13947 084 94 B- 05 


.Jb= £. 1000 OB 00 
KO= 5 .753 35 B 03 


HD* 1 FREQUENCE 2.38425B 0 GRAB -'SEC AMPLITUDE* 1 .39567D-06 NIT* 10 


LST 

SYSTEM- 

NUMERICAL SOLUTION OF 


1 +N+G 1 +N*N*’G2= 0 . 

GAMMA* 1 .38 

00 

OB 07 TO* 1 . 

00 0 0 OD- 

01 JV= 1.00000D 

KP* 

2.1600 OB 

02 K I = 9 . 

7000 OD 

o 

3 H* 6.00 0 0 OD 

Kl = 

1 .37102 

B 

03 K2= 1 . 

00 00 OB- 

01 K3= 5.0000 OD- 

ITERATION 


FREQUENCY 



AMPLITUDE 


0 

0 

. 2384254 0 08B 

01 

0 

. 1 39566676 OD- 05 


1 

0 

.25:j560yySSD 

01 

0 

. 16271 9281 7D- 05 


d 

0 

.263 01 8974 ID 

01 

o 

. 1 90988457 1 D- 05 


o 

0 

.28 063669 02B 

01 

0 

. 2 1 79832984 B- 05 


4 

0 

.297931 87 18B 

01 

o 

. 3 4 6 3 6 5 4 £ 1 3 D — L J 5 


5 

0 

.293701 0575D 

01 

0 

. 25 1 26 98 384 D- 05 


6 

0 

.253933775 IB 

01 

0 

. 17739529831'- 05 


*7 

1 

o 

. cl 6 3 cl 3 6576 3 B 

01 

0 

. 1941443335D-05 


8 

0 

.279 0768 007B 

01 

0 

.21 0567793911-05 


9 

0 

. 2 8 1 6 c‘ 3 U 9 8 6 B 

01 

0 

.2292847674H-05 


1 0 

0 

.•3214 144945B 

01 

Ij 

p 7 3 855657 6 OD— 05 


1 1 

o 

.30 1 5994 1 59D 

01 

0 

. £575 097 1 96Ii— U5 


12 

Ij 

. 1819751304D 

01 

0 

. 3275667 48 IB— u6 


13 

0 

. 18441 00614B 

01 

o 

. 5327 1 49873B- 06 


14 

o 

. 1 997997856D 

01 

Ij 

. 78451 7 0266B- 06 


15 

o 

.2160355299B 

01 

n 

. 1 0 3 6 8 8 5 4 0 7 B - 0 5 


16 

o 

. 63 1 6487 477B 

01 

o 

. 1 2 3 165 8 9 3 0 B - 0 5 


17 

0 

.24724571 63D 

01 

0 

. 152 0 15031 4D- 05 


18 

0 

. 25 034 93328 B 

01 

0 

. 1 3 0 7 0 3 4 1 3 0 B - 0 5 


19 

o 

.26433651 18B 

01 

0 

. 195391 048 OB - 05 


20 

0 

.281 9 03697 6B 

01 

o 

.2217526558H-05 


JG= 2.10 0 0 OB 00 
K0= 5.75835B 03 


NB= 1 FREQUENCY* 2.81909B OORAB-'SEC 


AMPLITUDE* 2.21753B-06 HIT* 20 


Figure 1 -2 1 b . 
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1ST SYSTEM-NUMERICAL S 
GAMMA* 1 .3300011 07 TD= 
KP= £.1600 OH 02 K I * 
Kl = 1.37102D 03 K£= 


GLUT I ON OF 1+N^G1+N*N*G£=0. 
1.00000D-01 .JV= l.OOOOOD 05 

9.70000D 03 H= 6.00000D 02 

i .oo o o on oo K3= l.oooo on - 03 


ITERATION 

FREQUENCY 

0 

0 .43000000 OOD 

1 

0. 43 0528 13 86 D 

2 

0 .43 05283350D 

3 

0 . 43 087443? OD 

4 

U .43036 3 4 3 5 3 D 


AMPLITUDE 

0 1 0 . 520 0 0 0 0 0 0 0 D - 0 5 
01 0 . 52233 039 0 1 D— 05 
01 0 . 5222453439D- 05 
01 0 .52457437230-05 
01 0.5243241 015D- 05 


JS* 2. 10 GOOD 00 
K0= 5.75835D 03 


MB* 1 FREQUENCY* 4.30865D OORAD'SEC AMPLITUDE= 5.2418511-06 NIT= 4 


LST SYSTEM-NUMERICAL SGLUTICN OF 1 +N*G 1 +N*N*G2=0 . 

bAMMA= 1 .33 00 On 07 T0= l.OOOOOD-Ol JV= l.OOOOOD 05 JG= 2.10000D 00 
KP= 2.16GG0B 02 K I = 9.700QOD 03 H= 6.00000D 02 K0= 5.75835D 03 
Kl= 1.37102D 03 K2= l.OOOOOD 00 K3= 1.00000D-02 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0.43 086451 93D 

01 

0 . 524 1 85 0££9D- 05 

1 

0 .429859 05 14D 

01 

0 . 522 0253 0 1 3D- 05 

. lL 

0 . 428729 4 964D 

01 

0. 5 1 94 4 726 04 D— 05 

3 

0.4 29191 35 05D 

01 

0 . 5 1 94433798D— 05 

4 

0 . 429 1 7624 1 ?D 

01 

0.51 9386774 9D— 05 


NB= 1 FREQUENCY* 4.29161D OORAD^SEC AMPLITUDE* 5.19344D-06 HIT* 4 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1+N*N*G2=0 . 

GAMMA* 1 .33000D 07 TO* 1.0G00GD-G1 JV= l.OOOOOD 05 JG= 2.10000D 00 
KP= 2.16000D 02 KI= 9.70000H 03 H= 6.00000D 02 K0= 5.75835B 03 
Kl= 1.37102D 03 K2= l.OOOOOD 00 K3= l.OOOOOD- 01 


ITERATION 

FREQUENCY 


AMPLITUDE 

. 0 

0 .429 16 0758 OD 

01 

0 . 5 1 934 4 06 9 9D- 05 

1 

Li • 4 c! d. y 6 4 1 < 4 y J B 

01 

0 . 5 0262 0 0492D- 05 

o 

L. 

0 .417709477011 

01 

0 . 4895733868D- 05 


0.41 49 095 0 04 0 

01 

0 .48 14 07154 OD- 05 

4 

0.4 131 739544 D 

01 

0 .47646 078? OD- 05 

5 

0.41262777 ?Sn 

01 

0 .47468227550-05 

h, 

0.41281 02 0 03 H 

ill 

0 .4743764 169D- 05 


ND= 1 FREQUENCY* 4.12S99D OOPAD-'SEC AMPLITUDE* 4 .742690- 06 NIT* 


Figure 1-2 2a. 
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LST SYSTEM-NUMERICAL S 
GAMMA= 1 .38000H 07 TO= 
KP= 2.16000B 0£ KI= 
K 1= 1 .37 1 0£B 03 K£= 


GLUT I DM QF 1+N*G1+N*N^G£=0 . 


1 . 00 00 OB- 01 
9. 7 000 OB 03 
1.00 GOOD 00 


JV= 1.00 0 0 0 D 05 
H= 6.00000D 02 
K3= £ . 0000011-01 


JG= 2.10000B 00 
K0= 5.75835B 03 


ITERATION 

0 

1 

£ 

3 

4 

5 

6 


FREQUENCY 
0 . 4 1 £ y 9 3 6601 D 01 
0 .403833831 SB 01 
0 .397707 8 0 4 £ B 0 1 
0 .3938 14 0399 D 01 
0.391 3 0 5 £ 8 3 5 B 0 1 
0. 39094 17445D 01 
0 . 3 9 0 3 686365 D 01 
0.391 0473948B 01 


AMPLITUDE 
0.4742691 1 03B-05 
0 . 4534458733D- 05 
0 .43699 1 6734D- 05 
0 . 426456592 6' D — U 5 
0 . 4 1 98877967B— 05 
0 .4183451 036D-05 
0 .4176817051D-05 
0 .4 17557537 0B- 05 


ND= 1 FREQUENCY= 3.91085B OORAD^SEC flMPLITUBE= 4.17559D-06 NIT= 7 


LST 3 YSTEM-NUMER I CAL SOLUTION OF 1+N*G1+N^N*G£=0 . 

GAMMA= 1 .38 00 OB 07 TD= 1.00000B-01 JV= 1.00000B 05 JG= 2.10000B 00 
KF'= £.16 0 0 OB 02 KI= 9. 7000 OB 03 H= 6. 00 GOOD 02 K0= 5.75835B 03 
Kl= 1.37102B 03 K£= 1.00000D 00 K3= 3.00000B-01 


ITERATION 

0 

0 . 

FREQUENCY 
391 0350693B 

01 

AMPLITUBE 
0 .417553635SB-05 

1 

0 . 

7 9 6 U £ 3 0 9 6 B 

01 

0.3926734444B-05 

■- 

L_ 

0. 

37 128 03 03 5B 

01 

0 .372054 1479B- 05 

O 

0 - 

3654 1 7353 4 B 

01 

0 . 3577 6 6594 IB — 05 

4 

0 . 

36 1 3985577B 

01 

0 . 3485 1 87376B- 05 

5 

0 . 

361942 049 7B 

01 

0 . 3483 1 883 £9B— U5 

»s 

0. 

361 1493 035B 

01 

0 . 3468979466B— 05 

7 

o.. 

36 15626 7 43 B 

01 

0.3468727258D-03 

p 

0 • 

3615619356B 

01 

0 .34674 1 56 19B- 05 


ND= 1 FREQUENCY** 3.61565B OORAB/SEC AMPLITUDE** 3.46667B-06 NIT= 3 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N^G 1 +N*N*G£= 0 . 

GAMMA= 1 .380000 07 TO= 1.00000D-01 .JV= 1.00000B 05 JG= 2.10000B 00 
KP= £.16 0 0 OB 0£ KI= 9.70000D 03 H= 6.00000D 02 K0= 5.75835D 03 
Kl= 1 .371 02 B 03 K£= 1.00 00 OB 00 K3= 4. 0000 OB- 01 


ITERATION FREQUENCY 

0 0 . 36 1 5649 056B 01 

1 0 .34 3 429794 9 D 0 1 

£ 0.384 06776.4 9 B 01 

3 0.1 3 09721 258D 01 

1 1.30972B 00 - 1 . 57257B-07 

IHC900I EXECUTION TERMINATING DUE TO ERROR COUNT FOP ERROR NUMBER 


AMPLITUDE 
0 . 3 4 6 6 6 ? 8 4 9 6 B — U 5 
0 . 3 0 9 3 0 8 9 0 1 £ D - 0 5 
0 . £ 6 699 8 0 1 8 £ B - 0 5 
-0 . 1 572573 1 93D- 06 


£61 


Figure 1 -22b . 
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L ST SYSTEM-NUMERICAL SOLUTION OF l+N*Gl+N*N*G2=0 . 
GAM fc, A» 1 . 3 80 COD C7 to = l.OCOCOD-Ol J V* 1. OOOOOD 05 


JG= 2.100000 00 


KP= 2 . 160 OOD 02 K I = 9.700000 03 H= 6.00000D 02 

K 1- 1.371020 03 K2= 1. OOOOOD 00 K3= 1.00000D-03 


K0= 5.758350 03 


ITERATION 

0 

1 


FREQUENCY 
0.21 OOOOOOOOD 01 
0.21058192930 01 


AMPLI TUOc 
0.5 8000000000-06 
0.57644082180-06 


2 0.21 106834450 01 

3 0.21134006370 01 

4 0.2 1 134585760 01 


0.57779394370-06 

0.57848853980-06 

0.57875519580-06 


N0= 1 FREQUENCY 55 2.113520 OORAD/SBC AMPL ITUDE= 5.788290-07 NIT= 4 


LST SYSTEM-NUMERICAL SOLUTION OF 1 + N*GH-N*N*G2 = Q. 
GAMMA= 1 . 380 f OO C7 T 0 = l.QnoOOO-Cl JV= 1.00C09D 05 


JG= 2.100000 00 


KP= 2. 160000 02 K 1= 9.70000D 03 H = 6.000000 02 

Kl= 1.371020 03 K2= l.QOOQOD 00 K3= 1.00000D-02 


K0= 5.75835D 03 


ITERATION 

0 

1 


FREQUENCY 
0.21135223630 01 
0.21180918 760 01 


AMPLITUDE 
0.57882942720-06 
0.5 850930 129D- 06 


2 0.21 1989521QD 01 

3 0.21209051680 01 

4 0. 21 21552444D 01 


5 0.21219679300 Oi 

6 0. 2 1 2226 1884D 01 


C. 58873250600-06 
0.5907279733D-06 
C. 5 5186241360-06 


0. 5526023104D-06 
0. 5530924560 D- 06 


N0= 1 FREQUENCY 55 2.12248D OORAO/SEC AMPLITUDE 55 5. 93446D-07 N I T= 6 

LST SYSTEM-NUMERICAL SOLUTION OF 
GAMMA= 1.330000 07 T0= 1.000000- 
K P= 2. 160 COD G2 K I = 9.700000 

1+NtG1+N*N*G2=0. 

01 J V= 1 . OOOOOD 05 JG= 2.100000 

03 H - 6.000000 02 K0= 5.75835D 

00 

03 

K 1= 1.371C2D 03 K2= 1.000000 

00 K3= 1 . OOOOOD- 01 


ITERATION FREQUENCY 

AMPLI TUDE 


0 0.2 l 22476046D 01 

1 0.21673833730 Cl 

2 0.21 9582 8C91D 01 

0.5 934460004D-06 
0.666285 8378 D- 06 
C. 7 11400 74880-06 


3 0.22128323420 01 

4 0.22248267900 01 

5 0.22276668260 01 

0.73721802100-06 
0.74964577550-06 
0.759949 10370-06 



7 0.22267340260 01 

8 0.2230362791D 01 


0.22313775530 01 


0.77067031190-06 
C. 7 71421 38 800-06 


C . 7 7 Qd 59 62860-06 


ND= 1 FREQUENCY = 2.231410 OORAO/SEC AMPLITUDES 7. 713320-07 N I T= 9 


Figure 1-23a. 
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LST SYSTEM-NUMFR ICAL SO LOT I DM OF 1 + N* G 1 + N »N* 02 3 0 . ' 

G4M^a= 1.38000D 07 T0 = I ,000000- OL J V= l.OOOOOD 05 JG= 2. 100000. 00 

KP= 2.16000D 02 KI= 9.700000 03 H= 6.00000D 02 K0= 5.758350 03 

Kl= i. 371020 03 K 2 3 l.OOOOOD 00 K3= 2.00000D-Q1 

iteration frequency amplitude 

0 0.2 2314055190 01 0.7 713 3 2 3 L 8 2D- 06 

1 0 . 2 2 9 1 8 9 65 9 50 01 C . 8 70378 57 240-06 


2 

0.23332686360 

01 

C .93602606470-06 


3 

0.2 3674924000 

01 

C. 97474662960-06 


4 

0.23 758376840 

OL 

0.1C01 50 89090-05 


5 

0.23 79792587D 

01 

0. 10134590340-05 


6 

0.23 822084950 

01 

0.1C21661878D-05 


7 

0 . 2302672 134 D 

01 

0. 1C260210140-05 


8 

0.23 84497 1350 

01 

0.1C265748260-05 


9 

0.23847699590 

01 

0. 1C279Q21240-05 


10 

0.23852289870 

01 

0.10283023410-05 



ND= 1 FREQUENCY 3 2.38548D OORAD/SEC AMPLITUDE 3 1.028880-06 NIT= 10 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N*N*G2=0. 


GAMMA* 1.38000D 07 TQ = l.OOOOOD-Ol JV= l.OOOOOD 05 JG= 2.100000 00 


K P 3 2. 1 60 OOD 
Kl= 1.371020 

C2 K I = 9. 

03 K 2= l.i 

71. GOOD 
DO 0 OOD 

03 H- 6 . 0000 0 D 02 

00 K 3= 3. OOOOOD— 01 

KQ= 5.75835D 03 

ITERATION 

FREQUENCY 


AMPLI TUDE 


0 

0.23854834610 

01 

0.1C2888 1516D-05 


1 

0.24678013500 

OL 

0.1 1650479630-05 


2 

0.25398528850 

01 

0.1 266 34 8 36 2D— 05 


3 

0.25578450900 

01 

0. 1 292 37 4 39 8D— 05 


4 

0.25776367730 

01 

0. 134420 1029D-05 


5 

0.2 5 8826 4 309D 

01 

0. 1369898918D-05 


6 

0.259401 1438D 

01 

0. 13843215460-05 


7 

0.25995488080 

01 

0. 1395250L81D-05 


8 

0.26029329 43D 

01 

0. 1402 7“5 89 71 0-05 


9 

0 . 260 56 5 33 5 ID 

01 

0. 14103999390-05 


10 

0.26 05424 268D 

01 

0. 1425499818D-05 


11 

0.26 124580090 

01 

0.14290249370-05 


12 

0.26151963730 

01 

0. 14293534980-05 


13 

0.26 15238207D 

01 

0. 14301339780-05 



ND= 1 FREQUENCY 3 2.615340 OORAO/SEC AMPLITUDE 3 1. 430250-06 NI T= 13 


A 


Figure 1 -23b . 
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LST SYSTEM-NUMERICAL SOLUTION OF l+N*Gl+N*N*G2«0. 


GA^ MA= I . 380000 07 to- I.QOOOQO-ql J V= l.OOOOOD 05 JG= 2.100000 OC 

W* D 'J 1 /! A O ^ i y T -r j— v /-i r-s r-s — . _ . . 7~" — ^ — _ — _ Z I " - • — — • ~— - — 1 '* 


KP= 2.160000 02 Kl* 9. 

700000 

03 H = 6 .000000 02 

K0= 5.758350 02 

Kl= 1.371020 03 K 2= 1 . 

000000 

00 K 3= 4. 000000-01 


iteration 

FREQUENCY 


AMPLI TUDE 


0 

0.26153383290 

01 

0. 14302522270-05 


1 

0.2749145 1770 

01 

0. 1 6576896760-05 


2 

0. 27155CC 606D 

01 

0. 1 8331919630-05 


3 

0.285686 70450 

01 

0. 19921476640-05 


4 

0.30181244280 

01 

0. 22431135990-05 


5 

0.3236273 1020 

01 

C. 26300592730-05 


6 

0.29855021500 

01 

0.22525799410-05 


7 

0.31768331180 

01 

0.25675390120-05 


8 

0.19 141 123C9D 

0 1 

0.42841346750-07 


9 

0.20833868460 

01 

0.643026 84600-06 


10 

0.22882432030 

01 

0 . S 35 74 6 1 00 8 D -06 


11 

0. 24575102290 

01 

0. 1 1954391140-05 


12 

0.25991500320 

01 

0. 14232 643080-05 


13 

0.27197696170 

01 

0.16358287280-05 


14 

0.28244343390 

01 

0.1 E2099 9738 0-05 


15 

0.29220335750 

01 

0.19938006230-05 


16 

0.28804182890 

01 

0.2C716253490-05 


17 

0. 30 153749000 

01 

0.22590807780-05 


18 

0. 32966725560 

01 

0.27652305310-05 


19 

0. 30 0 5244 254 D 

01 

0.22179562590-05 


20 

0.32077926C70 

01 

0. 25743091210-05 



ND= 1 FREQUENCY* 3.20779D OORAO/SEC AMPLITUOE= 2. 574310-06 NIT* 20 


LST SYSTEM-NUMERICAL SOLUTION OF 1 + N*Gl + N*N* G2 = 0 . 

GAMMA= I.33CCCD 07 TT= l.OOOOOP-QI J V= l.OOOOCD C5 JG= 2.1000CD 00 

KP= 2.160000 02 KI= 9.700000 03 H- 6.000000 02 K0= 5 ■ 75835D 03~ 
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Figure 1 -23c . 
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Figure l-24a. 
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Figure 1 -24b . 



57 
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Figure 1 -25a 
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Figure 1 -25b . 
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2. Prediction by Numerical Methods of Self-Sustained Oscillations 
in a Two-Axis Model of the Nonlinear Sampled-Data LST System 

2-1 . Introduction 

It has been demonstrated in [2] that the stability equations as 
a result of the describing function analysis may be solved by a 
numerical -iterative technique. In Chapter 1, the numerical technique 
is applied to the prediction of self-sustained oscillations in a two- 
axis model of the nonlinear LST system with continuous-data. In this 
chapter the same numerical technique will be extended to the prediction 
of self-sustained oscillations in a two-axis model of the nonlinear 
LST system with sampled-data . 

In the two-axis model of the LST system, the system contains two 
nonlinearities, and the general form of the stability equation in the 
sampled-data case is 

1 + G B (T,n)N(A,n,<f>) + G c (T,n)N 2 (A,n,c|>) = 0 (2-1) 

where G B (T,n) and G^(T,n) are transfer functions which describe the 
linear elements of the coupled systems, and N(A,n,<fp) represents the 
discrete describing function of the CMG nonlinearity. The amplitude of 
the sinusoidal signal at the input of the nonlinearities is denoted by 
A; T is the sampling period in second, and n is the integer which relates 
the period of the sustained oscillation to the sampling period; 4> is the 
phase angle of the input sinusoid. Therefore, 
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Tj, = nT = period of sustained oscillation 

. or 

2tt 

0 ) c = = frequency of sustained oscillation in rad/sec 

2-2. The Sampled-Data Two-Axis LST System Model and Its Stability 

Equation 

Figure 2-1 illustrates the signal flow graph representation of 
the sampled-data single-axis LST system. The nonlinear friction of 
the CMG is modeled by the branch with the gain N, where strictly, N 
denotes the discrete describing function N(A,n,<|>). The sample-and-hold 
at the input of the nonlinearity is added for analytical purposes, so 
that the overall system may be treated by the 2 -transform method. 

Since the actual system does not have the sample-and-hold preceding N, 
the system configuration of Fig. 2-1 should be considered as an approxi- 
mation to the sampled-data single-axis LST system. 

Figure 2-2 shows the two-axis model of the LST system, with the 
coupling made at the output stages. The coupling coefficients and 
K 3 may be varied for various degrees of coupling between the two single- 
axis systems. 

For the study of self-sustained oscillations, we set x-j = X 2 = 0- 
From the signal flow graph of Fig. 2-2, the following cause-and-effect 
relationships are written for the variables at the outputs of the four 
samplers: 

E* = -A*E* + A*N*X* + A*E* - A*N*X* (2-4) 

E* = -A*E* + A*N*X* + A*E* - A*N*X* 


( 2 - 2 ) 

(2-3) 


(2-5) 



K s + K t 

6-| (s) = K q + K ]S 

V s) ' J r s 

" Ks + K, 
P I 

U G 

G 7 (5) = jV 

V 

1 n' TS 

«lio< s >" s 


Figure 2-1. The single-axis 

sampled-data LST system model. 
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X* = -B^N*X^ + B*E* 
X* = -B*N*X* + B*E* 

where 


A-, = 


G ho G l G 2 G 3 G 4 G 6 HA 1 


G. G-. G^G^G.Gj-HKgA-, 

_ ho 1 2 3 4 6 2 1 _ « r „ 

A 2 A A l b 4 K 2 

o 


A. = 


G ho G l G 4 G 6 HA 1 


6, G-, G / ,G^HK 9 A 1 
a - ho 1 4 6 2 1 r K 

A 4 A A 3 G 4 K 2 

o 


B, = 


G ho G 6 A 2 

SA_ 


Bo = 


G ho G 2 G 3 G 6 A 2 _ P r r 
sA^ B 1 G 2 G 3 


An = 


1 + G 3 G 6 


A_ = 


1 ♦ 2G 3 G 6 - K |g| + G| G Z 6 - Zl/foKj - g|g|g|k| 


A 2 = 1 + G 3 G 6 - G 3 G 6 K 2 G 4 ‘ K 2 G 4 


( 2 - 6 ) 

(2-7) 

(2-8) 

(2-9) 

( 2 - 10 ) 

( 2 - 11 ) 

( 2 - 12 ) 

(2-13) 

(2-14) 

(2-15) 

(2-16) 


In the last nine equations all functions are functions of s. The 
sampled signal flow graph of the system is drawn in Fig. 2-3 using 
Eqs. (2-4) through (2-7). 

The determinant of the sampled signal flow graph is obtained from 
Fig. 2-3 by use of Mason's gain formula. The expression is given in 
terms of the z-transform functions: 
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A(z) = [1 + 2A-| (z) - A 2 (z) + A^(z)] + 2N(A,n,<t)[B 1 (z) - B 2 (z)A 3 (z) 

+ A 2 (z)B 2 (z)A 4 (z) + 2A 1 (z)B ] (z) - A ] (z)A 3 (z)B 2 (z) + A^zjB^z) 

- A 2 (z)B ] (z)] + N 2 (A 5 n,(J))[B^(z) - A 4 (z)B 2 (z) - 2B ] (z)B 2 (z)A 3 (z) 

+ A 3 (z)B 2 (z) + 2A 2 (z)A 4 (z)B- 1 (z)B 2 (z) + 2A 1 (z)B 2 (z) 

- 2A 1 (z)A 3 (z)B- | (z)B 2 (z) + A 2 (z)B 2 (z) - A 2 (z)B 2 (z)] 


The components of A(z) are determined as follows: 


A-, (z) 


0 - Z_l )3 


K I HJ v (K o tK 1 s)(J G s^K p S + K I )(s^ + K 3 s + Kg) 


SA 


op 


A 2 (z) 


= (1 - 


z_1 >* 


K,K 2 H(K 0 + + K s + Kj ) 


op 


A 3 (Z) 


= (1 - z' 




HJ v {K q + K-jS) { J^s 2 + K S + Kj)(s* + K 3 s + Kg) 


op 


V z) 


B-j (z) 


B 2 (z) 


(1 - z-')J 


HK 


2 (K q + K ]S )(J g s^ + K s + Kj) 


= (1 - z" 


- (1 - Z 


5 a 




op 


(V 2 + K P s+ Kj )[J^(s 2 + K 3 s + K 2 ) 2 - Kg] 


SA 


op 


W 2 + K p S + K I^Jy(s 2 + K 3 S + k/ 


K 2 ] 


S 2 A 


Op 


A cp = [J G s4 + ^pV* + < K p * 2 K I J G )S ‘ + ZK p K l S + + V 


.3 . , „Z 


2i r t 2 / 2 


+ Kg ) 2 - K 2 ] 


(2-17) 

(2-18) 

(2-19) 

( 2 - 20 ) 

( 2 - 21 ) 

( 2 - 22 ) 

(2-23) 


(2-24) 
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The expressions for the discrete describing function of the CMG 
non! inearity, N(A,n,<|>), may be found in [2]. In the present case, 
since the two axes are considered to be identical, and the couplings 
are symmetric, it is assumed that the amplitudes of the input signals 
at the two nonlinearities are identical. 

The stability equation is obtained by setting A(z) of Eq. (2-17) 
to zero. It is apparent that A(z) = 0 is of the form: 

A(z) = 1 + G a (z) + Gg(z)N(A, + G c (z)N 2 (A,n,4>) = 0 (2-25) 

which can be written as 


1 + Gg(T ,n)N(A,n,(J)) + G c (T,n)N 2 {A,n,<j>) = 0 


(2-26) 


by dividing both sides of the equation by 1 + G^(z), and setting 


z 


e J27r/n 


(2-27) 


Equation (2-26) shows that in the sampled-data system, there are 
four variable parameters in the stability equation in T, n, A, and <f>, 
whereas in the continuous-data system there are only two variables in 
A and to. 


2-3. Exact Solution of the Stability Equation of the Two-Axis Sampled- 
Data LST System by Numerical -Iterative Techniques 

The numerical-iterative method described in Chapter 1 is now used 

b 

to solve the stability equation of the two-axis sampled-data LST system, 


Eq. (2-26). 

In Eq. (2-26) we define 
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Gg(T»n) - G^-j + jGj-j 

(2-28) 

G c (T,n) = G r2 + jGj 2 

(2-29) 

N(A,n,4>) = N r1 + jNj-j 

(2-30) 

N 2 (A,n,(j)) = N r2 + jN I2 

(2-31) 


where G R1 , G^, G r2 , Gj 2 , N r -j , N^, N R2# and Nj 2 are all real quantities. 

When Eqs. (2-28) through (2-31) are substituted into Eq. (2-26), 
and after grouping the real and imaginary parts, we have, 

A = A R + jAj = 0 (2-32) 

where 

A R = 1 + G R1 N R1 " G I1 N I1 + G R2 N R2 ‘ G I2 N I2 " 0 (2-33) 

h ' G I1 N R1 + G R1 N I1 + G I2 N R2 + G R2 N I2 ‘ 0 < 2 ‘ 34 > 

These two expressions represent two equations in four variables 
in A, T, <p, and n. It is necessary to assign values to two of these 
variables and solve for the remaining two using the numerical -iterative 
method. In general, it is practical to assign values to T and n, and 
Eqs. (2-33) and (2-34) are solved to yield A and <|>. Alternatively, we 

may assign values to n and <j>, and solve for A and T. The latter case 

/ 

has been carried out in the current study. 

As in Chapter 1, we define 

T 
A 


x 


(2-35) 
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F = 



Then, Eqs. (2-33) and (2-34) can be written as 


(2-36) 


F(x) = 0 (2-37) 

The Newton-type quadratically convergent numerical method described 
in Chapter 1 can now be directly applied to this two-variable system. 

The following system parameters are used: 


"^GFO ft-lb 

H = 600 ft-lb-sec 

J G = 2.1 ft-lb-sec 2 

K = 5758.35 
o 

K-j = 1371.02 

Kp - 216 ft-lb/rad/sec 

Kj - 9700 ft-lb/rad 

J v = 10 5 ft-lb-sec 2 

y = 1 .38 x 10 7 


Three sets of coefficients of coupling are used for this study: 

(1) k 2 = 0.1 k 3 * 0.1 

(2) K 2 = 1.0 K 3 = 0.1 


(3) K 2 = 5.0 


k 3 = 0.1 
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For each integral n, and a fixed value of <}>, a solution of Eq. (2-37) 
is attempted with an initial guess on T and A. A convergent solution for 
T and A together with the assigned n and (j> describes a sustained oscillation 
in the system. One of the characteristics of the numerical method is that 
the solution depends on the initial conditions set by the analyst. As 
shown in the earlier studies, it is quite possible that a system may have 
a multiple number of solutions for the same set of n and <J>. An exhaustive 
search can only be conducted by a systematic scan of all the possible 
combinations of the initial values in the iterative process, unless one 
has general information at the outset on the number and the proximity of 
the solutions. 

Figure 2-4 shows the plots of the two solutions of T versus A for 
K 2 = K 3 = 0,1 ’ ancl for rela tive large values of n. The results are shown 
only for <j> = 0°. When all the values of <f> are used, the solutions of T 
and A trace out a loop in the T versus A domain. For large values of n, 
these loops are very small, as the solutions are relatively insensitive 
to n. As n decreases, these loops become more pronounced, as indicated 
in Fig. 2-5. As n decreases, one of the solutions disappears, as the 
amplitude of A becomes exceedingly small. 

Figure 2-6 gives the computer results for l< 2 = S = 0.1 that 
correspond to one set of solutions for <j> = 0°. Figure 2-7 tabulates the 
other set of results for K 2 = K 3 = 0.1 with <j> = 0°. Figure 2-8 gives the 
tabulations of results for K 2 = = 0.1, n = 4, and variable <j>. 

Figures 2-9 and 2-10 give the results with 1C, = 1.0 and = 0.1. 

Figure 2-11 shows the solutions for K 2 = 5.0 and K 3 = 0.1 with <f> = 0°. 

In this case, the amplitude of oscillation decreases as n decreases, so 
that there is less chance for the low-frequency oscillations to occur. 



73 


The results indicate that harder coupling (larger values of K« and K^) 
tends to reduce the chance of self-sustained oscillations. 



(sec) 


1.38 x 10 


Amplitude ( xlO ) 

Figure 2-4. Amplitude and frequency (27r/nT) of self-sustained oscillations for various sampling 

7 

periods in the discrete-data two-axis LSI system; y - 1.38x10 . 
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LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE, K2= 1.000D-0I* K3= l. 0000-01 
GAMMA* 1 .33 GOOD 07 N= 100 PHI* 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0.33200000000-01 0 .60000000000-06 

1 0.3342££9250D-01 0.59506793070-06 

£ 0 . 3359786 0230- 0 1 0.57479764490-06 

3 0 .33663515850-01 0.566821 07 130-06 

4 0.33695665120-01 0.56253290430-06 

5 0 . 33695372 15D- 0 1 0.56274 0569 ID- 06 

6 0 . 3369532634D- 0 1 0 . 5435 03 029 0D- 06 

7 0.3364 0466280- 0 1 0 .35716994390— 06 

8 0 . 336 7637 1 £3D- 0 1 0 . 5593 1 £54 54 D- 06 

9 0 . 33698558920- 01 0 . 56 09S 058S6D- 06 

10 0.336684 06 04D- 01 0.56294 2754 0D- 06 

11 0. 33693 032230- 01 0 .563 09639 18D- 06 

NO* 1 SAMPLING PERIOD* 3 . 369720- 02SEC AMPLITUDE* 5.63342D-07 NT* 11 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE > K£= 1. 0000-01 > K3= 1.000D-01 
GAMMA* 1.330000 07 N* 96 F'HI= 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0.33697203490-01 0 .563342355 0D- 06 

1 0 .347429 153SD- 01 0 . 5443359999D- 06 

£ 0.35 023754 02D- 01 0 .5470937221D-06 

3 0.351 151 2195D-01 0 .551 29091 1 7D- 06 

4 0.351 5743863D- 0 1 0 . 5536 0397990- 06 

5 0 .35160814550-01 0 .554694361 7D-06 

6 0 .35159674660-01 0 .556299 1373D-06 

7 0.35136354620-01 0 .5566.5054490— 06 

NO* 1 SAMPLING PERIOD* 3 .51S49D-02SEC AMPLITUDE* 5.56679D-07 NT* 7 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE. K£= 1 . OOOD-Ol > K3= 1.000D-01 
GAMMA* 1 .330000 07 N= 92 PHI* 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 . 35 1 6492 122D- 0 1 0 .55667903310-06 

1 0.36326837390—01 0 .536.64324870— 06 

£ 0 . 36.6.:3329 0 1 6>D— 0 1 0 • 53:9 37 05576:0— 06 

3 0. 3673534 759D- 01 0 .54355435680-06 

4 0.36782941 9 0D- 01 0 .54585945650-06 

5 0 . 368 08974 07D- 0 1 0 . 54694 1 £6530- 06 

6 0.36759004 020-01 0 .54940398840-06 

7 0. 36798860440-0 2 0 . 549225 18200-06 

8 0.36S0191726D-01 0.56134149770-06 

9 0 .2:686631 8660- 0 1 0 .55061 744 070-06 

10 0 . 368451 36330- 0 1 0 . 54952525890- 06 

11 0 . 36332233860- 01 0 . 549 1 3 1 56960- 06 

18 0.2:6934369930—01 0 .5346.25466:90— 06* 

13 0.36915434580- 01 0.53789922420-06 

14 0.36846533300-01 0 .54396828000- 06 

15 0.36727252040-01 0 .551 2047322D-06 

NO* 1 SAMPLING PERIOD* 3 .67273D-02SEC AMPLITUDE* 5.512050-07 NT=15 


Figure 2-6a 
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LST $ V S TE ft- NUMEP I C fiL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE > K2- 1*00 CUf— 0 1 j* K3- J.000D-01 
GAMMA* 1.3300011 07 M* 83 PH 1= 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 . 3673 7353 04 D- 0 1 0.551 £ 04 7333 O- 06 

1 0 .33 04351 071 D- 01 0.537 65 586 63H - 06 

3 0 .384 0576574 B- 01 0 . 52979547470-06 

3 0 . 33 53 4 684440- 0 l 0.53394057230-06 

4 0.38530 0 03 1 9 0 - 0 1 0 . 5363532 059 D - 06 

5 0 • 386 O 9783 920- 0 1 0 . 537346 1 5390- 06 

6 0 . 33809 076390- 0 1 0 .53359670330-06 

7 0 . 33633698630- 0 1 0 . 53737 1 33640-06 

8 0 . 3 3 6359 £3 9 0 0 - 0 1 0 . 5 4 0 77508670- 06 

9 0 . 33649 4 65 84 D - 0 1 0 . 53 7 66687 63 0 — 06 

1 0 0 . 386 544 33 970-01 0.5 36 9 97347 00 - 06 

11 0 .38658376630-01 0.53673330890-06 

NO= 1 SAMPLING PERIOD* 3 .366370- G3SEC AMPLITUDE* 5.364390-07 NT=11 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMFLITUDE CASE • K£= 1 . 0000- 01 . K3= 1.0000-01 
GAMMA* 1.380000 07 N= 84 PHI* 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 . 3 866 36643 10-01 0.5 36 4 385 33 1 D - 06 

1 0.40 033649760- 01 0.515 03497450-06 

£ 0.4 039 1001 340- 01 0.51 3 07 077470- 06 

3 0.4 05223392 0D- 01 0 . 533 1 6 09357D- 06 

4 0 . 4 05 8 49 03750- 0 1 0 . 5343 7174330-06 

5 0.4061 86 Q4450-01 0 . 53 54 0 338 1 9 0 - 06 

6 0.4071 0843310-01 0 .52393744390-06 

7 0.405 9 6 6 3 0 2 3 D - 0 1 U .53 0 999347 6 0 - 0 6 

8 0.4 06373 386 1 0 - 0 1 0 , 52643 039 6 1 0 - 06 

9 0.4 063 7335430- 01 0 . 53643 143750-06 

ND* 1 SAMPLING PERIOD* 4 . 063730- 03SEC AMPLITUDE* 5. 364140-07 NT* 9 


LST SYSTEM-NUMERICAL SOLUTION CP DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE * K2= 1.0000-01* K3= U000D-01 
GAMMA* 1 . 38 0 0 OD 07 N= 3 0 PH I * 0.0 

ITERATION SAMPLING PEP I OH AMPLITUDE 


0 

0 . 4 063 73141 8D- 0 1 

0 .58641439040-06 

1 

0 .4319 079343D- 01 

0.50147044720-06 

£ 

0.436 13495 ISO- 01 

0 .5035421 1610-06 

•Hi 

•J 

0. 43766 07503D- 01 

0 .507310192 30 - 06 

4 

0.4283 S 9 5 03 4 D - 0 1 

C . 5 0 9 3 6 3 92270- 0 6 

5 

0.428779 04 1 ID- 01 

0.51 032917430-06 

6 

0.429416 329 2D- 0 1 

0.5 09741 12270- 0'6 

7 

0.42957395760-01 

0 .5034 22766 8 0— 0 6 

8 

0 .429959951 00-01 

0.50613523770-06 

9 

0.42932332410-01 

0.51206561 170-06 

10 

0 . 4 2 8 9 9 7452 9 0 - 0 1 

0 . 5 1 3 J 56 4 094 0 - 0 6 

1 1 

0.4 29723 4 92 U 0-01 

0 . 5 C> 6 3 0 3 6 4 9 0 D - 0 6 

12 

0 .42 123535220-01 

0 .531 05 03269 D— 06 

13 

0.4259 04 67670-01 

0.5 18693641 2D- 06 

14 

0.4 2772539 180-01 

0 .514706569 ID- 06 

15 

0.42341 138280-01 

0.51374 1 4959D- 06 


ND* 1 SAMPLING PERIOD* 4 .334 12D-02SEC AMPLITUDE* 5.13741D-07 NT* 15 


Figure 2-6b. 
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LST SYSTEM-NIJMEPICflL solution of discrete describing function 
TWO AXIS-EQUAL AMPLITUDE CASE * K£* l.OOOD-Oli K3= 1. GOOD' 01 


GAMMA* 1.380001’ 07 N= 76 
ITERATION SAMPLING PERIOD 

0 0.48841 18G28D- 01 

1 0.44641 18167D -01 

8 0.45 134 0469 01'- 01 

• 3 0.4531456953D-QI 

4 0. 4540 097269 D- 01 

5 O'. 454469 024 0D- 01 

6 0.4550158575H-01 

7 0.4550362167D-01 


PH 1= 0.0 

AMPLITUDE 
0 % 5 1 374 1 4959D- 06 
0 • 484 1 942659B— 06 
0 . 4 3 5 0 8 0 0 6 2 3 B — 0 6 
0 ♦ 4 3 83325 027 1'- 06 
0 .4902121 173B-06 
0.481 056 1 5 02 D- 06 
0.4? 08 03 1 82 0 B- 0 6 
0 » 4? 06 1 0858 OB- 0 6 


ND* 1 SAMPLING PERIOD- 4 .55056D-02SEC AMPLITUDE* 4.9050SD-07 


NT* 7 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE * k2* t.OOOD-Oi* K3* 1 .00 0D- 01 

PH I = 0.0 

AMF'L I TUBE 
0 .4905 0 778 0 2D- 0 6 
0.46 02676 06 OB- 06 
0 .4611 052328D-06 
0 .463 9 223 4 4 4 D - 06 
0 .46542791 02D- 06 
0 . 4 66 1 0 08 946 D - 06 
0 . 466 0 1 26242D- 06 
0.4 6599596 17D- 06 

ND= 1 SAMPLING PERIOD* 4 .84594D- 02SEC AMPLITUDE* 4.65983D-07 NT* 7 


LST SYSTEM-NUMERICAL SOLUTION DF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE i K£= l . 00 OB- 01 > K3= 1 .000D-01 


GAMMA* 1.38 00 OB 07 N= 63 PHI= 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 



0 

0.434594 23 2 8D - 0 1 

0 .46593331940-06 


1 

0 .50771 233 1 9 D - 0 1 

0 . 43 0554 1 22 OD- 06 


c. 

0.5141 496277B- 0 1 

0 .430 03 16467B- 06 


3 

0 .51676971 85D- 0 1 

0.4 3 1961 27330-06 


4 

0 .5 1 3 06 593 9 5 D- 01 

0 .432362&233D-06 


cr 

0 .518751991 Oil- 01 

0 .43337333 OSD- 06 


6 

0.51 9£”'3£33SIt- u 1 

0.4 333 5 556 9 3D- 06 


7 

0.519S794877D-01 

0 . 4 8 8 3 3 8 1 2 o 5 D — U 6 

ND* 1 

SAMPLING PERIOD* 5.192 

300- 022 EC AMPLITUDE* 4 .333230-07 


GAMMA* 1.38000D 07 N= 72 
I TERH T I ON SAMPL I NG PER I OD 

0 0.455 0564878 D - 0 1 

1 0.4 74 8 4 92 58 2D- 01 

2 0 . 4 3 0 3 27 07 3 5 D - 0 1 

3 0. 482469 0167D- 01 

4 0.4335 15844 ID- 01 

5 0.4340710023D-0i 

6 0 .4345336 027D-0 1 

7 0.43459 138230- 01 


LST SYSTEM- 

NUMERICAL SOLUTION 

OF DISCRETE DESCRIBING FUNCTION 

TWO AXIS-EQUAL AMPLITUDE CASE 

? K£* 1 . OOOB-O 1 > K3= 1.000D-01 

GAMMA* 1.38 

0000 07 N= 64 

X 

»^-i 

ii 

Q 

4 

ITERATION 

SAMPLING PERIOD 

AMPLITUDE 

0 

0.51 92797 1 4 2D- 0 1 

0 . 4 3332 78 53 3 D - 0 6 

1 

0 . 54 6 7 72 4 4 98 D- 0 1 

0 . 39 106 7289 0 D - 0 6 

£ 

0.5544 07 16 OSD- 01 

0 . 3 8 8 9 3 8 0 3 7 8 D - 0 6 

3 

0.55767301 12D-01 

0 . 38965 0 1 9580- 06 

4 

0.5 59323721 80- 01 

0 . 3 39 9 3 4 3 0 4 6 D - 0 6 

5 

0 ■ 56- 0 1 *9 4 6 ij 4 ij D— 01 

0 .38996842170-06 

6 

0.56 074 4 17690-01 

0 . 3899899 1 26D- 06 

ND* 1 SAMPLING PERIOD* 5 .60744D-02SEC AMPLITUDE* 3.9001 0D-07 


Figure 2-6c 
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LST SYSTEM— NUMER ICRL SOLUTION □ F DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE ChSE > K2= 1. 0000-01? K3= 1. 0000-01 
GAMMA®* 1 -38 0 0 01* 07 N= 8 0 PHI= 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 .56074407240-0 1 0 . 38 0 0 1 043 0 1 B- 06 

1 0 . 5943855 1 84 B- 0 1 0 .3377828 184 D- 06 

£ 0,60343573990-01 0 .334 09255050-06 

3 0 . 607564 96 150-01 0 . 3332972694 B- 06 

4 0.6 09 6 So 0899 D- 0 1 0 . 3326 4 5 1 74 4 D- Ob 

5 0.61 OS 1 4 36 OSD- 0 1 0 . 332 1 7 74 34 70- 06 

6 0.61139 9 86 020-01 0 . 3 32 1931 0 1 3D- 06 

7 0.61 139934770- 0 1 0 . 32 S 9 0 3 4 3 5 0 D - 0 6 

8 0.61091761 63 D - 0 1 0 . 332 1 8569640- 0 6 

9 0.611361 629 40-01 0 - 332 4 4 252 050-06 

10 0.61136 064650-01 0 - 33773 1 3 7260 - 06 

11 0.61201 S6 225 D- 0 1 0 . 33 1 35842940- 06 

12 0-61203975190-01 0 .33 016222660-06 

13 0.6121 3374360- 0 1 0 . 3297677 1 550- 06 

14 0.612181 16190-01 0 -3295921224H-06 

ND= 1 SAMPLING PERIOD- 6 . 12248D-02SEC HMFL I TUDE= 3.293510-07 HTM4 


LST SYSTEM-NUMERICAL SOLUTION DF DISCRETE DESCRIBING FUNCTION 
TUD AXIS-EQUAL AMPLITUDE CASE ? K2= I . 00 00-01 > K3- 1.0000-01 
GAMMA = 1 ,380000 07 N= 56 PHI= 0.0 


ITERATION 

SAMPLING PER I CD 

AMPLITUDE 

0 

0 -6122479469 D - U 1 

Ci .329 3 5 0 9 S 7 8 D - 0 6 

1 

0.654 33 01 76 ID- 01 

0 . 26395775 080- 06 

2 

0.66451732420-01 

0 . 2 6 0 8 2 5 3 3 8 2 D - 0 6* 


0 . 6 6 957 9 794 0 1* - 0 1 

0 . 258 94942440- 06 

4 

0.67 2 2 6 0 4 4 0 6 Is - 0 i 

0 .25743746 050 - 0 6 

5 

0 . 673 7 0 1 9 84 8 0 — 0 1 

0.25643975950-06 

6 

0.674 3 8 9 6 9 0 3 0 - 0 1 

0 . 8 5 6 8 8 7 9 3 0 6- D - 0 6 

— * 
r 

0.6743 9 328 470-01 

0,25571570620- 06 

3 

0.67440923880-01 

0.25590943330-06 


ND= 1 SAMPLING PEP I OB= 6 . 744 1 7D-G8SEC AMPLITUDE®* £.558430-07 NT= 8 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE? K2= 1 - 0 0 OB- 01 ? K3= 1.0000-01 
GflMMA= 1 .380000 07 N= 52 PHI= 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 . 6 7 4 4 1 669 7?D - 0 1 0-25584 3 0 5 79 D - 0 6 

1 U . 738 85 077270-01 0 .15921 66 0 050- 0 6 

2 0 - 75 1 88299830- 0 1 0.1 18 03763520- 06 

3 0.75441970060-01 0 . 14 148552240-06 

4 0 . 75433931870 - 0 1 0 . 1 53988 11110- 06 

5 0 .7540 0 8 72 84 D - 0 1 0. 15779657470-06 

6 0.7534 0 5 23 23 It- 01 0.16 0 0 97146 3 B- 0 6 

7 0.7534 6 75354 0 -01 0.1 589 1 1 338 40- 06 

8 0 - 765636 1 7 1 30-0 1 0.1 0522686390- 06 

9 0.760 9634 6760-01 0.129753 788 1 0 - 06 

1 0 0 .75726425640-01 0 .14431236980-0 6 

11 0 . 755 1 35924 1 D- 0 1 0.1 53 1 4975930- 06 

12 0 . 754 1 6 1 3 1620-01 0.1 5737 053870- 06 

13 0.75 353 5 3854 0 -01 0.1 59 7 0 03 993 0-06 

14 0.7537343 09 8 D -01 0.1 5*79 0774 9 3 D - 06 

15 0.75330 86 1 350 - 0 1 0.1 6 0 0 0554 09 D- 0 6 

Nrr= 1 SAMPLING PERIOK= ?.53309D-0£$EC fiMPLITUIiE= 1 .60006D-07 NT=15 


Figure 2-6d. 
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LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 


TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .330000 07 N= 100 
ITERATION SAMPLING PERIOD 

0 0 .15300000000-01 

1 0-15131 £9565Ii- 0 1 

2 0 - 15 123 6333 9D- 0 1 

3 0.1 5 09 5334 98 D- 01 

4 0.150743 09980- 0 1 

5 0. 15075774 19D- 01 

6 0.150S043353D-01 

7 0.1505980197D-01 

8 0.1505978776D-01 


, K2=K3=0. 1 
PH I = O.U 

AMPLITUDE 
0.500000000 OB- 05 
0.5123194566D-05 
0.51S6988252D-05 
0 . 52 1 38 088 09D- 05 
0 .524 0422 02 4 D- 05 
0.5231 844 04SP- 05 
0. 5225 764 83 0D- 05 
0.5&53128272D-05 
0 .5253562257D-05 


NB= 1 SAMPLING PERIOD- 1 .5059SP-02SEC AMPL I TUDE= 5.253S7D-06 N=* 8 


LST SYSTEM-NUMERICAL SOLUTION 
TDD AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .330 GOD 07 N= 98 
ITERATION SAMPLING' PERIOD 

0 0. 1505977 177D-01 

1 0.1524686473D-01 

2 0. 1 545879 048B- 01 

3 0.1 5568563 06D- 01 

4 0 . 156269746 ID- 01 

5 0.1565694349E-01 

6 0.1 566 1 06863D- 0 1 


□F DISCRETE DESCRIBING FUNCTION 
■ K2= K3=0 . 1 
PHI = 0.0 

AMPLITUDE 

0 .5253S65498D- 05 . 

0.542421 094QD- 05 
0.535 063779 ID- 05 
0.53 0654 7 0 99D - 05 
0.523255 009 ID- 05 
0.527 G912439D-05 
0.5263612074D-05 


ND= 1 SAMPLING PERIOD* 1 -56629D-02SEC AMPLITUDE* 5.26612D-06 M* 6 


LST SYSTEM-NUMERICAL SOLUTION 
TUG AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1.38000B 07 N* 92 
ITERATION SAMPLING PERIOD 

0 0.1 5 66289358 D- 0 1 

1 0.1 5369283 G7D- 01 

2 0.161 026951 OH- 01 

3 0.1 622686675B - 0 1 

4 0.1689 157791 D- 01 

5 0.1 63250 0326 B- 0 1 

6 0.1 633 08464 ID- 01 


OF DISCRETE DESCRIBING FUNCTION 
* K£= K3-0. 1 
PHI- 0.0 

AMPLITUDE 
0 .52661 1 679 3D— 05 
0. 5447574799H-05 
0. 53679933 14D- 05 
0.532 04 0788 ID- 05 
0. 5294671 135D-05 
0. 528 19 0568 OH- 05 
0. 527968 1922D- 05 


NP* 1 SAMPLING PERIOD* 1 .63327D-02SEC AMPLITUDE* 5.8766 4L-06 N= 6 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .330 ODD 07 N= 88 
ITERATION SAMPLING PERIOD 

0 0 . 1 6332726 06D- 0 1 

1 0.1 655487224 D-Oi 

2 0.1681 022992D- 01 

3 0,1 694 596 0 3 4 D - 0 1 

4 0.1701 655283 P-01 

5 0.170531004 OB- 0 1 

6 0 . 170605961 6D-01 


OF DISCRETE DESCRIBING FUNCTION 
> K2= K3=0. 1 
PHI* 0.0 

AMPLITUDE 
0.5276636633D-05 
0. 54691 592440-05 
0.538474393 ID- 05 
0.533445651 ID- 05 
0 . 53 073 7484 ID- 05 
0 .52937937270- 05 
0 .529 1479667D- 05 


ND= 1 SAMPLING PERIOD* 1 . 70625D- 02SEC AMPLITUDE* 5. 288 04 D- 06 N* 6 


Figure 2-7a. 
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LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE* K2=K3=0.1 
GAMMA* 1 .33 GOOD 07 N= 84 PHI= 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0.17062493980-01 0 .52880353050-05 

1 0.173021691011-01 0.54929651860-05 

£ 0.1 75826422 0D- 01 0 .54 032 1827 ID- 05 

3 0.17731596290-01 0 .53499600750-05 

4 0.1 780891 45 ID- 01 0 .53213975851-05 

5 0.1 73489897 0D- 0 1 0 .53069635220-05 

6 0.17858572150-01 0 .53043933190-05 

ND = 1 SAMPLING PERIOB= 1 . 7S605D-02SEC AMPLITUDE* 5.30055D-06 N= 6 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE t K2=K3=0.1 
6AMMA= 1.380000 07 N= 80 PHI* 0.0 

ITEPATION SAMPLING PERIOD AMPLITUDE 

0 0.17860458440-01 0.53003534160-05 

1 0. 18119822270-01 0 .55 1935 16 1 3D- 05 

£ 0. 1842922889D-01 0 .54237 09952D-05 

3 0 . 1 85934 1236D- 01 0.53671654640-05 

4 0. 1 8678473920- 01 0 .53369647790-05 

5 0.1 372256983D- 0 1 0 .53216139150-05 

6 0 . 1873473955:0-01 0.5318625613D-05 

HD* 1 SAMPLING PERIOD* 1 .S7366D-02SEC AMPLITUDE* 5.31440B-06 N* 6 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE? K£=K3=0.1 
GAMMA* 1.380000 07 N= 76 PHI* O.n 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 . 1373659666D-01 0 .531 44 04685D-05 

1 0.190182081411-01 0 .55487497270-05 

£ 0.19361 16657D-01 0 .54465681 160-05 

3 0 . 19543 00373D- 01 0.53863750320-05 

4 0 . 19637038700-01 0.53543550620-05 

5 0.1 96857552 1 0- 0 1 0 . 5338 0 1 62280- 05 

6 0.1 97 01 065830-01 0.53344121320-05 

NO* 1 SAMPLING PERIOD* 1 .97029D-02SEC AMPLITUDE* 5.32987D-06 N* 6 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE > K2=K3=0.I 
GAMMA* 1.33000D 07 N= 72 PHI* 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0.197 02898 1 7D - 0 1 0 .532987432 1 0-05 

1 0 .200 09833 75 D- 0 1 0.55816925820-05 

£ 0.20391 980820- 0 1 0 .54722286910-05 

3 0.2059442799P-01 0.54079689800-05 

4 0 .20698935900- 01 0 . 53739204 000-05 

5 0.2 0753 0 08 190-01 0 . 53565 078490- 05 

6 0.20772052880-01 0.53520941850-05 

ND= 1 SAMPLING PERIOD* 2 . 07739D-02SEC AMPLITUDE* 5.34734D-06 N= 6 


Figure 2-7b 
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LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE. K2=K3=0.1 
GAMMA= 1 .3800011 07 N= 68 
ITERATION SAMPLING PERIOD 

0 0 . £ 077383646 D- 0 1 

1 0 .21 1 09907160-01 

2 0.21538191750-01 

3 0.21764388040-01 

4 0.218317144 00- 0 1 

5 0 .2 1342 05 08 ID- 01 

6 0.21965454620-01 

7 0.213673871 10-01 

ND= 1 SAMPLING PERIOD* £ . 1 36980- 02SEC AMPLITUDE* 5. 366330-06 N= 7 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE. K£=K3-0.1 
GAMMA* 1.380000 07 N= 64 PHI* 0.0 


ITERATION 

SAMPLING PERIOD 

AMPLITUDE 

0 

0. £1 96979444 D- 01 

0.5366325 1661'- 05 

1 

0 . ££33797565 D - 0 1 

0.566 086 4 35 7 D- 05 

o 

0 .££820653170-01 

0 .55341 3962 3D- 05 

3 

0 .£30759 0 01 5D- 01 

0.546 03238 16D- 05 

4 

0 .2320721 9 GOD- 01 

0.542I483097D-05 

5 

0 .£327489 029D-01 

0.540161 £733 D- 05 

6 

0. £33 03257930- 01 

0 . 5395 08 7 322 D- 05 

7 

0 .2330541216D- 01 

0.539 0 33 1951 B- 05 


ND= 1 SAMPLING PERIOD* £ .33080B-02SEC AMPLITUDE* 5. 339320-06 N= 7 


PH I = 0.0 

AMPLITUDE 
0.53473332900-05 
0 .561S843722D-05 
0 .55 012438970-05 
0.54324326790-05 
0.53961 123750-05 
0.53775230770-05 
0.53721205190-05 
0.53672837760-05 


LST S VS TEM-NUMERI C A L SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1.380000 07 H= 60 
ITERATION SAMPLING PERIOD 

0 i} • £ 33 08 0 03761!- 0 1 

1 0 . £ 3715464671*- 01 

£ 0 .2 4263 36 869 B- 01 

3 0. £455357339 D- 01 

4 0.2470238392D-0! 

I 5 0.347733755 ID- 01 

6 0 • £4 8 1 3 0 0 1 77 B - 0 1 

7 0 .3431 56 033 IB- 01 


OF DISCRETE DESCRIBING FUNCTION 
, K£=K3=0.1 
PH 1= 0.0 

AMPLITUDE 
0 .538931 7 0 £ £ D ~ u 5 
0 . 57 09c 7 0 1 £5B- 05 
0.55722739 ODD- 05 
0.549£6E4723D-05 
0.545 0394 726D- 05 
0. 5429575 152D- 05 
Q. 5421 7353 19D- 05 
0.541 72395551'- 05 


ND= i SAMPLING PERIOD- 2.4S1S5P- 02SEC AMPLITUDE* 5.41611D-0<$ N= 7 


LST SYSTEM-NUMERICAL SOLUTION 
TWD AXIS-EQUAL AMPLITUDE CASE 
GANNA = 1 .330000 07 H= 56 
ITERATION SAMPLING PERIOD 

0 0. £4318457751'- 01 

1 0 .£5272££59 0D-0I 

£ 0.25900371 16D- 01 

3 0.2623180446D-01 

4 0 .26401 77621 D- 01 

5 0 . £64 88 9 0 4 7 £ E - 0 l 

6 0. 265296 4 524 D- 01 

7 • 0 ,£65330774 ID- 01 


OF DISCRETE DESCRIBING FUNCTION 
, K£=K3“0* 1 
PH I - O.U 

AMPLITUDE 
0.54 161 09 1 9 OB- 05 
0 .5765341 9£ 3 D - 05 
0.561661 68 1 OB- 05 
0 . 553 0395 04 3D- 05 
0 .548538 14? OB- 05 
0 . 546 £ 4 3 4 53 1 B- 0 5 
0.54532S6148D-05 
' 0 ,54492 1 1 264D-05 


ND= 1 SAMPLING PERIOD* £ . 65363D- 02SEC ANPLITUDE= 5.44773D-06 N= 7 


Figure 2-7c 
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LSI SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE* K2=K3=0,1 
GAMMA® 1.38 00 OB 07 N= 5£ PHI® 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 * £653626294 B- 0 1 0 . 5447726 057D- OS 

1 0 . £7 045 1 63 65B- 0 1 0 . 583 1 0664 7 4 D- 05 

£ 0 - £777 1 26829 D- 01 0 . 56639293 05 B- 05 

3 0 . 23 1 53301 32 D- 01 0.557520 1 503D- 05 

4 0 .£83496564711-01 0 .552644 016SD-05 

5 0.2344 975378 Ii-01 0 .55 01 65337 ID- 05 

6 0 . £3498 171 5£B- 0 1 0.5491 072646D-05 

7 0 - 285 0 3 0 95 6 7 D - 0 1 0 . 5 43 7 5 268 44 D - 05 

NB= 1 SAMPLING PERIOD* 2 .85067D-02SEC AMPLITUDE- 5.48557D-06 N= 7 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS -EQUAL AMPLITUDE CASE ? K£= K3-0.1 
GAMMA- 1 , 3 3 0 0 0 D 07 N= 48 PH I = 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 . £85 066 533 2D- 0 1 0 .548556876SD-05 

1 0 .29 03 1 82 1 9 OB- 0 1 0 .5909 1 954 09D- 05 

£ 0.29 929 9 76 6 3 D - 0 1 0 .5731615644 D- 05 

3 Ci .3 03753691611-01 G . 5629273377H- 05 

4 0.30 6 0375 8£3 D - 0 1 0 .55762247640- 05 

5 0 - 3 0 7 1 9875 17D- 0 1 0 . 554934 89771' - 0 5 

6 0 .3 077733 074D-0 1 0 .553736322 OP- 05 

7 0 . 3 0 784896 5SB- 0 1 0 .553425253 IB- C'5 

NB= 1 SAMPLING PERIOD- 3 . 07S88D- 02SEC AMPLITUDE- 5.53I63D-06 N® 7 


LST SYSTEM-NUMERICAL SOLUTION DF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE * K£=K3«0.1 
GAMMh= 1 .3S00GD 07 N- 44 PHI- 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0-3 078877 1 4 1 D- 0 1 0 . 553163 1 333D- 05 

1 0.31444 39 05 3 D - 0 1 0.6 0 03648 0 1 0 D - 05 

2 0.3244 70986 1 D- 0 1 0 . 58 08 1 558 040- 05 

3 0 .3297263 07 0D- 01 0 . 56959 141730- 05 

4 0 . 3324 066 9 35 D- 0 1 0 . 5637937 0520- 05 

5 0 . 3 3 3768 0368 D - 0 1 0 .56086816970-05 

6 0 .33445 1 293 ID-01 0.5595031 0 £ 00 - 05 

7 0 .3345727 029B- 0 1 0.55923618380-05 

ND= 1 SAMPLING PERIOD* 3 . 346 1 3D- 02SEC AMPLITUDE- 5.588890-06 N® 7 


Figure 2 -7 d 
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LST SYSTEM- MUM EP I CRL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA- 1-33 GOOD 07 N* £ 0 
ITERATION SAMPLING PERIOD 

0 0.6000000 0 0 OB- 0 1 

1 0.648 08799 34 D - 0 1 

£ 0 . 6673 026 0 1 9 B - 0 1 

3 0.6769888 12 OH- 01 

4 0.631 3 04334 1 1*- 0 1 

5 0 .633425 04290-01 

6 0 .6843 0 89 59 4 D - 0 1 

7 0 .63446455430-01 

3 0 .6846078 1 730-01 


□F DISCRETE DESCRIBING PUNCH ON 
K2=fC3=0. 1 
PHI* 0 . U 

AMPLITUDE 
0.100000 0 0 0 OB- 04 
0 . 761 7689854 D- 05 
0 .71 72493 189D- 05 
0 . 6996828 138D- 05 
0. 69 124331 5 0D- 05 
0 . 668853 1 1 39B- 05 
0.68 63 1 9 06 1 2 B - 05 
0.63581 689 06D- 05 
0 .685474493 OB- 05 


ND* 1 SAMPLING PERIOD* 6 . 34679B- 02SEC AMPLITUDE* 6.S5169B-06 


M= 3 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .33 00 OB 07 N= 84 
I TEPAT I ON 8 A NFL I NG PERI DU 

0 0 . 6 8467379 0 4 D - 0 1 

1 0.6 78 9 6 889 08 B - 0 1 

8 0 . 63 373 6 1 6 58D - 0 1 

3 0 . 533 1 598 1 93B- 0 1 

4 0.57851 4 7883 D - 0 1 

5 0.57 36 9 72 579 D - 0 1 

6 0.579838 3375 B - 0 1 

7 0 . 58 £ OS 8998 4 D - 0 1 

8 0.583 565475 IB- 01 

9 0 .564 8629 65 8 B - 0 1 

10 0 . 5 8483 7 0888 D- 0 1 


DF DISCRETE DESCRIBING FUNCTION 
K2=K3*0.1 

FHI= o:o 

AMPLITUDE 
0 . 6 6: 5 1692 6 6. 3 D - 0 5 
0 . 3 84 9 7 9 7 3 5 3 P - 05 
0 . 4653 6444 34 D- 0 5 
0.51 1355551 70- 05 
0 .61365 08163D- 05 
0 .643 7S 8545 3 B- 0 5 
0.6435582S88B-05 
0 .64 06 344573 D- 0 5 
0.63 9 0 8 0 8 0 9 6 D - 0 5 
0 . 6 3 8 5 6 8 1 0 9 5D- 05 
0 . 6 3 8 1 2 38 2 6 5 D - 0 5 


NB= 1 SAMPLING PERIOD* 5 . S4322D-02SEC AMPLITUDE* 6 .331 06D-06 N=10 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .330 00H 07 H= 88 
ITERATION SAMPLING PERIOD 

0 0.5 84 381 55 89B- 01 

1 0. 532 0396 04 7D- 01 

a 0 . 55 088 0467 7B- 0 1 

3 0.531 82428 04D- 01 

4 0.581 08 1 9 143D- 0 1 

5 0.51 528529 OSD- 0 1 

6 0.5 1243443 OSD- 01 

7 G .51 1 1 741517B-01 

S 0.51 06 04657811-01 

. 9 0. 5 095 8 08-379 D- 01 

10 0. 5 095684 164 D- 01 

1 1 0.50951204430-01 

12 0 . 5 096 9 373 8 5 B - 0 1 

13 0.5 09694 0058D- 01 


GF DISCRETE DESCRIBING FUNCTION 
K2-K3=0. 1 
PHI* 0.0 

AMPLITUDE 
0 .6331 06312 IB- 05 
0 . 4 079 36 1 1 3 4 D- 05 
0 .477033 873 9 D- 0 5 
0 . 529893735 1 D- 05 
0.565 04 67287 D- 05 
0 ; 58 5891 a 9 6 9 D - 0 5 
0 . 59 7 26 04748 D - 05 
m , 6 03 0 0 58 9 09 D - 05 
0 . 6 0568 0 1 388 B- 05 
0.6100 9 3 71 0 1 D - 05 
0 . 6 0388 73 9 42 D- 05 
0 . 6 0 74 0 84 692B-05 
0 .60805 01 70 OB- 05 
0 . 6G77955 158B-G5 


MB* 1 SAMPLING PERIOD* 5 . 09695B- 02SEC AMPLITUDE* 6.07377B-06 N*13 


Figure 2-7 Q. 
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LST SYSTEM- NUMERICAL SOLUTION Of DISCRETE BESCR I B INb FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE K2-K3-0.1 


Gfthf1R= 1 .38 

00 CD 07 N= 38 

PHI® 0.0 

ITERATION 

SAMPLING PERIOD 

AMPLITUDE 

0 

0 . 5 0969458 0 1 D- 0 1 

0,6 07377 0S68D-05 

1 

0.51 03497983D-01 

6.41 1916414 011-05 

3 

0 .434 3473884 D- 0 1 

0 .4757019 016D-05 

3 

0 . 4 639633 6 0 3 D - 0 1 

0 . 5 2 363 3 4 63 9 D- 0 5 

4 

0 . 46 04 573883 D- 0 1 

0 . 5527345 1 S6B- 05 

5 

0 -456 03 081 OSD- 01 

0.570221 0'?4 70-05 

6 

0.45333561 13D-01 

0 .57951 55652D- 05 

7 

0.45£79£1£9?B-01 

0 .584 155 134 2D- 05 

S 

0 . 4 53 1 94 466 5 D - 0 1 

0. 53662 12535P- 05 

9 

0.454730 08 7 9 D— 0 1 

0. 57279 16338D- 05 

in 

0 .45863 156 08 D- 01 

0.5314044 1 94 D- 05 

n 

0. 4537433954 D-nt 

0 .58541 4 7556 D- 05 

18 

0.4534631 043H-01 

0-5381 19S353D-05 

13 

0.45644 13 164 D- 01 

0 - 56 037 0937 0 p - 0 5 

14 

0.4633 4 65 36 3D- 01 

0 .5387431 6 1 6 D- 0 5 

15 

0.4573 083770D-01 

0 .587 044391 ID- 05 

16 

0 .4545788389H-01 

0 - 537 3738 757 D- 05 

17 

0 .4530483 077D- 01 

’ 'I' . 5 3 S 3 6 ' 9 9 8 7 9 D - 0 5 

18 

0 .4519358759D-01 

0 . 59 0839 0757D- 05 

19 

0 .4513646 046 D- 0 1 

0 . 5385845 069 D- 05 

3 0 

0 .451 375651 3D- 01 

n . 5 878 83 0 0 96 1*- 05 


HD® 1 SAMPLING PERIOD® 4-51 376B- 08SEC AMPLITUDE® 5.37228B-G6 N=20 


LST SYSTEM-NUMERICAL SOLUTION DF DISCRETE DESCRIBING FUNCTION 
TWO AXIS -EQUAL AMPLITUDE CASE K2-K3=0.1 
GAMMA® 1.38 GOOD 07 M= 38 PHI= 0.0 

ITERATION SAMPLING PER I CD AMPLITUDE 

0 0-451 37585 1 3D- 0 1 0 . 587833 0 096P-05 

1 0-45434413 05 D - 0 1 0 . 4 079539734D- 05 

3 0.431 9483 1 35D- 0 1 0.470431 3647 D- 05 

3 0-41 3753533 7 D - 0 1 0 . 5 1 53 7 5 35 4 3 D- 0 5 

4 0.41161 39664 D- 0 1 0 . 5433353 03 00- 05 

5 0.40794 009640-01 0 - 55934333 76D- 05 

6 0.406151 63 88 D - 0 1 0 . 56 78 091 43 OP- 05 

7 0.4 05453 0647B- U 1 0 - 57 1 9556 0 1 3D- 05 

8 0.4 05 1715401 D- 0 1 0 . 57334 37457D- 05 

9 0.4 05808451 4 P- 01 0 . 57334671 74B- 05 

NP- 1 SAMPLING PERIOD® 4 . 05834D-03SEC AMPLITUDE® 5.78936P-06 N® 9 


Figure 2-7f 
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LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS- EQUAL AMPLITUDE CASE R2=K3=0.1 


GAMMA* 1.38 

00 OP 07 N* 4 0 

PHI* 0.0 

ITERATION 

SAMPLING PERIOD 

AMPL I TUBE 

0 

0.4 05884 1794 P-01 

0. 5789359 1 39 D- 05 

1 

0.41 05375766D-01 

0 . 3 9 9 3 9 4 5 09 9 B - 0 5 

2 

0 . 3886 4 0 7364 P - 0 1 

0.46425147350-05 

3 

0.3754 96435 OD- 01 

0 . 5 096797 069D- 05 

4 

0.379 4 399 22 OB- 01 

0.5485304 1 03 D- 05 

5 

0.37266 0291 ID- 01 

n .559624 65 070-05 

6 

0 . 3 69664 72 89 D- 0 1 

0 .56296388760-05 

7 

0 . 368 1 4 66 5 79 D- 0 1 

0 .56466 01267P-05 

8 

0 .3671 742546P-01 

0 .5669053 6 0 7 B - 0 5 

9 

0.36718362550- 0 1 

0 . 563776388 1 B- 05 

10 

0.3674531 1250-01 

0 • 56 1 36 6 49 39 P— 05 

11 

0 .36712907790-01 

0 . 56 3588 1 029 P - 0 5 

12 

0.36740227180-01 

0 .561 6 05831 SB- 05 

13 

0.367094069 OD- 01 

0 .56366344 69P — 0 5 

14 

0.36724 104 1 3D- 01 

0 .5624374 99 1 P - 0 5 

15 

0 . 3669 1 7262 IP- 0 1 

0 . 5645 1 32773P - 05 

16 

0 . 3669257664 p- 0 1 

0 . 5 6 4 0 1 3 0 3 4 3 P - 0 5 


NB= 1 SAMPLING PEEIOI*= 3 .£6*40D-0£SEC AMPLITUDE* 5 .G3392D-06 N=16 


LST SYS TEN-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA- 1 .380 0 OD 07 N= 44 
ITERATION SAMPLING PERIOD 

0 0.36694 04 01 9D- 01 

1 0.374 013853 011-01 

£ 0.3491 01 0532D- 01 

3 0 . 348 1 1 3 6 75 3D- 0 1 

4 0 .341 1 19590 OD- 01 

5 0 . 33 794 £ 56 9 7D - 0 1 

0 . 3o8o8tz93*i> ID - U 1 

7 0.33543674 7 OD- 01 

S’ 0 • 33546 78748P— 0 1 

9 0.3 36937893 1I<- 01 

10 0. 3358 04236 ID- 01 

1 1 0. 33541 434 35D- 01 

1 2 0 . 3348 766 7 36 P - 0 1 

13 0 . 334 3 753 55 7D- 0 1 

14 0-33527 19 03 ID- 01 

15 0 . 3 3 5 2 33 7 6 3 6 D - 0 1 

16 0.33521 19786D-01 

17 0.3352224 71 3D- 01 

13 0.3352487997D-01 


OF DISCRETE DESCRIBING FUNCTION 
K2=ft3=0. 1 
PHI' 0-0 

AMPLITUDE 
0 . 56339 1 9976 P- 05 
0 . 3 9 0 8 91553 3 B - 0 5 
0.45389346 3 4 B - 0 5 
0 . 53367 0382 ID- 05 
0 . 54 3 34986 04 B- 0 5 
0 . 55 4188 0 3 0 0 D - 0 5 
0 . 556 82 1 1 972D- 05 
0,559348 17 04 P-05 
0 . 556 0 5 335 74 D- 05 
0 .54571 1 02 1 4D-05 
0 .552477554 ID' 05 
0.555586 3 0 0 9 D - 0 5 
0.5539 3 98 1 8 OD- 0 5 
0 . 5 678 6 68 056 B- 0 5 
0 ,55263433490-05 
0.55 96966 8 64 D - 0 5 
fi .5575 070 0990-05 
0.55660 197 32D- 05 
0 . 55 63343 3630- 05 


HP- 1 SAMPLING PERIOD* 3 -35304D-02SEC AMPLITUDE* 5.56041D-06 N=13 


Figure 2-7 g 
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LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1.3SG00D 07 M- £0 
ITERATION SAMPLING PERIOD 

0 0.6000000000B-OI 

1 0-643 0379984 D- 0 1 

a 0.6678 026 0191'- 01 

3 0.67692321 SOD- 01 

4 0 .6813 04334 ID- 01 

5 0.6S34250429B-G1 

6 0.6843 083 5 94 D- 01 

7 0 . 6 34 4 645 5 4 3 P- 0 1 

8 0 -6846 072 1 73P- 0 1 


□F DISCRETE DESCRIBING FUNCTION 
K2=K3=0. 1 
PHI* 0.0 

AMPLITUDE 
0 . 1 0000 0 0 0 0 QD- 04 
0.7617629854B-0S 
0.71 72493 1 89D- 05 
0. 6996322 183P- 05 
0 . 69 1 8483 1 5 OD - 0 5 
0 .6888531 139 D- 05 
0 . 6868 1 9 06 1 8D— 05 
0.6 8 58 l 629 OSD- 05 
0 . 685474493 OD- 05 


ND* 1 SAMPLING PERIOD* 6 .84679D- 02SEC RMPLITUDE* 6.85169D-06 N = 8 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL RMPLITUDE CRSE 
GAMMA* 1.38 GOOD 07 N* IS 
ITERATION SAMPLING PERIOD 

0 0.6846 7879 04 D- 0 1 

1 0.7 0432567 07D- 01 

2 0.7274274026H-0 1 

3 0.7391 2561 51H-01 

4 0.7449677 6 09 D - 0 1 

5 O . 747896721 6 D— 0 1 

6 0.7494 03586 1 D - 0 1 

7 0.7500155196 D- 0 1 


OF DISCRETE DESCRIBING FUNCTION 
K2~K3-0. 1 

PHI* 0.0 

RMPLITUDE 
0 . 685 1 69266 3 D- 05 
0 . 74 94 237996 D- 05 
0. 73535732 7SD- 05 
0.727 08 074 4SD- 05 
0 . 7224 62 0882D — 05 
0 . 72 0121 9539D- 05 
0.71 3 926 9 55 6 D- 05 
0.7137 93 6 3 5 5 D- 0 5 


NB= I SAMPLING PERIOD* 7 .50031D-02SEC AMPLITUDE* 7.13146D-06 N= 7 


LST SYSTEM-NUMERICAL SOLUTION 
TWO MX IS- EC UAL AMPLITUDE CASE 
GAMMA* 1 . 33 0 0 OD 07 N= 16 
I TERR T I ON SRMPL’I NG PER 1 0 D 

0 0. 75 0 031 3376 D- 01 

1 0. 77572 07225 D- 01 

2 0 . 3 03 6526528 D - 0 1 

3 0.61771 0168411-01 

4 0 . 8247 1 7 578 5 D - 0 1 

5 0 .82822 03857 D- 0 1 

6 0 . 83 0 0 1 8 1 793 D- 0 1 

7 0.8309 1001450-01 

8 0 . 83 0 9249 7 0 7 D - 0 1 


□F DISCRETE DESCRIBING FUNCTION 
K2*K3*0. 1 
PHI* 0.0 

AMPLITUDE 
0.7181 455222D- 05 
0.79041 02 56 5 D- 05 
0 . 778 355 6 632 D- 05 
0 .76993299 1 2D- 05 
0 . 7654 03 1 462D- 05 
0 .7631 0 0342 ID- 05 
0.761 9 06 2 536 D- 05 
0.76 1567974 ID- 05 
0.761 06293 02D- 05 


ND= 1 SAMPLING PERIOD* 3 .30950D- 02SEC AMPLITUDE* 7.61094D-06 N= 3 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GhMMh* • 1 .3300 OD 07 N* 14 
ITERATION SAMPLING PERIOD 

0 0.83 0 9 5 0 1 0 7 3 D - 0 i 

1 0 • o6< 7 1 y 3 9 U 2 9 D — 0 1 

2 0 . 9 0 1 764 7394 D - 0 1 

3 0.91 9 0 6 913 S 3 D - 0 1 

4 0.927 688 133 ID -01 

5 0 . 93 1 336 1 2 04 D- 0 1 

6 0.934 175 0696D-01 

? 0 . 93 5 3 4 9 69 42 D - 0 1 

8 0. 9353941 6 07D- 01 


OF DISCRETE DESCRIBING FUNCTION 
K2*K3*0.1 

PHI* 0.0 

AMPLITUDE 
0 . 76 1 0937 6 09D— 05 
0 .8403642551 D- 05 
0 .3313391 1 06 D- 05 
0 . 32463 074 9SD- 05 
0 .82 055 1 0559D- 05 
0 . 3 1 8 4 55 8 6? 9 B- 05 
0. 3 1 73538624 D- 05 
0.31 68 7 2384 9D - 0 5 
0 .8 1 67792323D- 05 


ND= 1 SAMPLING PERIOD* 9 . 35426D- 02SEC AMPLITUDE* 3 . 1 663 9 D- 06 N= 8 


Figure 2-7h 
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LSI SYSTEM-NUMERICAL SOLUTION 
TWO AX I S- EQUAL AMPLITUDE CASE 
GAMMA* 1 * 33 0 GOD 07 N* . 12 
ITERATION SAMPLING PERIOD 

0 0. 9354264 7740- 01 

1 0 . 991611 * 22211-01 

2 0.1 035439723 D 00 

3 0.1 05 732536 3D 00 

4 0.106S237479B 00 

5 0.1 073668 119D 00 

6 0.1 0764 09935D 00 

7 0.1 07 7* 3993 7 0 00 

8 0 . 1 0781605240 0 0 


OF DISCRETE DESCRIBING FUNCTION 
K2=K3=0. 1 
PH I = 0.0 

AMPLITUDE 
0.816688*57990-05 
0 . 899 0 1 56255D - U5 
0 .8984 1 62 086D— 05 
0.89441 19276D- 05 
0 . 89167459650-05 
0.89 02 207 03 4 D- 05 
0.38945136160-05 
0 .889 00917720-05 
0.8894539016D-05 


ND= 1 SRMPLINtS PER 100* 1 . 073 180- 01 SEC AMPLITUDE= 3.889440-06 N* 8 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1.38000D 07 N* 10 
ITERATION SAMPLING PERIOD 

0 0.1 07817811 3D 00 

1 0.11771 00597D 00 

2 0.1 £323362 1 90 00 

3 0.18605038430 00 

4 0. 127463968 OD 00 

5 0 . 12© 167 038 OD 0 0 

6 0. 12851 8898 OD 0 0 

7 0 . 1287 044 745D 00 

8 0 . 1 288 1 9.6*3330 00 


OF DISCRETE DESCRIBING FUNCTION 

K2-K3-0.1 

PH I = 0.0 

AMPLITUDE 
0 . 38 8 9444 8 3 ? D - 0 5 
0.962 14 15477D-05 
0 . 93 0 1365517 D- 0 5 
0 . 984 06989220- 05 
0 .98492 1 1 9260- u5 
0 - 9 85224 5354 D- 05 
0.985361 16990-05 
0 . 9 8 5 3 3 6 6 9 0 2 D - 0 5 
0 - 9855 £355 5 3D- 05 


N0= 1 SAMPLING PERIOD- 1 . £882 00-01 SEC AMPLITUDE- 9. 854370-06 N* 3 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .330000 07 N* 8 
ITERATION SAMPLING PERIOD 

0 0 . i8o8195375P 00 

1 0,14973713370 00 

2 O.15532OS5320 00 

3 0 . 1 5937384570 00 

4 0. 16120071880 00 

5 0.162113562 60 00 

6 0. 16856771900 00 

7 0 . 1628001760D 00 

8 0.16894753050 00 


DF DISCRETE DESCRIBING FUNCTION 
K2-K3=0. 1 
PHI* 0.0 

AMPLITUDE 

0 .98543657740-05 
0. DO 1988965 00- 04 
0 . 1 0804674160-04 
0,11028858850-04 
0.1 11 86355 7 8D- 04 
0.11174860510-04 
0.11 19856757D-04 
0.1181 0886330-04 
0.1 1213291 150-04 


ND= 1 SAMPLING PERIOD* 1 .629480- 01 SEC AMPLITUDE* 1.123230-05 N= 8 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .380000 07 M* 6 
ITERATION SAMPLING PERIOD 

0 0 . 1629 4 83 72 2D 00 

1 0 . 22 1 89675 06 P 0 0 

8 0.22367417390 00 

3 0.82751 031880 00 

4 0.82959492490 00 

5 0 . 8 3 0 6 5 7 35 1 5 P 0 0 

6 0 .33 1 1 896 0280 00 

7 0.2314570519D 00 

S 0.23 1596? 73 3D 00 

9 0.231 71 9375 ID 0 0 


OF DISCRETE DESCRIBING FUNCTION 
K2=K3=0,1 

PHI* 0.0 

AMPLITUDE 
0.1 18222 9409 D- 04 
0.1 0 6 5 9 7 1 9 6 0 D - 0 4 
0.1 198790559D-04 
0.1256521 1270-04 
0 . 1286281 146D-04 
0.I30 15727550- 04 
0 . 1 3 093743390^04 
0 .13133136390-04 
0.13152567790-04 
0.1 31 582 17280-04 


ND= 1 SAMPLING PERIOD* 2 .316950- 01 SEC AMPLITUDE* 1 .315730-05 N= 9 


Figure 2-7i 
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LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE) K£= 1 .OOOD-Ol » K.3= 1. OOOD-Ol 
GAMMA 3 1.38000D 0? N= 19 PHI= 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 .69000000000- 0 t 0 .68500000000-05 

1 0.6938008868D-01 0 . 71 04879736D-05 

3 0.7067198796D-01 0,706'i^*1387.i7D— 05 

3 0.7111818 1 7£D- 0 1 0 . 70828 063000-05 

4 0 . 713477636 ID- 0 1 0 .7015467044B-05 

5 0.7 1469586 OSD- 01 0 .7008501 293D- 05 

6 0 . 7 149620972D- 0 1 0 .7007741 175D-05 

ND= 1 SAMPLING PERIOB= 7 . 14990D-02SEC AMPLITUDE 3 7.00315D-06 N* 6 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE? K2= 1. 000D-01? K3= 1. 000D-01 
GAMMA 3 1 . 33000D 07 N= 17 PHI= 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0.71 49397 424D- 0 1 0 .7003143764D-05 

1 0. 735962 !££• ID- 01 0 . 770253 1 492D-05 

£ 0 .7612:67904 3D-01 0 . 7565776232D- 05 

3 0.7741 36 035 1 D- 0 1 0 .747791 0174D-05 

4 0.7S04933352D-01 0 .743175949SD-05 

5 0 .78367 1 532 1 D- 0 1 0 . 74 03489682D- 05 

6 0.7853062384D-01 0 .7396573734D-05 

7 0.736 038523 0 D - 0 1 0 .739448891 0D-05 

ND= 1 SAMPLING PERIOD= 7 .36051D-02SEC AMPLITUDE 3 7.3382 ID- 06 N= 7 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE? K£= I. OOOD-Ol) K3= 1. OOOD-Ol 
GAMMA 3 1 .38000D 07 N= 15 PHI- 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0.786 05 1 3935D- C 1 0 .73382 051 25D- 05 

1 0.813135972SD-01 0 .81 9 1990949D- 05 

£ 0 .8442897636D-01 0 .8070950868D- 05 

3 0 . 8598 06595SD- 0 1 0 .793725639 ID- 05 

4 0.8675165745D-01 0 .7942301 5 07D- 05 

5 0.3713578826D-01 0 . 7919630257D-05 

6 0.S733187734D-01 0.7907942364D-05 

7 0.8743345623D-01 0 .79 03856553®- 05 

8 0.3743533238D-01 0 .7900274383D-05 

ND= 1 SAMPLING PERIOD 3 3 .74330D-02SEC AMPLITUDE 3 7.90034D-06 N= 8 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE? K£- 1. OOOD-Ol? K3- 1. OOOD-Ol 
GAMMA 3 1 .38000D 07 N= 13 PHI 3 0.0 

ITERATION SAMPLING PEPIOD AMPLITUDE 

0 0.8743786585D-01 0 . 79 0 0333? 08 D- 05 

1 0.91 17515473D-01 0 .83236 11 04 OB- 05 

2 0.95 03524 084D- 01 0 . 87397921 SOD- 05 

3 0 .9708350587D-01 0 .866899023 00-05 

4 0 .97935968600-01 0 .86837309730-05 

5 0 .98464 0797 1 D- 0 1 0 .86 0333-1 827D- 05 

6 0.38705916420-01 0 .85978566870-05 

7 0.9S33613977D-01 0 .85927613510-05 

S 0.9S84468557D-01 0 .35942296 1 8D-05 

ND= 1 SAMPLING PERIOD 3 9 .334 79D- 02SEC AMPLITUDE 3 3.59222D-06 N= 3 


Figure 2-7j 
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LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE, K£= l.OOOP-Ol* K3* 1.000B-01 


GAMMA* t .’33 00 0D 07 N= 11 
ITERATION SAMPLING PERIOD 

0 0. 9334794 386D- 01 

1 0-1 044 74 74'? AD 00 

2 O . 1 095187279B 00 

3 0,11201431 45D 00 

4 0.1 1 32569 073P 00 

5 0,11 337264 15D 00 

6 0.11413114 15B 00 

7 0 . 1 1 4 3 45293 8 D 0 0 

3 0. 11 433 3333 6 D 00 

9 0 . 1 143991 £3?B 00 

10 0.11440 0643 OB 00 


PHI- 0.0 

AMPLITUDE 
0. 35928 19278D-05 
0 . 93403004500-05 
0 . 9 345574 3 68 D- 05 
0.93 0935 1 2 54 B— 05 
0 . 9533723 4 73 D - 05 
0 .957015535411-05 
0 .95333253 04D- 05 
0.955927351 OB- 05 
0.95345271 0211-05 
0 . 95 54 8 7923 OD- 05 
0.955331 0577D- 05 


NP* 1 SAMPLING PERIOD* 1 . 144 01 D- 01SEC AMPLITUDE* 9.55743D-06 


LST SYSTEM-NUMERICAL SOLUTION DF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE , K2= l.OOOB-Ol, K3= 1.000E-01 
6AMMA= 1.330 OOP 07 H= 9 PHI* 0.0 

AMPLITUDE 


ITERATION 

0 

1 

c! 

3 

4 


SAMPLING PERIOD 
0.1 14401 3031 D 00 
0.1 £41494U17D 00 
0 .130353773911 0 0 
0 . 13395 0732 OP 00 
0 . 1 35304 01 4£P 00 
0,1334233133 D 00 
0 . 1 33331 1 41 3D 00 
0.1 3704 133 19H 00 
0 . 1 371 7313 79P 00 


0 . 9 5574252 97 D - 0 5 
0. 1 037 754 3 92 D- 04 
0 . 1 091 3 0 1385E- 04 
0 . 1 09742 0953 D -04 
0 . 1 099 04 493SB- 04 
0.1 09974379 ID- 04 
0.1 100 1 3553 ID- 04 
0.11 00234373D-04 
0.11 002394 04 D- 04 


ND= 1 SAMPLING PERIOD® 1 .371 PSD- 01 SEC AMPLITUDE® 1.1004SD-0J 


LST SYSTEM-MUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 


TWO AXIS- 

EQUAL AMPLITUDE C 

:ase , 

K2= 

1.0 OOP- 01 p K3= 

GAMMA* 1 .33 00 OP 07 N* * 7 

ITERATION SAMPLING PERIOD 

PH I = 

0.0 

AMPLITUDE 

0 

0.137178132311 

0 0 

0. 

1 1 004 199 29 D- 04 

1 

0 . 153894 21 £6 D 

00 

0 . 

1 1 8339 0 4 33 D- 0 4 

2 

0 . 1 63 1 359 3 33 D 

0 0 

0, 

1 274324 753B- 04 

3 

0 . 17 04 183 526 D 

0 0 

0. 

13 07433534 B- 04 

4 

0.1 72 3332 833 D 

0 0 

0, 

1328051 O03B- 04 

5 

0.1 737443 075D 

0 0 

0. 

1 32 9 2 0 0 5 ? 7D - 04 

6 

0.1 74292797 ID 

0 0 

0. 

1332 379 123IJ-04 

7 

0.1 745694433B 

0 0 

0. 

1334713151 B- 04 

2: 

0.1 74 7839 923 D 

00 

0. 

13354 02053D- 04 

9 

0. 174558 1227B 

0 0 

0. 

1 33 7 0 783 75 D~ 04 

10 

0.17474 05935D 

00 

0. 


11 

0.1 7474 0795 IB 

no 

0. 

1 33325 1723B- 04 

NB* 1 SAMPLING PERIOD* 1 

.7474 

2D- 01 SEC AMPLITU 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE > K£= 1.000D-01> K3= 1 . 0 ODD- 01 


GAMMA* 1 .33 

00 OB 07 N= 

«r 

PH I = 0.0 

ITERATION 

SAMPLING period 

AMPLITUDE 

0 

0 . 1 747423SG5D 

0 0 

G.1336S17555B-04 

1 

0 .2495296658B 

0 0 

0 . 1 £7592445 4 D- 04 

o 

0.24670G9249D 

0 0 

Cl. 149433245 ID- 04 

3 

0 .251 08 192 OSD 

0 0 

0 . 1 595232373B- 04 

4 

0.2535581 ?25» 

00 

0 . 164974 7219D-04 

5 

0. 2548453 134D 

0 0 

0.1 67831 642 ID- 04 

3 

0 .25549592 1 SB 

00 

0 . 1 693 025549D— 04 

7 

0 . 2558226854 B 

00 

0 .1700500328D-04 

8 

0 .2559394 *: 1 6D 

00 

0.1704250683D-04 

9 

0 » 2 58 Go y 1 83 2 D 

00 

0 . 17060153660- 04 

10 

0. 853 033 333 4 B 

00 

0 . 17077 1 89G7D- 04 


ND= 1 SAMPLING PERIOD® 2 .56084D-01SEC AMPLITUDE® 1 .70779B-05 


N=10 


N= 8 


N=1 1 


N=1 0 


Figure 2- 7k 
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LST SYSTEM-NUMERICAL .SOLUTION DF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE > K2~ 1.000D-01* K3= 1.000B-01 
GAMMA= 1.33000P 0? N- 4 PHI- 9.000D0D 01 
ITERATION SAMPLING PERIOD AMPLITUDE 

0 0-50000 0 0 0 0 OD 00 0 - 1500000000B-04 / 

1 0.5154240450D 00 0 . 1377904361 D-04 

£ 0 * 5395'? 4 3 0 7 0 D 00 0 . 1446567295D-04 

3 0 .5480033734D 00 0 . 1474731 509B-04 

4 0 .5501232936D 00 0 . 14S6946298B-04 

5 0 .54983533 05D 00 0 . 1493377S98D-G4 

6 0.54361532561' 0 0 0 - 1496763755D-04 

7 0 .54955191 1 ID 0 0 0 . 1 4985359 1 4D- 04 

3 0.5495572657B 00 0. 1498 126 09 OD- 04 

ND= 1 SAMPLING PERIOD- 5 -49637D-0 1SEC AMPLITUDE* 1 .4991 3D- 05 N= 8 


LST SYSTEM-NUMERICAL SOLUTION DF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE > K£= 1.000D-01* K3- 1 .000D-01 
GAMMA* 1 .38 GOOD 07 N= 4 PHI- S..00000D 01 
ITERATIDN SAMPLING PERIOD AMPLITUDE 

0 0 . 5496365534 D 00 0 . 1499127478D-04 

2 5.49637D-01 -1.91174D-05 

2 5.50136B-01 -1 .95688D-05 

1 0 -4234 186796 D 00 0 .84545623 03D- 04 

2 4.23419D-01 -7 ,331 SOD- 05 

2 4.23S42D-01 -7.4581 ID- 05 

2 0 . 373 0698891 D 0 0 0 , 1 692428 1 45D- 04 

3 0.41 83434 144D 00 0. 1637625694D-04 

4 0.4312546179D 00 O . 1654907656D-04 

5 0. 43974 01 524 D 00 0 . 16651229 19D-04 

6 0.44 383 0 1 0 96 D 00 0 - 1 6 7 04 042 94 D - 04 

? 0 .44584331 95 D 00 0 . 167306484 7D- 04 - 

Q 0 .4463524 79 3D 00 0 . 16743908330-04 

9 0.44 739 0682 2D 00 0 . 167501 07 1 90-04 

10 0 -44752 056 OSD 00 0 . 1675 1 57086D- 04 

NB* 1 SAMPLING PERIOD- 4 .475 1 2D- 01 SEC AMPLITUDE* 1 .674S8D-05 N-10 


LST SYSTEM-NUMERICAL SOLUTION DF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE > K2= 1.000D-01? K3* 1.000D-01 
GAMMA* 1.3S000D 07 N= 4 PHI= 7.00000B 01 
ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 .4475120 178D 00 0 . 1 67488 1 54 OD- 04 

1 0.4216082457D 00 0 . 1861 170429D-Q4 

2 0.4051047902D 00 0.19000821440-04 

3 0 .3933394 0790 0 0 0 . 19 1331 8637 D- 04 

4 0.33676 04 154D 00 0 . 19281 80740D- 04 

5 0. 33254 075020 00 0 . 1933585459D-04 

6 0 * 33 0 1 1 77764 D 0 0 0 • 1 9363364 77D— U4 

7 0 .3787622 054 D 00 0 . 1937772905D-04 

8 0 .378G1602S8D 00 0 . I933492295D- 04 

9 0.3776 07655 OD 00 0 . 1938848 190D-04 

10 0 . 3773823844B 00 0 . 19390234 13D-04 

ND- 1 SAMPLING PERIOD* 3 .77255D-01SEC AMPLITUDE- 1.93911D-05- N»10 


Figure 2-8a. 
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LST SYSTEM 'NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO RX IS- EQUAL AMPLITUDE CASE* Kc= I.OQOB-Ot. K3= 1 .OOGIi-Ol 
GAMMA® 1 .33000D 07 N= 4 PHI= 6.00000D 01 
ITERATION SAMPLING PERIOD AMPLITUDE 

0 0.3?7’£54994£D 0 0 0. 19391 10S75D-04 

1 0.359G882449D 00 0 .21223 13759D- 04 

£ 0 .347443 1 1 5 1 D 00 0 .2199493076B-04 

3 0.33S9745293D 00 0 .££4639471 8D- 04 

4 0.333644S753D 00 0 .2275171224P-04 

5 0.33045S7295B 00 0 .2S92U5436B-04 

6 0.3236099425D 00 0 .2301675458D-04 

7 0 .3275734451 D 00 0 .2306894923D-04 

S 0.3270001735D 00 0 .2309677767D-04 

9 0.3266852982D 00 0 .231 U34342D-04 

ND= 1 SAMPLING PERIOD® 3 .2651 1D-01SEC' AMPLITUDE® 2.31189D-05 N= 9 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE f K£= 1.000D-01. KG- 1.000D-01 
GfiMMA® 1.3S000D 07 N= 4 PHI= 5.00000D 01 
ITERATION SAMPLING PERIOD AMPLITUDE 

0 0.32651 1468 OD 00 0 .231 1335352D-04 

1 0.31 34437 038D 00 0 .24969591 09D-04 

2 0 .3043395 095D 00 0 .261 56385 1 0D- 04 

3 0 .39 836 08 07 1 Ii 00 0,86? 86538 36D- 04 

4 0.29464 7 3 52 9 0 0 0.274 0556 0 6 OD - 0 4 

5 0 . 29844 0829 31? 0 0 0 . 8 7688 5 922 3D- 04 

6 0.291 171 7454 D 00 0 .8734392 157D- 04 

7 0 .2904579878D 00 0 .27937 12255D- 04 

8 0.29006 038 1 2 D 0 0 0.2793493 0400-04 

9 0 . 289836 0556 D 00 0 .83011339660-04 

MD- 1 SAMPLING PERIOD* 2 .89699B-01SEC AMPLITUDE* 2.S0285D-Q5 N= 9 


LST SYSTEM- NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE* K8= I.OOOD-Oli K3= 1 .OOOD-Ol 
GAMMA- 1,380001 07 N* 4 PHI* 4.00000D .01 
ITERATION SAMPLING PERIOD AMPLITUDE 

0 0. 889699 0149D 0 0 0 .23 023473 030- 04 

1 0 .879489 0678B 0 0 0 .3001 1331 03D-04 

2 0 . 27299 X 72471) 0 0 0 .31464 1 031 4D- 04 

3 0.86 6 9976 36 9 D 00 0 .3242627 184 D- 04 

4 Cl. 2666235 1 14D 0 0 0 . 330 15 1 346 1 D- 04 

5 0.2652524 683 D 00 0 . 3 3 3557 4996B- 04 

6 0 .26*4 4 7 7843 OD 00 0 . 3354573 933 D- 04 

7 0 . 264 048 1 1 88 D 00 0 . 3365 027538D-04 

8 0.2637896353D 00 0 .337 1 1 08549D- 04 

9 0 .263593384 2D 00 0 .33771443660-04 


10 

0. 26353723 1 OD 

0 0 

0 . 33734 29577D- 04 

1 1 

0 *263469986 0.D 

00 

0. 338835376811-04 

12 

0 . 8 6- c >• 6 4 cl b 6 9 D 

00 

0 . 3 43522 8752 D - 04 

13 

0 . 2633 1 0 1 399 D 

00 

0.33944665? ID- 04 

14 

0.263573 037 ID 

00 

0 .33704942450-04 

15 

0 . 8635730 069D 

00 

0 . 3368 04 0 1 7 1 If- 04 


HD= 1 SAMPLING PERIOD* 2 .63573D- 01SEC AMPLITUDE* 3.368040-05 N=15 


Figure 2-8b 
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LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TUD AXIS-EQUAL AMPLITUDE CASE ? K£* 1 -000D-01p K3* 1 -GOOD-Ol 


GAMMA* 1.33 

0 0 0 D 07 N= 

4 

PHI* 3 . 000 Cl Oil 01 

ITERATION 

SAMPLING PERIOD 

AMPLITUDE 

0 

0 .863573 0 069D 

00 

0. 3363 04 01 7 IB- 04 

1 

0. 25572 008 6 3D 

0 0 

0 .35682924230-04 

2 

D.2512315177D 

0 0 

0.37 0993692 OD- 04 

3 

0.248801 1337D 

00 

0. 37995265! 3D-04 

4 

0 .£4?429?53£D 

0 0 

0 .3S51949943D-04 

5 

0 .24 66 77 384 3D 

00 

0 .3381 190135B-04 

6 

0. £4626234090 

00 

0 . 38972568 1 7D— 04 

7 

0.246 01 7296 OD 

0 0 

0 . 39 07 0 37428 D - 04 

S 

0 .£45611493511 

0 0 

0 . 393 485 0 3 6 6 D - 04 

9 

0.245971 336 OD 

00 

0 .3364 233 696 D - 04 

1 0 

0.245331 0373D 

0 0 

0. 3391 65933 OD- 04 

1 1 

0 . £4 574 83 3 56 D 

0 0 

0 . 3 9 048 014 34 D - 04 

12 

0 .2460 09327 OD 

00 

0 .391 1 03742 ID- 04 

13 

0.24 57639 8 7 9B 

0 0 

0.3924392 727 D- 04 

14 

0 . 2453306847D 

0 0 

0.398139791 ID- 04 

15 

0.245643 042 ID 

0 0 

0.3930359796D-04 


HD- 1 SAMPLING PERIOD* £ .45643D-01SEC 'AMPLITUDE* 3.93086D-05 


LST SYSTEM-NUMERICAL SOLUTION OF III SORE TE DESCRIBING FUNCTION 


TWO AXIS-EQUAL AMPLITUDE CASE ? 
GAMMA* 1 -32 GOOD 07 N= 4 
ITERATION SAMPLING PERIOD 

0 0. £45643 042 ID 0 0 

1 0 ,£4034 15437D 00 

a 0. £3822626 14D 00 

3 0 ,23743300350 00 

4 0. 2371315275H 00 

5 0 .2369731 329 D 00 

6 0. 236859 16G9D 0 0 


K2* 1.00 CD- 01 > K3* 
PHI = 2.00000D 01 
AMPLITUDE 
0 . 3 93 0 859796 D - 04 
0.4 0398445 02D- 04 
0 . 4 0 686823 33 D- 04 
0 .4 06443532 0D- 04 
0 . 4 0 54 9 597 08 D- 0 4 
0-405 0 046427D- 04 
0.4052036897D-04 


1 . 0 0 OB- 01 


ND* 1 SAMPLING PERIOD* £ .36794D-01 SEC AMPLITUDE- 4.05412D-05 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE ? K2= 1 . 000 D' 01 , fc.3* 1 .000D-01 


GAMMA* 1 .38 

GOOD 07 N= 

4 

PHI* 1 . 0 00000 01 

ITERATION 

0 

JAMAL I Mb PERIOD 

AMPLITUDE 
0 . 1 674694 09 3D- 04 


1.; u 

1 

0. 32 1 7 32 7 9 85 D 

0 0 

0.2144 6 0841 6 D - 0*4 

2 

0.34793041 33 D 

0 0 

0 . 1476851 326 D- 04 

3 

0 .3560 04489 OD 

i j 0 

0.146951 3 7 06 D - 0 4 

4 

0.359343 00 03 D 

0 0 

0 . 1 4 68 4 87 728 D- 04 


HD* 1 SAMPLING PERIOD* 3.6 0729D- 0 1 SEC AMPLITUDE* 1 .46358P-05 


H=15 


N= 6 


NT* 4 


Figure 2-8c. 
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3. Study of Unequal -Ampl itude Oscillations in Two-Axis Coupled 
Nonlinear Systems 

3-1 . Introduction 

In the preceding two chapters conditions of self-sustained 
oscillations in the two-axis LST system have been studied by a 
numerical -iterative technique using the describing function. It 
is demonstrated that the numerical-iterative technique can be applied 
to the continuous-data as well as the discrete-data LST system. The 
difference between these two types of systems lies only in the 
stability equations, and these made no difference to the iterative 
algorithm. 

It is important to note that since the dynamics of the two 
axes are identical, it is assumed a priori that the amplitudes of 
oscillations, if there are any, are of equal amplitude. Certainly, 
all the solutions of the self-sustained oscillations obtained in the 
previous chapters are sinusoidal signals with equal amplitude in 
both axes. However, in this chapter we shall show that, in general, 
it is possible to have unequal magnitudes of oscillations in the 
coupled system even if the dynamics of the two axes are identical. 
Furthermore, when the amplitudes of oscillations are equal, the 
phase can still be 180 degrees apart. 

Since we have reasoned that the numerical-iterative technique 
does not distinguish a continuous-data system from a discrete-data 
system, only models of the former version are discussed here. 
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3-2. Conditions of Unequal-Amplitude Oscillations in The Continuous- 

Data Two-Axis Coupled LST System 

Figure 3-1 shows the signal flow graph of the continuous-data 
two-axis model of the LST system. The inputs to the two nonlinearities 
are designated as E-j (jto) and E 2 (jw), respectively. Notice that these 
two variables are complex quantities in general, and each represents 

o 

a magnitude and phase in the sinusoidal steady state. Since the 
describing function of the CMG nonlinearity is a function of the 
input amplitude and frequency. Fig. 3-1 shows that the describing 
functions which correspond to E-j(jw) and E 2 (jw) are N(E^,w) and 
N(E 2 ,u 3), respectively. For simplicity, the following notations are 
used: 

E-| ( Joj) = E-j 

( Jco) — £2 

N(E r w) = N 1 

N(E 2 ,w) = N 2 


The transfer functions in Fig. 3-1 are defined as 


G, - 


(K o + K lS )Kj 


1 K p s + Kj 


(3-1) 


K s + K t 
G„ - 


(3-2) 


G~ = 


3 J G s 


(3-3) 
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“I 



Figure 3-1. The simplified continuous-data two-axis 
LST system model. 
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G 


4 



+ K 3 s + K 2 


(3-4) 


The determinant (characteristic equation) of the coupled system of 
Fig. 3-1 is written 

A = 1 + G A + G b1 N 1 + G b2 N 2 + G c N 1 N 2 = 0 (3-5) 


The last expression is also known as the stability equation of the 
system. 

The linear transfer functions in Eq. (3-5) are 


G 


A 


2G 2 G 3 


+ EG^G^G^G^H - K^G^ + GgG^ 
+ 2G 1 G 2 G3G 4 H - 2G 2 G 3 G^K^ - 


+ (G ] G 2 G 3 G 4 H) i 


G 2 G 3 G 4 K 2 


(3-6) 


G 



G 3 G 2 G 3 , G 1 G 2 G 3 G 4 H 

+ 4 * 

S S S 


2 ? ? ? ? 
G 3 G 4 K 2 G 2 G 3 G 4 K 2 


(3-7) 


G 


C 



r 2 rV 

G 3 G 4 K 2 


(3-8) 


Equation (3-5) represents a necessary condition for self- 
sustained oscillations in the system of Fig. 3-1. However, in general, 
if the amplitudes of E-j and are not equal, Eq. (3-5) represents a 

set of two nonlinear equations with four unknowns in E-j , E^, the 
frequency w, and 4), the phase between E-j and E^. Therefore, in 
addition to the two equations which are obtained from the real and 
imaginary parts of Eq. (3-5), a third complex equation must come from 



the relation between E-j and Eg which is derived as follows. 

The following equations are written directly from Fig. 3-1: 


E i ■ k 


x 4 = G 3 X 3 


x 3 = ~ N 1 E 1 

+ GgX 

x 2 = ' x 4 - 

S l e v 

0 v = G 4 X 5 


x 5 = Hx 4 - 

K 2 0 H 


c _ 1 
E 2 s y 4 

y 4 = G 3 y 3 

y 3 = “ N 2 E 2 + G 2 y 2 

y 2 = ” y 4 " G 1 G H 

9 H = G 4 y 5 

y 5 = Hy 4 " K 2 0 v 

From these equations it is observed that 


(3-9) 


(3-10) 



n 2 g 3 

1 + + G 2 6 3 h 

1 + 4” + G 2 G 3 


(3-11) 


Eliminating E 1 , Eg, e v> e^, x 4 and y 4 from Eqs. (3-9) and (3-10), we have 
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N 1 G 3 

x 3 " s X 3 + G 2 X 2 
x 2 = " G 3 x 3 " G 1 G 4 X 5 
x 5 = HG 3 X 3 “ K 2 G 4 y 5 

n 2 g 3 

y 3 = s~ y 3 + G 2 y 2 
y 2 = " G 3 y 3 " G l G 4 y 5 ■ 
y 5 = HG 3 y 3 " K 2 G 4 X 5 

These equations are solved to yield 

x 5 1 + + B + C 

= -k 2 g 4 (i + a 2 + b) 


(3-12) 


(3-13) 


(3-14) 


or 

x 5 _ -k 2 g 4 (i + a 1 + B) 
y 5 1 + A-j + B + C 

where 


A 


A 


^3 

1 s 


N 2 G 3 
2 “ s 


B = G 2 G 3 


(3-15) 


(3-16) 

(3-17) 

(3-18) 


C = G^G^H 


(3-19) 
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Thus , from Eq. (3-11 ) , 

Ei (1 + Ag + B)Xg 1 + Ag + B + C 
E^ = (1 + A 1 + B)y 5 = -K 2 G 4 (1 + A ] + B) 

or 


E 1 “ K 2 G 4^ + A 2 + B) 
Eg 1 + A-j + B + C 


(3-20) 


(3-21) 


It is interesting to point out that the two equations for in 

Eqs. (3-20) and (3-21) should be equal to each other. Thus, 


E-j -KgG^O + Ag + B) 1 + Ag + B + C 
Eg 1+A^+B + C - KgG^ ( 1 + A-j + B) 


(3-22) 


When the last equation is expanded, it can be shown that it is 
identical to the condition of A = 0 given in Eq. (3-5). 
Multiplying Eqs. (3-20) and (3-21) yields 


'E (1 + A 2 + B)(l + A 2 + B + C) 

= (1 + A 1 + B) (1 + A 1 + B + C) 

When |E-|| = | E^ | » N-j = Ng, and A-j = Ag; Eq. (2-23) gives 



(3-23) 


(3-24) 


Equation (3-24) implies that E-| and Eg can be of the same 
phase or they can be 180 degrees out of phase. Instead, when the 
numerical -iterative technique is applied to the | | = jEgj = |E| 
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situation, it is found that the case, E-j = - E ^ is true. 

Returning to the stability problem, Eqs. (3-5) and (3-23) form 
the necessary conditions of self-sustained oscillations for the 
general case of the system in Fig. 3-1. 

Since Eqs. (3-5) and (3-23) are independent mathematical 
relations, it is conceivable that in general there may be E-| f E 2 
solutions which may satisfy these equations simultaneously, thus 
signifying the possibility of self-sustained oscillations with un- 
equal amplitudes in the two axes of the coupled system. 

Now since there are three unknown parameters in the magnitudes 
of E-j and E^, and w, it is difficult to arrive at a sustained-oscillation 
solution via the numerical-iterative method for the complex LST system. 

For the given system with the specific nonlinearity, even if a 
sustained oscillation exists, one has to be able to give a proper set 
of initial guesses for E-j , E 2 » 4>, and w in order to arrive at the 
convergent solution. This has been proven to be quite difficult for 
the LST system. In order to gain further insight to the problem, . 
two systems of lesser complexity are considered in the following sections. 
Through these simpler systems we are able to draw general conclusions 
on the possibility of having unequal -ampl itude oscillations with the 
CMG nonlinearity. 

3-3. A Simple Two-Axis Coupled Nonlinear System 

Consider the simple two-axis coupled system shown in Fig. 3-2. 

The characteristic equation of the system is 

A - 1 - k 2 G 2 N 1 N 2 


(3-25) 
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A simple two-axis coupled 
nonlinear system. 


The input signals of the nonlinear elements are written from 


Fig. 3-2. 

E-j = kx 2 = N 2 GkE 2 (3-26) 

Thus, 


-A- = N 9 Gk (3-27) 

E 2 2 


Also, 


E 2 = kx 1 = N-jGkE-j 


(3-28) 


from which we have 


E 2 I^Gk 

Equations (3-27) and (3-29) are multipled to give 


(3-29) 



(3-30) 


In order for the system to exhibit self-sustained oscillations, it is 
necessary that A = 0 be satisfied, along with any one of Eqs. (3-27), 
(3-29), and (3-30). These provide two complex or four real equations 
with four unknowns in |E^|, |E 2 |, to, and <J>. 

In general, if appropriate initial guesses are available for 
|E-]{, |E 2 j, to, and <j), the numerical-iterative method will seek out 
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the solution to the system of nonlinear equations represented by 
the equations discussed above. This is useful for the case when 
|Ej| = jE^I = E, and <j> = 0, since the initial guess for these 
solutions is available from the single-axis system solutions. 

In the case when | ] f | [ and <j> f 0, a general procedure 

is needed to solve the nonlinear equations. 

One method of solution is to span the to, |£^|, and | ] space 
and plot A (or A - 1) or Eq. (3-25) in magnitude versus phase co- 
ordinates, and seek out the A = 0 solutions. These solutions then 
become candidates for possible solutions of Eq. (3-27), (3-29) or 
(3-30), with the variable parameter 4> free. 

An alternative is to solve Eq. (3-27), (3-29) or (3-30) 
first and then use the solution or solutions to check with Eq. 
(3-25). Along this direction, a necessary condition for the 
existence of limit cycles is given by Eq. (3-30) which is inter- 
preted as follows. 

Equation (3-30) is written 


2 



E 1 


N 2 

E l 


n i | 


Since = N(E 2 ) and N-j = N(E-j), Eq. (3-31) will be satisfied 

2 

if E and K/N(E), where K is an arbitrary constant, have the same 
magnitude for two different values of E. This means that if the 
curves of E and K/N(E) are plotted with respect to E, there must 
be two intersections between these two curves in order to satisfy 
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Eq. (3-31). Figure 3-3 shows the plots of E 2 and |1/N| in db versus 
E in db, with N being the describing function of the CMG nonlinearity 
with y = 1.38 x 10 6 . It is apparent that these two curves intersect 
at at most one point for any K. Therefore, we may conclude that 
with the CMG nonlinearity in the system configuration of Fig. 3-2, 
self-sustained oscillations of unequal amplitudes cannot occur with 
any transfer function G. 

Let us define 

f * |N 


Also, let 

9 = ' t °g 1 o h (3-34) 


(3-32) 

(3-33) 


x = log 10 E (3-35) 

Then, the slope of the curve for j!/N| in Fig. 3-3 is given 
by 


dg _ 1 dh df dE _ r . df 
dx h df dE dx " tn dE 


(3-36) 


or 


4a = 

dx 


-Eh 


d|N 

dE 


(3-37) 





Ill 


With 


f = 


(3-38) 


Eq. (3-36) gives 


da- 2 

dx 


(3-39) 


Thus, the necessary condition for Eq. (3-31) to be satisfied is that 
the slope of the 1 og-j q j 1 /N | versus 1 o9-|qE curve exceeds the value 2 at 
at least one point. This will ensure two intersections between the 
E^ and |1/N| curves. 


3-4. A Simplified Multiple-Loop Two-Axis Coupled Nonlinear System 

Figure 3-4 shows the signal flow graph of a two-axis coupled 
system which is simpler than the LST system and yet still retains 
the essential characteristics of the system. 

The transfer functions of the system are defined as 


Go = ^ (3-40) 

s 

H = s + 10" 3 


These transfer functions and parameters are selected so that 
the single-axis system has two limit cycles, one stable and one un- 
stable, with y = 1-38 x 10^ and T q = 0.1 for the CMG nonlinearity. 
The determinant of the coupled system is 
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“I 



"I 


Figure 3-4. 


+ N 1 N 2 G 1 2 (1 - k 2 g 2 ) 


(3-41) 


The relationships between E-j and E 2 are determined from the following 
equations: 


x 2 

y l 

y 2 


HG-|X-| - K-jG 2 y 2 

”^2^2 - ^2^2 
HGiy i - k.jG 2 x 2 


Thus, 


(3-42) 


E-j X-| 1 + N^G-j + HG-jG^ 

= ^7 = -K 1 G 2 C1 + N-j G-j ) 


(3-43) 


or 


E 1 -K-jG 2 (1 + N 2 G 1 ) 

= 1 + N-jG-j + HG^ 


(3-44) 


Again, if we equate Eq. (3-43) to Eq. (3-44), we would get A as in 
Eq. (3-41). 

From Eqs. (3-43) and (3-44) we have 

e } ) 2 (1 + n 2 g 1 )(i + n 2 g 1 + hg 1 g 2 ) 

E^ = (1 + N-j G-, )(1 + Mi G-j + HG 1 G 2 ) 


(3-45) 
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The numerical-iterative method was first used to solve for the 
conditions of limit cycles for | E-j J = |E 2 | given by Eqs. (3-41) and 
one of the equations in Eqs. (3-43), (3-44), and (3-45). It turned 
out that when |E-j) = |E 2 ), 4> = 180°; that is, the two sinusoidal 
signals at the inputs of N-j and N 2 are 180 degrees out of phase. 

Using appropriate initial guesses for |E^| = |E 2 |, to, and with <f> = 

180°, the numerical-iterative method yielded limit cycle solutions 
shown in Fig. 3-5 for various degrees of coupling coefficient K-j . 

Figure 3-6 shows the numerical iterations for the various values of 

K r 

A graphical procedure was carried out for the solution of limit 
cycles when E-j f E 2> since the numerical-iterative method yielded no 
convergent solutions. 

Figure 3-7 shows the plots of A of Eq. (3-41) for various com- 
binations of [E-|| and |E 2 |. For limit cycles to exist the plot of A 
must go through the (0 db, -180°) critical point. Figure 3-7 shows that 

— 'X 

the |E-]| = |E 2 | = 0.9163 x 10~ J locus does go through the critical point. 
However, most of the |E-|| f |E 2 | trajectories are far from intersecting 
the critical point. Although, it is difficult to span the |E-j|, (E^j 
space completely, these plots do give indication of possible existence 
of unequal-amplitude oscillations. 

The other necessary condition that needs to be checked is given 
by Eq. (3-45). Let 

f = (1 + Nfi 1 + HG 1 G 2 )(1 + NG-j ) (3-46) 

9 

Then, for a limit cycle solution, it is necessary that the E versus 



E curve intersects the K/f (K = arbitrary constant) versus E curve 
at two points. Figure 3-8 shows the trajectories of these curves 
in the decibel coordinates. Since the curves for 201og^K/f is raised 
or lowered for different values of K, the following possible solutions 
are indicated from Fig. 3-8: 


10 

Range of ]E-| | and 

10' 3 

10~ 5 - 10" 4 

10" 2 

10‘ 3 - 10" 5 

10" 1 

-5 -1 

10 3 - 10 1 

1 

10" 4 - 10" 2 

10 

10" 4 - 10" 3 


However, checking with the A curves in Fig. 3-6, they show 
that A cannot be zero in the ranges of amplitudes indicated above. 

The conclusion of this study is that although theoretically 
sustained oscillations with unequal amplitudes can occur in a two- 
axis coupled system with identical dynamics, the nature of the CMG 
nonlinearity most likely precludes this possibility. Another 
interesting discovery is that while a great majority of the limit 
cycles are of equal amplitude in the two axis, the phases are 180 
degrees apart. 




Figure 3-5. 
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0 . 1 792237 D 

03 

4 

0 .791492 0D-02 

0. 1 04 08 0511-01 

0 . 1 0593 04 D~ 0 1 

0 . 1 797632D 

03 


0 . 8 1 6£4 1 3D- 02 

0 . 1 034 69 OD- 01 

0.1 036335D- 0 1 

0 . 1 799723D 

03 

« 6 

0 .3 19313 3D -02 

U . 1 03 :j? boD- 0 1 

0. 1 03377311-01 

0 . 1 8 0 001 5D 

03 

i 

0.8 19358811-02 

0. 1 0337461*- 01 

0 . 1 033746D- 0 1 

0 . 1 $00 01911 

03 

IND= 1 FREQUENCY* 8.19858D- 

03R AD-SEC El= 

1 .0337511-05 E2= 

1 . 0337511- 02 



THETA= 1 .S0005IJ 05 H I T = 


Figure 3- 6b 
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TEST SYSTEM-NUMERICAL SOLUTION OF PELTfiCU »Ei »E£ »THETfi>»Q . 

GAMMA- 1 .SSOOOIi 06 TD= 1.000000-01 hSM= 1. 000000-03 Kl= 8 . 0G000D-04 


ITERATION FREQUENCY El E£ 

0 0.81 985850- 0£ 0. 1 0337460-01 0.1 0337460-01 

1 '0.34824850-02 0. 18332740-01 0 .218081 70-01 

£ 0 .31219260-02 0 .£7207 1 00-0 1 0.3050192D-01 

3 0:36945730- 0£ 0. £4404670-01 0. £7379730- 01 

4 0 . 4333 1 091'- 03 0 . £££43 OSD- 0 1 0 .£3343150- 01 

5 0 .49 03.5050- 0£ O. £0437650-01 0 .£1 075480-01 

6 0 . 533599 9 D- 02 0 . 1 943£S 1 0-0 1 0 . 1 963£90D-01 

7 0. 5546392B-0S 0 . 1399701 D- 01 0.19031360-01 

3 0 .559 02 060-02 0 . 13913630-01 0 . 1 891 4 1 80- 01 

9 0.55914650-02 0. 33911210-01 0.13931130-01 

10 0.55914500-02 0.13911230-01 0 . 1891 1230-01 


PHASE-Et 
0 . 13000190 03 
0.17535310 03 
0.1 7539420 03 
0.17718120 03 
0.17347430 03 
0.17922090 03 
0.17972690 03 
0.1799537H 03 
0.1300006D 03 
0.13000190 03 
0.1300019D 03 


IN0= 1 • FREQUENCY® 5 .59145D-O3RA0.'SEC El- 1.891120-02 E2= 1.391120-02 

THETA= 1.80002D 02 NIT- 10 


TEST SYSTEM-NUMERICAL .SOLUTION OF BELTACW jEI >E£ »THETR>=G . 

GRMMA= 1.330000 06 TD= 1.000000-01 ASM- 1.000000-03 Kl= S.50000D-04 


ITERATION 

FREQUENCY 

El 

£2 

PHASE-E1 


0 

Ci .559 145 ID- 02 

0. 1391 1230-01 

. 0.1891 123Ii- 0 1 

0 , 1 8 0 Ci 0 1 9D 

03 

1 

n ,456 025 1H’ 02 

0.23995170-01 

0.2316S68I'-01 

0 . 1 797 862 Ei 

03 

c. 

0 .49612990-02 

0.2368 9920- 01 

0 -2378736 D- 01 

0 . 1 7992970 

03 

o 

Ci ,47 V;: 07511- 02 

0 .23546090-01 

0.2354633D-01 

u , i 

03 

4 

0 .47154530-02 

0 .23533090-01 

0 .£35380211-0 1 

0 . 1 S 0 Ci 0 1 9 D 

03 

5 

0.471 542.1 D'0£ 

0.23538 130- 01 

0.235381SD-01 

0,180001 9D 

03 


IMP- 1 FREQUENCY- 4.71 5420- 03R:hB.'SEC El- £.353320-02 E£= 2.353320-02 

THETA- i .3 00020 02 NIT- 5 


TEST SYSTEM 
GAMMA = 1 .38 

-NUMERICAL SDLUTIDN DF DELTAS U ? E 1 
Ci 0 Ci Ii 09 T D = 1 . 0 Ci 0 C» 0 £- 0 1 A S M = 

9 E2 ? THETA > = Q - 
1 .000000-03 

Kl= 9 . 0 00 000- 

04 

ITERATION 

0 

FREQUENCY 

0.47154220-02 

El 

0. £3533180- 01 

E2 

0 . 83533 1 SB- 0 1 

PHASE-E1 
0 . 130 00190 

03 

1 

0 .34350920-02 

0 . 3 1 4 2 7990- 0 1 

0 . 31 075950- 0 1 

0 .1791 95 1 D 

03 

c 

0 .35329830-02 

n ,32546880-01 

0.383 05390-01 

0.1 7977280 

03 

■ 7 ; 

0 . 3659 1 7 00- 02 

0.3172 7350-01 

0.31 733340-01 

0 . 1 7996600 

0,3 

4 

0 . 36343 1 00-02 

0 .31559890-01 

0.31559950-01 

0 . 1 3 0 0 0 1 0D. 

03 

5 

0 . 36354230- 02 

0 .31556280-01 

0 .31558130-01 

0 . 1 8 0 Ci 0 1 90 

03 

6 

0 . 36854 1 00-02 

0.3 155 93 in- 01 

0.31559310-01 

0 . 1 8000190 

03 


I N0- 1 FREQUENCY- 3 .6354 1B-03RAD-'3EC El- 3.155630-02 E£= 3.155630-02 

THETA- 1.300020 02 NIT- 6 


TEST SYSTEM-NUMERICAL SDLUTIDN CF 0ELT A < U ? E 1 • E2 <• THETA > = 0 . 

GAMMA- 1.380000 06 TD= 1.000000-01 ASM- 1.000000-03 K1- 9.500000-04 


ITERATION 

FREQUENCY 

El 


E2 

PHASE-E1 


0 

0.368541 OD- 02 

0.31556300- 

ill 

0.3 15593 OP- 01 

0 . 1 8 0 0 0 1 9 D 

03 

1 

0 . 19769880-02 

0 .41593360- 

01 

n . 55529c.' 90- U 1 

IJ • 1 <■ cl 6* c 2 8 D 

03 

2 

0 . 16574160-02 

0 . 6* 5 9 6 9 6 6' D — 

01 

0.70 019210-01 

0 . 1 77 06*0 00 

03 

3 

0 . 18003270-02 

0 .6581 8770- 

0 1 

0 . 9 9 9 9 8 5 3 D — 0 1 

0 . 1 790 1 01 O 

03 

4 

U . 2 02 01 040-02 

0.59547120- 

01 

0 .58787540-01 

0 . 1 7952880 

03 

iC 

U .21967 040-02 

0,549 14390- 

0 1 

0 .5484 0-^50-01 

0 . 17970690 

03 

6 

0 » c 3 2 W c. O — 0 cl 

0.52179190- 

01 

0 .53190990-01 

0 . 1 79889 1 O 

03 

7 

0.23815360-02 

0.51 073 080- 

01 

0 .5 1 07 0250- 01 

0 . 1 7998390 

II 8 

o 

0 .239 1 5580-02 

0 .50893350- 

01 

0 .5 08Y26UD-I.it 

0. 18000190 

03 

9 

0.23914990-02 

0 .5 0894 760- 

01 

O.5O334770-OI 

0 . 1 3 Ci 0 0 1 9 D 

03 


I HD- 1 FREQUENCY- £ .391500- 03* A D-- SEC El- 5.039470-02 £2= 5.089470-02 

THETA- 1.800020 02 NIT- 9 


Figure 3-6c 



(n,tn) represents 

|,|e 2 |) - uo~‘\io~" 



DEGREES 


Figure 3-7- 
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DISTRIBUTION LIST 


Final Report 1-75 NAS8-29853 

AS21D 5 
ED12/S/Seltzer 10 
EDll/G/Nurre 1 
EM34 1 
AT01 1 



